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EFFECT OF CERTAIN TYPICAL NONLINEARITIES ON THE 
ADJUSTMENT OF A MECHANICAL PILOT 


M. E. Salukavadze 


(Tbilisi) 


The power balance method was used for studying the effect on the mechanical 
pilot adjustment of control nonlinearities, due to a restricted aileron turning angle in 
list stability tests, of the considerable drive nonlinearities and the mechanical pilot 


mixer relay quasi-nonlinearities in longitudinal plane stability tests. Similar problems 
were studied in[1, 2}. 


Effect of Restricting the Aileron Turning Angle on the Mechanical Pilot Adjustment. 


The free movement of an airplane, at a constant speed, with a list can be approximately expressed in the 
form 


Jt = + 
where J, is the airplane's moment of inertia in kG -msec*, Mg st is the natural damping factor in 


kG * msec (moX < 0 is the degree of banking damping, y => V* is the aerodynamic thrust in kG m?, p is the 
v is the speed of the plane in m*sec™’, 


is the wing area in is the wing span in m), is the aileron efficiency in kG-m < 0 


is the degree of aileron efficiency), y is the angle of list (y » 0 when listing on the right wing), 5a is the aileron 
deflection angle (6, » 0 when the starboard aileron is turned*down). 


constant time delay T,=const > 0 and the time constant T». 


Mechanical pilot equations; summing-telay equation 
q = Ky %, (2) 


where Ky is the angle transfer ratio, Kj is the angular velocity transfer ratio in seconds; the equation of the ampli- 
fier and power drive 


+3=q(t—T)); (3) 


the kinematic equation of the aileron drive 
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where A is the aileron turn-restricting angle. 


Let us now find in the plane of two mechanical pilot parameters (Ky and Ky) the limits of the aeroplane 
stability. 


In analyzing the airplane periodic oscillations (if they exist), it is possible with an accuracy sufficient for 
practial purposes to deal with the fundamental only, i.e., to consider that 


= 1 8in wt + 7, cos of, (5) 


where y, and 7, are Fourier coefficients of the first (fundamental) harmonic with an angular frequency of w. This 
assumption is justified in certain cases by flying experiments. 


If (5) is fulfilled, we have 
= o,sin w(t — ty). 


Since the system is self-contained, it is possible to assume that tg=0 and 


= d98in wt. (6) 


By using equations (2) and (3) it is possible to express the Fourier coefficients of the function of y by means 
of amplitude og: 


A, = cos wT, — oT, sin oT ;, A, = oT, cos wT, + sin wT’, 
rR K+ (K;, 


The amplitude of the fundamental 7 is expressed in terms of its Fourier coefficients according to formula 
to= To). (8) 


Let us determine the relation Y, and y, to A by means of the power balance method. For this purpose let 
us multiply both sides of equation (1) by jand integrate them with respect to t between the limits of 0 and 


T (T=4F, the period of oscillations). We obtained the well-known conservation of energy equation 
T T 


+ Mia = 0. 


The first term represents the work of the dissipation forces for the period, the second represents the work of 
the aerodynamic flow forces. With the assumption made in (5) the last equation takes the form: 


+- (3,c0s ot — a} = 0. 


Aat 
s<—A, 
where 
0 


The integrals of the equation are Fourier coefficients of an odd function of a= 6a(a). Hence the first one 


disappears and the second one is easily calculable: 


where 


tain 


T 
(9) 
Substituting the integrals in the last equation by their values just obtained we have 
(10) 


Muitiplying both parts of equation (1) by coswt and integrating with respect to t between 0 and T we ob- 


+ 0. 


Inserting the value of ¥, from (10) we have: 


Je 4% 
(11) 
x 


From relationships (7), (8), (10) and (11) we have system of equations in terms of Ky and Kq? 


J, 
+ = — [1 + 


(12) 
AK, — AwK. = [1 + xw. 
Since the determinant of the system 
elt 
is always distinct from zero at w > 0, it follows from system (12) that 
= + cos oT: + (M—o*T,) sinew? 
(13) 


K, = sin —(M — 
mu 


where M= ~ = and is determined from formula (9). 


The parameter og @ A. The equality sign between oy = A corresponds to a linear case when X= 1, 
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Relationship (13) determines the boundaries of the stability region in plane{Ky, Ky} corresponding to the 
oscillation amplitudes 


Yo = %h, (14) 
where 
h = {{1 + / + 


_ From equation (13) it is possible to calculate Ky and K} corresponding to various values of w > 0, and to 
plot a curve on plane {k, Kj}. For the linear case, i.e., when *= 1, the curve has approximately the shape 
shown in Figure 1. 


The shaded region (Figure 1) represents the part over which the airplane list is stabilized in the plane of 
two mechanical pilot parameters Ky and K}, 


From the above the following conclusions can be drawn: 


1) The stable region of the airplane in the plane {K,, Ky} is restricted by the inequality K, > 0, Ky> 0 
and the parameter defined curve r(w >0). Point {K,, Ky} of curve [specified by two definite numbers w > 0 
and 0 = A, corresponds to a stable state flying conditions with list oscillations of an amplitude yp and an angular 
velocity 


The latter should be understood as follows: each point on curve I corresponds to a stable state oscillatory 
condition of flying with an amplitude of ¥, and a frequency of w; each point inside the curve in region G at 
Ky >0, Ky > 9 can only correspond to a condition asymptotically approaching zero. The system's behavior out- 
side curve [' can be obtained by studying the effect of higher harmonics. 


Fig. 1. Region of stabil- Fig. 2. Stability boundaries with 
ity in the plane (Ky, Ky). different T,. 


2) If the region of stability of the system with A= dy is known, i.e., if the curve Ig) = 4 is known, the curve 
Tg, >A can be obtained from the first one by a geometric extension transformation with coefficient x: 


= 
3) It will be seen from relation (13) that the region of stability varies inversely with T, (Figure 2). 


With T,=0, quantity Ky is not restricted, but with an arbitrary small T; > 0, quantity Ky becomes re~ 
stricted, 


Logitudinal Stability of an Automatically Controlled Airplane . 


In analyzing the longitudinal stability of an airplane with a mechanical pilot, let us examine the consider- 
able nonlinearity of the power drive and the quasi-nonlinearity of the mixer relay. 


If the airplane speed V fluctuations are not taken into account, the longitudinal movement equation can 
be written as: 


Jb = + Mia + Mia + 


(15) 


where ® is the angle of pitch, 8p is the angle of elevator rotation from the neutral position, c= & ct is the 
deviation of the angle of attack & from its stable value a», determined from the conditionsof a stable state flight 


taney Cy) / Gay Gy resistance coefficient, Cy is the lifting force coefficient (a function of «), Jy is 
the aeroplane's moment of inertia about the transverse axis Oz from point O in the direction of the right wing in 
KG + msec?, = kGimsce (me is the degree of pitch damping), md y KG ts 
the degree of Mg = mg ySb, kG*m (mf is the degree of longitudinal static stability), Me = kG*m 
end B is the elevator efficiency). In the second equation (15) 


Fy 


where m is the mass of the airplane in kG-sec’+m™’, ba is the wing chord in m, CY =.9Cy/dc: is the angle of 
slope tangent of curve Cy at point a, and P is the tractive force of the engines in kG, It ts assumed that force P 
does not produce a moment about axis Oz. 


Let us examine the mechanical pilot with- 


8 7 kgd out feedback, whose block schematic is shown in 
Figure 3 in the same manner as for directed act~- 


According to the deviation of the airplane 


+, from a set course |9= $= $= 0} with respect to the 
155 angle of pitch 9, the speed § and acceleration 5, 
ye the sensitive elements unit transmits to the mixer 
Fig. 3. Block schematic of the mechanical pilot, 
vertical gyro, K§> from the precision or damped 
1) Precision gyroscope; 2) mixer relay; 3) amplifier; 
4) servomotor; 5) elevator; 6) angular acceleration Bysoscope ng 


meter; 7 vertical gyro. ay 


Constants Kg, Kg, and Kg are known as 
transfer ratios, and serve in this circuit as adjust- 
ing parameters. 


Mixer relay 2 sums up the signals algebraically. If $= 9 (t) the transmitted signal, due to delay T, in the 
relay, will be: 


q=K,9(t +7,) + K,d(t—7,) + Ky (¢—T7)). 


(16) 
The mixer relay signal g is amplified in amplifier 8, which has a time constant of Ty: 
(17) 
The speed of elevator resetting |is determined by the servomotor volocty characteristic: 
= (18) 


where f(a) is a continuous rising odd function (Figure 4). 


mV 
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In order to investigate the effect of the servomotor nonlinearity, let us study the velocity characteristic 
frequency spectrum assuming that o= 0, sint. It is assumed that in the interval |-7,+1| the velocity characteris- 
tic satisfies Dirichlet's system of sufficient conditions: 1) that the interval |-1, +] can be divided into a finite 
number of partial intervals in each of which the function is monotonic, 2) that all the points of discontinuity are 
regular, 

Hence due to the function being odd in the interval |-1,+1| the velocity characteristic can be represented 
as a Fourier series: 


(eosin t) = C,sint + C,sin 2t¢+...+ C,sint, 


where Fourier series coefficients C,,, Cyy..., Cy; 
dt 


are functions of the amplitudes 09. 
The amplitude~phase frequency characteristic of the seryomotor is represented by the following expression 


+n 
© = = tat. (19) 


It is said that the servomotor has a “filtering” property if | Cy4,| << || C,| (n=1,2,...), i.e., if it is possible 
to neglect the Fourier coefficients C,, Cy ...,Cn as compared with C;. Then we obtain approximately: 


For instance, the linear function f(0)=8o has “filtering properties, 
10) The amplitude-phase frequency characteristic (APFC) for it will be 
= 8 (6 is a constant), 


Let us assume that with certain values of Kg, Kg, and K§ an 
airplane with a mechanical pilot is oscillating periodically with re- 
spect to its angle of pitch 9. 


It is known that the airplane possesses sharply defined “filtering” 
properties, and the Fourier coefficients of the higher harmonics are in 
id their absolute values negligibly small as compared with those of the 
Fig. 4. Servomotor velocity first harmonic 9, and 9,, Hence we can assume approximately that: 
characteristic, 


§~9,sinet +%,coswt. (20) 


The amplitude components 9, and 9, can be expressed by means of the mechanical pilot parameters. Let 
us assume that the initial time of observation t, is chosen so that with condition (20) quantity o varies according 
to the law: 


= o98in wt. (21) 
Then from equations (16) and (17) we obtain: 
9, = [A(K, — Kgo*) + % =(B(K,— Kyo") — 40K, (an 


where 
A B= sinw7, + 07, cosw7,, 


= (K,— + (K,o)*. 


An expression for the amplitude of the fundamental will now be required re Se, 
%=+ 1+ . (23) 


Let us write the power equation for the airplane: 


K + F(a, a) — M°*8,$ = 0. (24) 


It is derived from the system of equations(15). The first two terms refer to the power (kinetic, potential, and 
dissipated), the third term refers to the power received from the aerodynamic stream; 


K = + 
F (a; a) = + cree + 
where 
ay = —(M*+ M*), = — M*) Fy — 
Qe, = — (MP Fy + M*) F,. 


On the basis of (24) it is possible to write down an energy equation considering that with periodic oscillations, 
the increment of kinetic energy for one period AK= 0: 


T t.4+T 
F (a; a) dt 8 ddt = 0, (25) 


where T= 22/w is the period of oscillations. 
From the second equation (15) we obtain an expression for a: 


— sin wt + + wt}. 


With the knowledge of a, it is. easy to calculate the first integral 
t+T 


| F (a; a) dt = nw93C (w), (26) 


where the frequency function C (w) = “#4 +08 only depends on the properties of the airplane . 


The second integral of (25) after integration by parts assumes the form: 


T {,+T 
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‘We obtain from equation (20) considering that function f(c) is odd: 


T t,+T 


Using expression (22) to obtain 9, we have: 
J A(K, —K,.) + BoK 
Tore 93 (a9) 2 (27) 


Hence the energy equation (25) can be represented in the form: 


A(K, — + BoK, = —ote(o) (25 
Equation (25°) is an equation of generating amplitudes. 


Let us now find the boundaries of the airplane stability in the plane of coefficients (Kg, K§). Let us multi- 
ply the first equation (15) by sinwt and integrate it with respect to t in the interval to ty to t+ T. We shall obtain 
a linear equation; 


a(K, — + boK, = 0, (28) 


a= AJo—BM, b=BJw+ AM, 


— M® + 


Let us combine equations (25°) and (28): 


A(K, — + BoK, = — °C (w) 1+? 


The determinant of this system does not become identically equal to zero 
a 
Then for every value of parameter w which does not equal zero determinant (29) can be written as 


K,—K,w'= (J — MT,) sin wT, + (JT + M) con 
(30) 


In the first quadrant of plane (Kg, Kg) curve I, determined by formulas (30) where K§ is fixed, is of 
interest, 


where 
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When (09) = const relation (30) becomes, as it can easily be shown, a limiting expression of Gurwitz's 
criterion for linear systems, Thus in the simplest case of an ideal mechanical pilot; 


= K,9 + K,$+ Kz 


and for an airplane with a small lifting forces coefficient Fy ~ 0, and a longitudial movement expresses by 
equation 


relationships (30 degenerate into a limiting Gurwitz condition for third power equations 


(My + + My Ky) (M+ 
at K,>0, K,>9, K;>0. 


in Example. Let the servomotor have a velocity characteristic 


[@=e+e, 


where the plus sign corresponds to a hard and minus to a soft characteristic (Figure 5). 


Considering that o=G,sint it is easily shown that the APFC expression for a hard and soft characteristic 
has the form: 


© =1 +702. 
If &o,4) is known, it is possible to plot from formula (30) a curve [ (w > 0) in the first quadrant of plane 


(Kg, K§). For this calculation Kg remains fixed, 
For the case of K§ = 0 these curves will have approximately the shape shown in Figure 6. 


Fig. 5 Fig. 6. Regionof stability in the plane 


Kg). 


At Kg #0 the region of stability changes; it will be seen from (30) that each point of curve I is displaced 
to the right or left by the amount of Kgw* depending on the sign of K§. 


The shaded region of G in the plane of (Kg, K§) comprises the airplane's region of stability in the long~- 
itudinal direction, 
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It will be seen from (30) that the stability region with a soft servomotor characteristic is larger than with a 
hard one. It can be obtained from curve I, , gg which corresponds to the hard characteristic, by means of a 
geometrical extension by a coefficient related to oo. 


Let us prove that the region G lying in the first quadrant of plane (Kg, Kg) and bounded by curve I (30) 
and straight lines Kg = 0, and Kg = 0, is a region of longitudial asymptotic stability of the aeroplane, For the 
proof let us examine the stability of the periodic solutions which correspond to the position of point (Kgs Kg) 
on the boundary curve Ir. 


Let us assume that a point (Kg, Kg) on the boundary curve is represented by a periodic solution 
Let us apply to the airplane such a disturbing moment about axis Oz as to make the applied power a con~ 
tinuous function of time W= W(t) kG-m-sec™*, Then the power equation will become 
K+F(a; «)— + W (t) = 0. (82) 


It can be easily shown that the disturbing function W(t) can be chosen in such a way as to make the dis- 
turbed movement of the airplane represented by: 


$~4,,sinw(t —t,), (31") 


where amplitude 9,9 # 
The work done by the disturbing moment in the interval (t,, ty+ T) is: 


A, = W (t)dt = {c () + (2,0) [1 + 


The following heuristic criterion of the stability of the periodic solution (31) can be formed. 


For the periodic solution (31) to possess amplitude stability it is necessary and sufficient to fulfill the fol- 
lowing inequality: 


Ay 9.0) > 0, (33) 


i.e., the difference , 9» in the work Ap performed by the disturbing moment must be positive. 


In the case when Ap(%_— Se)) = 0, i.e., when Ap= 0, the solution of (31) is statically unstable with respect 
to amplitude 9, (there is neutral dynamic equilibrium). As an example of the latter case one can cite the linear 
system (0g) = const. 


The transient process in the stable region can be approximately represented if it is assumed that the dif- 
ference 9, ~ $e, is sufficiently small in its absolute value and that 


= sin w(t — (32") 


where h disappears with the difference )~ 9,9. 
Now it is easy to calculate that 


AK = 2nJ,w9%ph + 0(h), (34) 


where O(h”) is a quantity of the same order of magnitude as h’, 


Let us assume that at instant ty the disturbing moment disappears, i.e., that W(t)= 0 atte ty. Then the 
energy equation can be represented approximately in the form 


(K, — + Bok. 


9,) It follows from relation (35) that region C (Figure 7) is the 


region of the airplane stability since for each of its points 
(Kg» Kg) quantity h is negative. 


At the boundary I(30) of region G the solution of (31) is un- 
stable with the respect to amplitude since h > 0 if Seq F Sp. 


If finally (Kg, Kg) is a point in the first quadrant in a region 
outside G, it has (at least) two periodic solutions of the type of 
Fig. 7. Amplitude-phase frequency (31) [3], one of them stable and the other unstable, This is ex- 


isti plained by the APFC of (0,4) having in the majority of cases a 
_ pansies ; shape resembling a normal resonance curve (Figure 7). 


In order toevaluate (09)< = max 4(0,) two values of o are taken Og; and Gg, and the corresponding 
Sq, and So, are determined for which: 


h<0 at rat Gq > 
h=0 at o,= 09, orat = dog, 


h>0O at < So < 


It can be concluded from the above that the mechanical pilot transfer ratios adjusting limits corresponding 
to a stable condition of flying have been obtained. 


A Few General Remarks for Constructing Stability Boundaries, 


It is easy to see that from equations (13) and (30) arise the following very important relations for the further 
study of stability: 


a= Ky = Ke (36) 


On the basis of these relations it becomes possible to reduce the search for the boundaries of stability in 


the plane of two parameters to the solution of one differential equation. In fact from the second equation (7) 
and (10) we obtain: 


2 
— Awk; + 


where Ay = Ay (Ty), Ag = Ag (74). 
A similar equation can be written for longitudinal stability. 


If it becomes possible to prove that relations (36) are sufficiently general, then in a general case the in- 


vestigation of the stability region will be reduced to the solution of a first order linear differential equation of 
the form: 


+ P(T)K; = Q(T), 


d 
where K3 =, The general integral of this equation will be found from the formula: 


K, =e SP gel + 


dK. dK. 
t 
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In this case 
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STUDY OF PERIODIC MOTIONS ARISING IN A SERVOMECHANISM 
OSCILLATORY CIRCUIT AT CONSTANT EXCITATION 


I. N. Krutova 


(Moscow) 


An oscillatory circuit with two parallel control channels is examined. The proc- 
esses in the circuit are represented by three first order equations. Dynamic investigations 
are carried out by means of a phase plane and the point transformation method, The sta- 
bility and purity of oscillations arising in a circuit with constant excitation, the parameters 
of oscillation and the static circuit characteristics are determined. 


The present work is a continuation of article [1]. It deals with a mechanical pilot servomechanism oscil- 
latory circuit [2], whose diagram and short description of operation principles are given in [1]. 


The oscillatory circuit consists of two amplifying tube stages and relays connected in the anode circuits of 
the last amplifying stage, which is separated into two channels, All these elements are affected by a feedback 


circuit which has a link providing a delay in the variations of the signal, when the feedback circuit is closed or 
opened. 


In the presence of excitations, which exceed the dead zone values, at the input of the ocillatory circuit, 
the circuit begins to oscillate. The oscillation parameters depend on the value of the input signal, Variations 
in the oscillation parameters lead to changes in the energizing and release times of the output relays, Thus, a 
relationship is established between the energizing and release times of output relays and the value of the input 
signal, this relationship represents the static characteristics of the oscillatory circuit, 


When the dynamic characteristics of the servomechanism were measured, it was found possible, due to the 
great difference between the speed of the servomechanism and the speed of the oscillatory circuit parameter 


variations, to neglect the transients arising in the circuit when changing fromone stationary condition to another 
and to consider its static characteristics only. 


When an oscillatory circuit is designed, it is necessary to determine the stability and purity of oscillations 


existing in the circuit under oscillatory conditions and to find the parameters of oscillations and the static charac~- 
teristics of the circuit. 


All these questions are solved in the process of a general analysis of oscillatory circuit dynamics at con- 
stant excitation, Since the oscillatory circuit has considerable nonlinearities, this problem is solved by means 
of a phase plane and the point transformations theory developed by A, A, Andronov (3, 4}. 


Posing of the Problem, Equations of Motiog 


In posing the mathematical problem, whose solution provides a dynamics analysis of an oscillatory circuit 
with constant excitation, a dynamic model of the circuit, whose schematic is shown in Figure 1, is examined, 
The dynamic model is made on the basis of the following assumptions. 


1) The first amplifying stage where the external signal is added to the feedback is nonreactive and linear. 
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Fig. 1 
2) The second stage is also nonreactive and possesses nonlinear static characteristics F, (¢) and F, (¢ )shown 
in Figure 2, 
3) The outputs of the second stage [1] have at ¢ > 0 


a) 110) and ¢ < 0 different inductances and resistances since under 
these conditions the outputs have different closed circuits 
in which current is flowing. 
4) The nonlinear functions 0;) and represent- 
ing the operation of the circuits,are shown graphically in 
Figure 3. 
Fig. 2. 


5) The feedback circuit has capacity and resistance, 
The following system of equations represents this 


dynamic model, 
ds 
T, ++, =€, 
bed 
+2 =9, 
+o, =0, 
at 
7, + (2) 
+= O(a, 2), 3) 
(4) 
i at @O(c,)=1, O(c.) —0, 
| 0 at D(s,)=0, =0, 


The system of equations (1) to (5) is written in a relative form with the following notations: 


i is ‘00 


#(0 4, O,) is a nonlinear function valid in the absence of a simultaneous energizing of the relays [1]. 


Weer 


Here [1] iy and i, are currents in the output circuits 
06) of the second amplifying stage, Ugg is the feedback voltage, 
i Auy is the voltage on the grid of the first stage tube, upy is 
| e, input signal yoltage, uo is the supply voltage, L and R are 
+7 z G % the inductanceand resistance when the tube is operating and 
and R’ when it is blocked, Cs and Ry are the capacity and 
—_ resistance of the feedback circuit and b is the feedback co- 
efficient. 


Pig. 3. Phase Planes, 


It will be seen from (1) to (3) that the presence of a 


constant signal at the input of the oscillatory circuit does not 
change its movement equations as compared with the condition of free oscillations, Hence the phase space Oy, 


Og. n preserves the same form as in free oscillations [1] and the problem can be reduced to a simultaneous study 
of two phase planes 0), and Oy, which are coordinate planes of this space. 


Excluding variable t from equations (1) to (3) and integrating them, we obtain equation of surfaces whose 
intersection forms phase trajectories: 


atg> 0 


To — 
< 0 
a= {— — (0)] + — 10 (0) +8} 
— [n— + 
(8) 
1 = S19 (9) 


where a= Ts/Ty b= Ts/T2, Og Ogg no are the initial values of (0) = 


Phase planes 04,7 and G,,n are multi-layer since each of the phase space volumes 04, G, and n corresponding 
to a definite value of &o) is projected onto a separate sheet. Phase surfaces comprise six sheets, Straight lines 
€y = E/y serve as transfer lines on plane and straight lines 0, = 0,= €, 1 on plane 
Oy. Phase planes 0y,n and 0, with sheet numbering are shown in Figure 4. With the existing distribution of 
volumes in the space 0, Oy» sheets I, II, IV and VI on plane 0,, » and sheets I’, III’, IV’ and VI" on planeo,,y 
are completely superimposible and have similar boundaries. The phase trajectory projections on plane 6,7 and 
0,7 are determined by equations (6) to (9). 


Let us see now what movements of a representative point on phase planes can arise with constant signals of 
different values. For normal operation of a servomechanism it is essential to have a stable oscillatory circuit, 
therefore, the study of its dynamics will be carried out with the assumption that at £ = 0 its motions are stable at 
all the initial deviations, This means that on planes ©, and 6,, n the phase trajectories twist and the move- 
ment of the representative points tends to a state of rest with increasing time. We shall study phase planes with 


positive signals only, since with a change of sign the picture of the movements on planes o,,7 has the same form 
as that on plane 0,,n at & > 0 and vice versa. 


In the state of free oscillations [1] two regions were separated in the phase plane; the region of attraction 
to the rest point and region I, where’ the representative points traverse in their movement all the six sheets of the 


543 


| 
| | 
| 


Fig. 4 


phase plane and then return to their initial sheet. The presence of a constant signal leads to a transformation of 
movements inside those regions to an alteration of their boundaries and to the appearance of new forms of move- 
ment. 


Let us examine region I, At g¢ = 0 the phase trajectories twist in this region. With an increasing signal the 
movement is transformed; it is necessary to establish the character of movement changes in this region. Figure 4 
shows the movement in region I with initial deviations o, and 0». 


At n > &/y phase trajectories vary in the plane o,,n. It is possible to determine from the phase trajectory 
equations for sheets I’, I", and III" that in a fixed range of variations of n anda, derivatives 0093/3E, 3 Og, /3E 
and 0f/d¢ are negative. The notations here introduced are: Gs and Gy, are the values of a, on sheets III" and 
Ii” and f is the ordinate of the crossing point of the trajectory of sheet f and the straight line g=€,. Thus with 
an increase in the signal strength the trajectories in these sheets tend to straighten out. This gives weight to the 
assumption that with an increased signal strength the general movement trajectories in the upper semi-plane 
GO», n twist even more. In order to confirm this assumption let us examine the movement of a point in the upper 
semi~plane consecutively over the three sheets with initial and end points lying on the straight line n =&/y. 


Having divided this movement into two sections one of which includes the trajectories of sheets III’ and I’, 
and the other the trajectories of sheet II’, let us find the expression for determining f and 04: 


Sill 


(0) 


on = teat (11) 


Whence derivatives 0f/0¢ and 60,,/d¢ can be found in the form 


4 y/ Taalteatsé 


+e: + &) 


‘ 


The analysis of expressions (12) and (13) shows that both derivatives are negative and the final value of o, 
decreases with a rising signal strenght, since the decrease in f leads to a smaller o,,, Hence in the upper semi- 
plane o,, n the trajectories twist with a rising signal more intensively, 


In the upper semi~plane 0,,n,trajectories do not change, but the variation of the transfer line n = &/y and 
the lowering of the transfer point F ordinate leads to a decrease in the final value of o;, 


At n <€/y trajectories do change in the plane oy. It is easy to determine from the phase trajectory equa- 
tions of sheets I, I, and II that derivatives d0,,/d¢ and d0,3/d¢ , where oy, and o;, are values of o; on sheets II 
and III, are positive, but derivative dg/€, where » is the ordinate of the crossing point of sheet I trajectory with 
line 0, = €», is negative. Hence in the lower semi~-plane oy,n the trajectories acquire a greater twist with a rising 
signal. By examining the movement in the lower plane consecutively over the three sheets, similarly to the man- 
ner adopted for the upper plane, we find the expression for determining g and 04): 


Whence 0g/d_ and 004, /¢ can be found in the form: 


etat+ 


An analysis of expression (17) shows that 804, /og is positive. Derivative 0g/d¢ can be positive or negative. 
By comparing (16) and (12) it can be seen that 4g/d¢ is less negative than Of/ dg even at Oy = Og. In fact 
Oy>Oy9 and this inequality increases with a rising signal, since the final value of o, decreases in the upper semi- 
plane (or that of ogo in the lower semi-plane), Hence the decrease in f is greater than in yg. Comparison of 
derivativesdo,,/d¢— and 90,,//¢ shows that they are close to each other in their absolute values. 


The final yalue of o, in the lower semi~plane o,,n can either decrease or increase depending on the sign 
and value of a, in the lower semi-~plane 0,7 is smaller than the fall of o, in the upper semi~plane 0,7. 


Thus, the analysis of phase trajectory changes in region I establishes that a fall in the signal leads to a 
More intense twisting of the trajectories in plare ¢,,n, but in plane 0,, n two cases are possible: the twisting of 


(16) 
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trajectories can be either more or less intense than in the case of free oscillations. In the latter case, however, 

the rise in o, is always slower than the fall in o,, Hence the presence of a constant signal leads to an increase in 
stability of region I. This conclusion is confirmed by a number of numerical calculations made for different 
relations between parameters (in the region of a < 1, b >a) and different signal values. The analysis of trajectories 
also leads to the conclusion that with a rising signal the boundaries of region I shrink. 


Let us now study the region of the point-of-rest attraction, With ¢€, > > 0 plane o4,n has only a special 
point of sheet I and plane o,,n a special point of sheet IV’, Hence with any initial deviations in this region the 
representative points inevitably fall on these sheets, moreover the movement can occur in a general case with 
switching of the relay in one of the channels or in both channels, Figure 4 shows both types of motion with ini- 
tial deviations o,***, o,*** and 0, **, o,**. Boundaries of R and S can be shown for these types of motion 
on straight line n = €/y in a manner similar to the one used for free oscillations [1]. It follows from the analysis 
of phase trajectory variations that with a falling signal the boundaries of R expand and those of S shrink. 


When € > ¢, the special point of sheet I falls outside the limits of plane o4,n. At the same time all the 
trajectories of sheet I cross the boundary 0; = €,, and the representative point in the process of movement always’ 
passes beyond sheet II and then returns to the sheet from which it started moving. Thus, the region of the rest~- 
ing point attraction is transformed into a region of movements of a new type during which the representative 
point traverses either four sheets, if it reaches the boundary n = &/y, or two sheets, if the boundary n = E/y 
rises sufficiently high and the trajectory of sheet II crosses boundary 0, = €, lower than line n = &/y. It is pos- 
sible to separate in the phase plane, region I, in which the representative pointtraverses four sheets, and 
region Il, in which it traverses two sheets. Figures 5 and 6 show movements in regions II and I, Periodic 
movements always exist in these regions. 


| 


At — > y + €; the boundary n = €/y rises higher than n «1, The special point of sheet Il appears on plane 
Oy n, this point is stable and all the movement tends towards it, since the representative point on plane G,n 
for any initial deviation falls on sheet II, In this region periodic motions do not arise in the circuit. Signal 


& = y + & corresponds to the end of the oscillatory condition since the input signal becomes higher than that of 
the feedback. 


From the analysis of movements existing on the phase plane at constant signals of different values it ap- 
pears that the solution of the problem of limiting cycles can be reduced to the study of the division of regions 
Iland Ill into trajectories, The peculiarity of region II and Il consisting in the fact that the behavior of the 
contour in these regions is controlled by processes in one channel only, since in the other channel the relay is 
not switched and therefore only one phase plane can be considered, 


Point Transformation 


The analysis of the phase plane division into trajectories in regions I and II can be reduced to the study 
of point transformation of a straight line into a straight line. In this case it is sufficient to consider the phase 


plane o,,n with positive signals, since with negative signals the picture of movement along plane 9,7 will be 
similar, 


When studying the division of phase plane o;, n into trajectories it is convenient to regard the transformation 


of a straight line in o, = €, into itself, since all the trajectories of sheets I, Il, IV and V cross th:. line. Point 
transformation can be composed cf two correspondence functions which determine the transition of points on 
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straight line o,= ¢€, to points on straight line o,= €, (function A) and points on straight line O,= €; to points of 
straight line o,= €, (function B). Due to the existence of points of contact of straight line o,= €, with the phase 
trajectory of sheet II and of straight line o,= €, with the phase trajectory sheet I, point transformation holds for 
the stretch of the straight line o, = €, from the point of contact (n = u*) downwards, and for straight line 9)= & 
from the point of contact ( n =v *) upwards. On the straight line o,= €, it is possible to mark- off the boundaries 
ny, 4nd na, which divide the three regions of movement, The boundary ny between regions I and II is determined 
from the condition that the transit point A* must lie on the plane 0,,n inside the region bounded by the trajectory 
of the touching line o, = €, (Figure 4). If the ordinates of the crossing points of the trajectory passing through 

A’ (q*) and the trajectory of the touching line o, = €, (n,,) with axis n, be compared, condition n*= np will be 
found to hold for boundary 73, from which it is possible to obtain. by means of phase trajectory equations: 


i—a 


The boundary condition between movements of the second and third type is provided by the case when the 
ordinate of the transition point D (Figure 5) assumes the value of n = €/y, i. e., whenn,=€E/y. 


Let us now examine point transformation in region Il, Let o,= €, and » = u at the instant the movement 
begins, The value of coordinates on sheet II will obviously be 0, = o,,, n = &/y, Om sheet V it will be o,= «, 
n = v, on sheet IV it will be o; = €y, » = &/y and on sheet I it will be o, = €,, 7 = uy. By means of the phase 
trajectory equations and excluding the intermediate values of coordinates, we find the relation between u, v for 
function A and uy, v for function B in the form: 


(18) 


Derivative dv/ du is represented by expression 


When u varies between ~@ and u*, quantity v varies between 1 and &/y. The derivative dv / du is neg~ 
ative atu<u* = and positive at u> u*, 


The expression for the derivative dv/ du, has the form; 


1 
It follows from (22) that the type of variations of function B are similar to those of function A, since both 
(21) and (22) containe the same polynomial ¢ ~ yu~ €,, which determines the sign of the derivative, With 
the variation of v from I to ¢/y the quantity uy changes in the limits which we shall denote m and n, On the 
basis of the analysis of the phase plane it is possible to determine that m> 0 and n< u®, since at u,= u® the 
phase trajectory of sheet I only touches the line o, = ¢, and the above transformation becomes impossible. 


Analyzing the mutual position of functions A and B graphs,we can write the derivatives in the following 
form; 


(22) 
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(e,)'” (E — — (=y"" (24) 


In order to obtain a transformation of the above form the following conditions should hold: oy > «, 
0 < Oy, <q. Moreover in the band of variations of uz; and uy in region II there exist inequalities 


tr < ee > 1. It follows from this that in the band of = [< bed and hence in region II func- 
i 


tions A and B can only have one crossing point or no crossing point at all, 


Let us now study the region of movements of the third type. The same point transformation as in the pre- 
ceding case will be determined by two correspondence functions A and B, Denoting the initial value of n on 
on straight line 0,= ¢ by u and the end value by uy, and the value n on the straight line o= « by v, let us find 


by means of phase trajectory equations the following expressions relating u and v in function A’ and u, with y 
in function B"; 


@— 1) = 0, (25) 
(Gm 

The expression for the derivative dv/ du has the form: 


Derivative dv/du becomes zero when u = u® and @ when v = v* = ¢ ~ &/y. The latter equality is, how- 
ever, impossible for function A’ since at v = v® the trajectory of sheet Il only touches the straight line o, = & 
and the transformation in question does not exist. Therefore y > v* always holds and the denominator in ex- 
pression (27) is greater than zero. Function A‘has the same form as function A since with u > u® the derivative 
dv/ du is negative, and with u < u® it is positive. The expression for dv/ du, has the form: 


dv 


wra—e \o (28) 


Fig. 7 


The first factor in (28) is the same as in (27), However, contrary to function A‘, in function B’ the derivative 
dv/ du can assume a value of ow at v = v*, but cannot assume the value of zero, since in this case at u,= u® the 
trajectory of sheet I only touches the line 4 = «, and the above transformation is impossible, 
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In region II quantity v is always larger than u or u, and makes valid the inequality: 


d dv |, 
lal 
Thus, the investigation of the behavior of correspondence functions A, B and A*, B" shows that in the region 
of the second and third type movements they can only have one crossing point which corresponds to the stable 
limiting cycle, Hence withexcitation lying in the limits of ¢«,< < €,+ y the oscillatory circuit will always 
be in state of stable oscillations which can be of the second or third type, depending on the value of the input 
signal, The existence of one crossing point in the curves of the point transformation diagram shows that for each 
value of the signal there can only be one form of oscillation. 


Figure 7 shows a diagram of point transformations, The diagram has been plotted for the case when the 
curves determined by the correspondence functions have a crossing point in region II, 


Static Characteristics of an Oscillatory Circuit. 


The static characteristics of oscillatory circuit tr, = fy (€) and ty = f, (€), where t,, and tp are the making 
and breaking times of the output relays, are calculated on the basis of the sustained oscillations parameters u° 
and vo where u° and v°are amplitude values of n (Figure 7). By means of the correspondence functions we ob- 
tain for u° in terms of sustained oscillations of the second type the following equation: 


| 


v° is determined from equation (20), u° and v°for sustained oscillations of the third type are determined ir 
from the following equations: 


The values of t,, and t, are determined by solving equation (3): 


i—u° 


tm= Tmin (32) 
th ‘ (33) 


The period of sustained oscillations T = t,, + tp. 
The static characteristics of the oscillatory circuit,calculated for certain parameters,are shown in Figure 8. 


The static characteristics obtained theoretically are, however, incomplete, since with the assumption made 
for representing the processesoccurring in the circuit, it is impossible to determine the zone of ambiguity whose 
presence may be partly explained by residual magnetism and inertia of the relay moving parts. Due to this 
phenomenon the operating current of the output relay with a rising control signal (¢ ,) is larger than with a de- 
creasing signal (¢;). In investigating the dynamics of servomechanisms it is necessary to consider the full static 
characteristics since the zone of ambiguity can produce a substantial effect. In Figure 8 the static characteristics 
Of tm f; and t, = f, in the region of < are shown by a dotted line, 


2) 
3) 


: 
Fig. 8 
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DETERMINATION OF TRANSFER FUNCTION COEFFICIENTS OF 
A LINEAR SYSTEM FROM THE INITIAL PORTION OF AN 
EXPERIMENATLLY OBTAINED AMPLITUDE-PHASE CHARACTERISTIC 


E. E. Dudnikoyv 


(Moscow) 


The article deals with a method of determining approximate values of transfer 
function coefficients of a linear system from an experimentally obtained amplitude~ 
phase characteristic, The coefficients are determined from the initial portion of the 
characteristic (at w-» 0). The method is applicable to delayed systems. A numerical 
example of determining coefficients is given. 


At present there exist several methods of determining approximate values of transfer function coefficients 
of linear systems from experimental data. Sufficiently well developed and effective methods of determining 
these coefficients from experimentally obtained transfer characteristics of systems [1, 2] are known. In analyzing 
linear systems it is, however, often necessary to deal with experimental amplutude - phase characteristics. Sev- 
eral papers [3 to 6) deal with methods of determining transfer function coefficients of linear systems from experi- 
mentally obtained amplitude-phase characteristics. In all these methods, calculations are preceded by selecting 
a structure for the analytical expression of the system's transfer function. It is assumed, for instance, that the 
amplitude -phase characteristic of the system has the form 


K W (jo) = + Bio) K 


K 
= Ga)" (1 + Ajo) * 


(1) 


where K, B and A are constants, and D(jw) is a polynomial in jw. 


In [5] the transfer function of the system is approximated by a transfer function of an equivalent system 
of the second or third order with delay. In [6] transfer functions of even higher order systems are used for ap- 
proximation with subsequent calculations of corrections for the coefficients; this method does not require labor- 
consuming calculations and provides good results even with the first approximation. In this method, it is also 
necessary, however, to assume in advance the form of the transfer function. This circumstance in many problems 
considerably lowers the effectiveness of application of all the known methods. 


In this paper an attempt is made to calculate the approximate values of transfer function coefficients of a 
linear system without assuming in advance the structure of the required function, Essentially this method is 
similar to M, P. Simoyu*s method [2] for transfer characteristics. 


In a general case a system with lumped parameters is represented by an equation of the form 


+ yp + 1) tour = (2) 
= (bap + + be) Bin 


4 
| 


where n > k, agby (s= 1, ..., nz 2 = 0, ..., k) are constant coefficients, and p is operator tae 
Such a system has a transfer function 


k 


Wo(p) = i+ (3) 
1+ 


Let us multiply the numerator and the denominator of the amplitude~-phase characteristic equation of the 
system under consideration by 1/(jw)" 


W, (jw) = by + br by, + by Gor 


4 
(jo) + a; (ja)! +----+ a, _, + 
Let us convert this expression to the form: 
Wo (jw) = Ay + ®, (jo), 
where Ag is a constant, and 
Cy (ja)! + +--+ + + (fay? +6, 
(jo) = . (8) 
+ (jo)! + @y_, (ja)! +4, 
Let us transform function 
= 
separating the term which does not contain jw and the one which contains (jw)”* 
W, (jw) = Ai + By (jo). (8) 


Function @,(jw) has a numerator of the order(2~n) and a denominator of the order (1—n); 


From function @,(jw) we find its inverse function W,(jw) and transform it similarly to W,(jw) separating 
function $5 (jw), etc. 


As result of these operations the initial amplitude-phase characteristic W,{jw) can be represented by the 
expression 


W, (jw) = Ay + 
B, (ja)? + Ai + (9) 
By, (ja)? + Astro, 


i 
+ Ay 


Constants Ag, By, Ay By ...» By - Ap are determined as follows, Let us represent function Ws ($= 0, ..., ) 
in the form: 


W, (jo) = Po (@) + (@), 
W, (jw) = P,(@) + 78, (@), 


We (jw): = Px + (@), 


where Po(w), Py(w), ..-» Pn(w) are the real parts of these functions, G9 (w), G,(W), ...» On(w) are expreisions for 
the imaginary parts of the functions, 
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Let us assume that w is tending to zero, Then | Wo(jw)| will tend asymptotically to | Py (w)|and the latter 
according to (5) to | Ag] =| bg| ; | Wy(jw)] will tend asymptotically to | @,(1)| and the latter according to (8) 
to | By(jw)™*] ; | Wy (iw) — Bo(iw)~*| will tend asymptotically to | Py(w)| i.e., to | Ay]; | W2(jw) will tend 
asymptotically to | @w)|, i. to |By(jw)™"]; | W,(jw) — By(jw)~"| will tend asymptotically to |P,(w)]| , t.e., 
to | Wn(jw)| will tend asymptotically to | , to | and | — Bp 
will tend asymptotically to | Pp(w), i.e., to | . 


Fig. 1 


The determination of the constant coefficients begins with the finding of Ap. For this purpose the asymptote 
to which | P,(w)| is tending at w > 0 is plotted; the length cut by the asymptote on the Y-axis is equal to | Aol. 


The next stage is to plot function ('7) 


4 4 


The graph of function | 6,(w)| is plotted on a full logarithmic scale and the asymptote to which | 6(w)| 
is tending at w -> 0 (the tangent of the asymptote’s slope angle with respect to the X-axis is equal to ~1), The 
asymptote cuts off a length on the Y~axis equal to log| By| . 


Next the graph of function | P,(w)| is plotted and from it coefficient | A,| determined by the same method 
as was used for obtaining | Ag] . 


Function #,(jw) is calculated and from it function W,(jw) is plotted and coefficients | B,| and | A,| are 
determined, etc. 


The signs of Ao; jBy/w; Ay ..-» JBy-;/ws An are taken as those of the corresponding imaginary parts of the 
function from which they were derived, 


The end of the calculation is indicated by the given function $p44(jw) becoming equal to zero, 


After the last coefficients have been determined the entire function is written out in the form of expression 
(9) and is reduced to a common denominator, As the result of this the required expression of the amplitude -phase 
characteristic of the system is obtained. 


In this case the initial portion of the amplitude~phase characteristic was used for determining the approxim~ 
ate value of its coefficients (the coefficients were determined at w > 0). It is possible to use another version of this 
method when the coefficients are determined at w-> @. 


w=40 
w*30 
w=0 
was 
w=06 
| 
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If the system under investigation has a delay r, the value of this dealy is first determined. The determina- 


tion of the delay can be carried out either from the experimentally obtained system's transfer function or from its 
phase -fre quency characteristic. In the latter case, an asymptote to which the phase-frequency characteristic is 
tending at w-»+o is plotted. The tangent of the slope angle of the asymptote to the X-axis is equal to the de- 
lay. After + has been determined an amplitude~phase characteristic is plotted without a delay, This character- 
istic is then operated upon by the method described above, The thus calculated transfer function is then multi- 
plied by the transfer function of the delay link, 


Example. It is easy to demonstrate the basic advantage of this method on a numerical example. Let 
us assume that the transfer function of a linear system is: 
0.0192 p + 0,385 
w 0 
= 5 (1 


where r = 0,5 sec, but that it is given in the form of a table of real P(w) and imaginary @(w) values of its frequency 


response calculated with an accuracy of three significant figures (Table 1), 
The amplitude-phase characteristic of the system: 


W (10) = =P + 18 (0) (11) 


is shown in Figure 1. 


It is required to determine the structure and value of the coefficients of the system's transfer function by 
using the data on Table 1 only. 


TABLE 1 


cs) 0 0.2 0.4 0.6 0.8 1.0 1.5 2.0 2.5 3.0 4.0 


P(w)|0.385| 0. 0.399 0.275 0.2 
@(@)| 0 |-0. 247 305-0. 2-0-0. 0.014] 0,078 


o 5.0 | 8.0 10 15 20 2 30 Ss) 40 50 


P (@) 0 ono 0.011|—82-10-*/75- 10-*/48 -10-*/36- 4 “1 28°10-4| —92-10-* 


6 (@) 0.018 58-10-4 peters 6-10-* | 16-40 


The value of delay is determined from the phase~frequency characteristic graph (Table 2 and Figure 2); 
in this case (Figure 2)r is equal to 0,5 sec, 


An amplitude~phase characteristic of the system without delay is then plotted (Figure 3). 
Wo (j@) = Po(w) + (@). (12) 


From the graph of function | Pg(w)| (Figure 4)* coefficient Ag is then determined; in this case Ag= 0.385, 


*For convenience of lay-out functions,| P,(w) ||, | Py(«9] and | P,(w)| are plotted on a logarithmic scale along 
both axes. 
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quurad Next it is necessary to calculate P,;(w) and 6;(w), the real 
part and the coefficient of the imaginary part of function. 


sk W, (jo) = 
arctyr 
From the graph of function | P;(w)] (Figure 4) coefficient 
Ay of function W, (jw) is then determined; in this case Ay= —1,81, 


t w te The graph of function | @,(w)| is then plotted on a loga- 
: rithmic scale (Figure 5and coefficient B, determined; in this 
Fig. case By= 3.61. 


Fig. 3 


Then it is necessary to calculate P,(w) and ©,(w), the real part and the coefficient of the imaginary part 
of function 


W;, (je) = - 


From the graphs of function | P,(w)| and | 6,(w)| (Figures 4 and 5), coefficients A, and B, of function 
W;(j) are determined; in this case A= ~1,27, and By= 4.75. 


This completes the calculations since the last function $3(jw) is equal to zero. 


TABLE 2 


w |o}os | o6 | o8 | 1.0] 15] 20] 3.0] 5.0] 80] 20 |» 


| 0. 0, 0.7 0. 1.24 | 1.79 | 2.32 | 3.23 | 4.60) 6,27 


All the calculation results are entered in Table 3.¢ 
Thus 


1 
We (jo) = 0.385 + 


4 (13) 


*Since at frequencies w =< 1.0 the differences | 6,(w)—By/w| and | P,Q@)—A,| are very small the plotting of 
function W,(jw) must be started from w =1.0in order to obtain accurate results. 
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After transformation and taking delay into account we obtain 


W (je) = — (14) 


Fig. 4 
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Fig. 5 
Thus the calculation results were obtained with the same accuracy as the initial data. 
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APPLICATION OF SIMULATION IN ANALYZING LINEAR PULSE 
SYSTEMS WITH VARIABLE PARAMETERS 


G. P. Tartakovskii 


(Moscow) 


Simulation of linear systems with discrete operation and variable parameters with 
the object of obtaining pulse response suitable for further application is examined. 


It is shown that a pulse response as a function of the instants when pulses are ap~- 
plied, can be obtained by simulating a pulse system conjugate to the original system, 
Methods of designing schematic circuits for conjugate pulse system models,similar to 
corresponding methods of simulating continuously operating systems,are derived, 


Posing of the Problem 


In analyzing discrete action systems, represented by linear difference equations with variable coefficients 
considerable calculating difficulties often arise. This especially applies to complicated interrupted control 
systems with feedback represented by higher order equations. Therefore, in a number of cases, analytical in- 
vestigation methods have to be replaced by numerical ones. In particular, simulation methods can be applied to 
discrete action systems in a manner similar to one used for continuously operating systems, In this connection 
it is possible to simulate either the pulse system, consisting of the pulse element and the linear continuously 
operating portion, or to simulate the difference equation directly. In the first instance it is possible to use a 
continuously operating simulator, connected to a pulse element. In the second instance the model can be made 
up of registers, which accumulate discretely changing quantities, adders and amplifiers with discretely changing 
gain. The choice of the simulation method depends on the type of the problem to be solved, 


It is convenient, for instance, to find by means of simulation the pulse response of system h(n, k) [1).* 
From this response it is possible to find either by means of further simulation or analytically ** all the prop~- 


erties of the system in question. 
In this connection the response to an arbitrary regular signal ujn({n) is found from; 


Wout (n) = k) (). 


The correlation function of a discrete random process at the output of a pulse system with regularly chang~- 
ing parameters can be written in the form: 


n n+m 
Kout (n,m) = h(n,t) >) h(n+m, 9) Kin 


*Function h (n, k) represents the response of a system at discrete instants tp to a pulse of unity amplitude and 
fixed duration r (t,), occurring imstant t,; n and k are dimensionless integers. 


** In the latter instance the pulse response should be approximated by a discrete function. 
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K (n, m) = M — M, (n)| [x(n + m) — M, (n+ m))}}, 


where M denoted the mathematical expectation, My the mathematical expectation of quantity x. 
In the particular case of a discrete uncorrelated random process at the output of a pulse system 


Kin (n, m) = Bm, 


1 at 
n= at m0. 


The dispersion at the output is 


(2) = Kout (n, 0) = D> 


Thus the analysis of both regular and random processes in pulse systems can be carried out, if the pulse 
response is known, 


In formulas (1), (2), and (6) the pulse response h(n, k) appears as a function of the instants when the pulses 
are applied to the system, When: this system is simulated, h(n, k) is obtained as a function of registration instants 
n with a k fixed for each experiment. Thus, the same difficulty is encountered as in the case of simulating con- 
tinuous systems with variable parameters. 


In a manner similar to the one adopted for systems represented by differential equations [2,3,4] this dif- 
ficulty can be-overcome by simulating a pulse system conjugate to the original system. The response of such 
conjugate system to a unit pulse coincides with h(n, k) as a function of variable k. 


The aim of this paper is to find a method of compiling equations and schematic circuits of conjugate pulse 
systems for both the above mentioned methods of simulation. 


Finding a Conjugate System for a Model Consisting of a Pulse Element and a Linear 
Continuously Operating Portion. 


In order to find a difference equation for a pulse response h(n, k) of a discrete action system with variable 
parameters as a function of k, it is possible to use a method similar to the one used in paper [3] for continuously 
operating systems, This technique is based on the phenomenon that the motion of a system, whose input disturb- 
ance starting from a certain instant is equal to zero, is completely expressed by the initial conditions and does 
not depefid on the method by which they are obtained. 


For aiy pulse system with signals in the form of constant quantities applied to appropriate points of its 
circuit, it is possible to choose such initial conditions that the output signal will be in the form of a pulse response. 
Moreover, these initial conditions can be established by applying a certain sequence of pulses to the system's in~ 
put. The response to such a sequence can also be expressed by means of the superposition formula to a pulse 
response, As the result of this, it is possible to obtain the required equation which corresponds to the conjugate 
system, 

Let us describe an example of the application of this method to an interrupted control system of the second 
order (Figure 1)*, which is typical of many technical problems. 


If pulses of amplitude A,(k) and A,(k+ 1) be fed to the input of the system at instants k and k + 1, it is 
possible to select these amplitudes in such a manner that at instant k + 2,at point 1,there will be zero and at 
point 2, unity initial conditions, 


* In Figures 1, 2, and 5 the transfer functions of the system's elements are shown inside the squares. 


: 
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pl, 


Fig. 1, Simulated continuous control system. { 


The response at the cutput of the system at be het hou by q 
formula 
(M) A, (k) h(n, k) + Ay + 1) h(n, + 1). 


(7) 


It is also clear that the same response will be obtained at the output of the system if a single impulse 
1(n—k~—2) is applied to point 2 at instant k+ 2. In order to have an equal response at instant k+ 1 as well, it is 
necessary to add to the output of the system pulse C(k+ 1)8n,k+ 1» Where C(k+ 1) is the value of the output signal 
at instant k+ 1 when a pulse of amplitude A,(k) is applied to the input at instant k, and 6p, , is Kronecker’s 


delta symbol. 
3 n=k, 
10 nk. 


The application of the impulse to point 2 is equivalent to feeding it to the output of the systems and ap- 
plying a negative impulse to the input. Hence 


Uout (n) = 1(n— — 2) + C + 1) 8n, — a(n, & +42), 


(8) 


where a(n, k) is the transfer function of the pulse system, 
Equating (7) to (8) we obtain: 


A, (k) h(n, k) + Ag(k+1)k(n, k+1)+ a(n, k4+2)= 
= 1 (n—k—2)+C (kh +1) 8n, 
Substituting k in (9) by k+ 1 we have 


A,(k + 1) h(n, k+ 1) + Ag(k+ 2)h(n, k + 2)+ a(n, k+3)= 
=1(n—k—3) nae. am 


Subtracting (10) from (9) and considering that 
Aya (n, k) = —h(n, k), 


where A, is the difference with respect to variable k, we obtaim 


[1 — + 2)J h(n, + 2) + (A, (A+ 1) + + & +1) + 
+ A, (hk) h(n, k) + 1) bn, — 2) Bn, 
By introducing a new variable 
E=n—k—2 
and completing simple operations we finally obtain: 


(n — h(n, — §) Ag [Ai (m — &) + h(n, + (14) 
h(n, n— 8) = b 


$e, (11) 
(12) | 
(13) 
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Equation (14) determines pulse response h(n,n ~ ¢) as a function of €. This function is an output signal of 
& pulse system whose input is fed at instant — = 0 by a unit pulse, Such a system according to the terminology 
introduced above is conjugate to the original pulse system. 


Equation similar to (14) can be easily formulated for systems of a higher order as well, 
Reverting to the example in hand, let us find the values of Ay(k), Ag(k) and C(k) 


By solying the differential equation of the links with transfer functions 1/ pT, and 1+ pTk/ pT,, knowing the 
duration of pulse r and the repetition period T;, determined by the pulse element, it is easy to find the value of 
response at points 1 and 2. Moreover at point 1, at the end of the pulse operating at instant k+1, there will be 
signal: 


According to the set conditions, this quantity must equal zero: 


+ 1) — (Pe — + + AV = 0. (15) 


In point 2 at instant k+ 1: 


C (k-+ 1) = (18) 


and at the end of the pulse the output signal will be: 
Ay (Te + Tr — + 
+ [Aah + 1) — (Te + — + (Te + ay 


From equations (15), (16) and (17) we find: 


+1 

2(k + 1) = (18) 


A, (k) = 


r 


By substituting (18) im (14) we obtain the equation of the conjugate system: 


(19) 
k t 
+ h(n, + Me (at +t —sr) 


In compiling a difference equation directly it is easy to see that equation (19) corresponds to an interrupted 
control system, whose schematic is shown in Figure 2, This system is obtained from the original one by inter- 
changing the inputs and outputs of all the linear portion elements, reversing the direction of all the connections, 
interchanging the input and output of the whole system, transferring the pulse element to the input of the newly 
formed system and reversing with respect to time (replacing n by n—¢) the operation of the element, which varies 
with time. 


ement 
pl, Pl, 


Fig. 2, Interrupted control system, conjugate to the system of Figure 1, 


Conjugate Systems for Simulating Linear Difference Equations with Variable Co- 
efficients. 


Let us assume there is a linear pulse system with variable parameters, whose difference equation has the 


form: 


Dan) A ‘tout (n) = by (n) A‘ (n), v<m. (20) 


A model of such a system can be made up of registers, adders and amplifiers with discretely changing gain. | 
A step function should be applied to the input of the model, Let us find the equation for, and the means of ob- | 
taining, a conjugate system from the same elements. j 


Let us first examine the case when bj(n) = 0 at i 0 and be(n) = 1. 
The system*s equation then has the form: 


> a; (n) = Uin (n). (21) 
If the pulse response of this system is denoted by hy(n, k) we have: 
(2) = y(n, (22) 
l==2—oo 


Having in mind the results obtained for the particular case of the preceding paragraph, itis natural to as~ 
sume that an equation for the pulse response will be obtained if it is assumed that ugy¢(M) ™ 6p. ye 


Then 
i=o 
It is easy to find that | 
= (— — i). 
Hence 
iat 
[Ay & —i) ay (k (25) 
By introducing a new variable n- k™£€ and considering that 
(— (A = (n —8), (26) 
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we obtain an equation for the pulse response h,(n, k) as a function of the instants when the pulses are applied: 
Attia (n, (n = (27) 
1=0 


This, in fact, is the equation of a system conjugate to the system represented by equation (21). 


6,(n) | 


Fig. 3. Schematic of a model corresponding to the difference equation (20); a is 
the element which registers differences of the order of i, Hj is the i-th storage 
device, b(n) and aj(m) are amplifiers with discretely changing gain. 


The system described by equation (20) can be represented as a series connection of systems with equations; 
= >) bi (nm) (28) 
i=o 
and systems with equation (21). According to the superposition formula; 


Bout (2) = >> hy (n, (4), 


l=—oo 


h(n, k)= Ay(n, (L) x. 
l=—oo 


Substituting (24) and (26) in (30) we find 


h(n, = di (n—8)]. (31) 


i=0 


Above equations provide a simple method of designing a conjugate pulse system model from the model of 
the original system. The original system could, for instance, be as the one shown in the schematic Figure 3, It 
follows from equations (2'7) and (31) that the conjugate system is represented on the schematic of Figure 4, 


Thus, for designing a conjugate system model from a model of a linear difference equation, it is sufficient 
to perform the following operations: to interchange the inputs and outputs of all the elements, to reverse the 
direction of all the connections, to interchange the input and output of the system as a whole, and to reverse the 
gain response of the amplifiers with respect to time.* A unit pulse is fed to the input of the new model at instant 


With such a building-up method at points where in the original circuit there is parallel feeding of signals to 
several elements, in the conjugate system circuit there should be an adder and vice versa. 


| 
| 
or the pulse response 


t = 0, where ¢ serves as the current dimensionless time. Then the output signal becomes identical with pulse 
response h(n, 


Fig. 4. Schematic of a model conjugate to the model in Fig, 3, 
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Fig, 5, Example of a simulated system with an interrupted control, 


Fig. 6. Model corresponding to the system in Fig. 5. H, and H, are 
storage devices,f(n) is an amplifier with a variable gain, Ky Ky Ky 
and K, are amplifiers and 1 is an inverter, 


For each given value of the counting time n the experiment is repeated, Thus, a complete analogy exists 
between this and the well-known method of building up conjugate systems for continuously operating systems, 
with the exception that integrators are replaced by registers of discretely changing signals, the amplifiers have 
a discretely instead of smoothly changing gain,and a pulse of unit amplitude (in dimensionless units),instead ofa 
delta function,is fed to the input of the system, 


The inclusion in the circuit of elements, which register differences between discrete functions is, of course, 
optional, By using series and parallel connections and feedbacks, it is always possible to construct a model cor- 
responding to equation (20) without using such elements, It is also possible to show that the rules compiled for 
building up a conjugate system model, will remain in force in this case also, 


As an example, let us examine the system, whose schematic is shown in Figure 5, Such a system is common 
to many practical problems, By compiling and solving the differential equations of the links of the system in 
question, it is easy to obtain recurring relations expressed by equations: 


Ax(n) = f(n)uout (n) + mn (n), 


Awour = K Ba (n) — [ — Tye wour + (82) 
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r is the pulse duration and T, is the repetition period. 


Fig. 7. Model of a system conjugate to the system shown in Fig. 6. 


If we denote 


then equation (32) will correspond to the schematic circuit of the model shown in Figure 6, By means of the 
rules outlined above, we can obtain from the circuit of Figure 6 a circuit of a conjugate system model (Figure 7), 
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INVARIANCE PRINCIPLE IN COMBINED AUTOMATIC 
CONTROL SYSTEMS 


I. Dunaey 


(Moscow) 


A further development is given of the complete invariance principle [1-9] ap- 
plied to combined automatic control systems with two motors driving a common load, 


Automatic control systems (ACS) should be designed to provide in practice high quality devices independent 


of load variations and random input signal disturbances. The systems should provide high tracking speeds and ac- 
celerations. 


Of the various methods of selecting optimum ACS circuits and paramenters, which provide top speed with 
maximum accuracy, those based on the complete invariance principle, or invariance up to ¢, chosen to ensure 
best quality for each particular system, have been recently acquiring the greatest prominence. 


The invariance principle, i.e., independence of the controlled quantity from interference,has been worked 
out by Academician N. N. Luzin in conjunction with P. I. Kuznetsov [1-3] and Academician V. S. Kulebakin [4-7]. 


Academician V. S. Kulebakin showed in several papers [4~7] that the combined closed and open loop 
systems provide the greatest dynamic accuracy; he listed a number of their most important properties and estab- 
lished criteria for determining their parameters at which the effect of extemal interference is eliminated (if it is 
not superimposed on the input signal) either completely,or up to an arbitrarily chosen component ¢. 


Later A. G. Ivakhnenko [3] and G. M. Ulanoy [8-9] showed the advantages of combined systems as com~- 
pared with normal ones and proved that the methods developed by several authors [10] for the analysis and syn- 
thesis of these systems represent particular cases of applying the invariance principle. 


In these papers G. M. Ulanov investigated combined systems with one motor. 
Combined systems with two motors driving a single load [10, 11] are, however, of practical interest as well. 


When two or more motors are used, the method of coupling them to the load becomes very important. It 
is obviously impossible to use directional mechanical transmission for these purposes, since it jams when the 
sign of the velocity is changed by the input signal (12). 


This circumstance necessitates the use of reversible mechanical transmission among the various types of 
which, only the conical or cylindrical differentials provide an algebraic summation of the input signals, This 
ensures a smooth control over a wide velocity range [11] and the possibility of building up a combined system. 
In such a case the motor,operating according to the open loop principle, provides the main drive without the as~- 
sistance of the closed loop motor, whose object is to introduce a correction to the output of the open loop motor. 


It is quite impossible to use normal reduction gears and friction drives in combined systems, since, if they 
are used, each drive is assisted in its operation by the other drive, their input signals do not add and the basic 
Principle of the combined ACS operation cannot be implemented. 
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the form: 


Figure 1 shows a typical combined system with two motors driving a common load, 


Let us introduce the notation: Ij; is the load moment of inertia; 1, is the transmission ratio of the differen 
tial and reduction gear of the first motor (of the open loop); 1, is the transmission ratio of the differential and 
reduction gear of the second motor (of the closed loopy, M, is the moment developed by the first motor; I, is the 
moment of inertia of the reduction gear and rotor of the first motor, M, is the moment developed by the second 
motor; I, is the moment of inertia of the reduction gear and rotor of the second motor, 
= Wt), S(t), amd O(t) = + cy. 


When two motors are driving in parallel, a common load their basic dynamic equation can be written in 


Ua t+ M,+ M, 


M,+M, 


0 


0 


by the roots of equation: 


A(D) = 


1 
1 


—Y¥ 22(D)Y 23(D)is 0 


1 
0 
(D) is 
0 


00 
10 
00 
—i—ii 


0 
1 
0 
—1 


M 
+1 


+ hid + 


The differential equation of the closed system in a matrix form is written as: 


01 
00 
10 
—ii 


Assuming that the motors and their transmission ratios are the same, i.e., that iy= i, = i, and My= M,= M 
expression (2) can be rewritten as: 


From this formula it follows that the apparent rise in angular velocity at the output of the system with two 
similar motors driving a common load in parallel does not occur, for the load becomes evenly distributed between 


Fig. 2 


Let us denote by Y;(D) the transfer function of the closed loop control link, by Yy(D) and Y,,(D) the 
transfer functions of the open and closed loop amplifiers, and by Yyg(D) and YD) the transfer functions of the 


¥(t) 


0 


0 


The stability of this system is determined by the roots of the principle determinant of equation (4), t.e., 


= 1 + Yea (D) (D) in = 0. 


(1) 


(2) 


(3) 


(4) 


(5) 
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It will be seen from (5) that the stability of the system only depends on the parameters of the closed loop 
circuit. 


According to V. S. Kulebakin [4~7] the criterion of complete invariance i.e., the condition of zero error 
= $(t)—O(t)=0 at $(t)0, 
is the identity with zero of the minor A * (D) in the expression 
A*(D) (t) = 4 (D)8 (2). (8) 


For the case in point 
A*(D) = 1 + Yu i; = 0, 


whence it is easy to find the transfer function of the open loop control links 
(ti): 
The transfer functions of motors and amplifiers are usually provided. Let, for instance 
= ¥3(D) = (9) 
By substituting the given equation (9) in (8) we obtains 
= —kyD(DTn 1 + 1)(DTy +1), ky = (10) 


Figure 2 shows the schematic of the correcting link by means of which condition (10) can be fulfilled. 


We find the correcting link transfer function by means of the voltage node method. Denoting the voltage 
at nodes 1, 2, and 3 by Uy, Uy, and Ug, respectively, we can write the following equations: 


i i\; 
+ —(Cu+ = 9, 
i 


where » is the gain of the amplifier and D = d/ dt. 
After transformation the link transfer function will have the form: 


Assuming that the gain p of the dc amplifier ts equal to (40-50): 10° and neglecting the term 2 the trans- 
fer function of the link can be written as: # 


Y (D) = —kD{D*T, + DT, +4], (12) 
where 
= 
Con 


(Ria + Res + Res) 


The required values of R and C can be obtained from equation (12) with the knowledge of the system ele- 
ment parameters. 
SUMMARY 


1) A practical application of the complete invariance principle to combined systems with two drives 
Operating a common load is given. 


2) The circuit of a correcting link in an open loop network,providing the application of this principle, is 
supplied. 
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CALCULATION OF STATIC CHARACTERISTICS OF REACTOR 
CONTROL SYSTEMS 


D. A. Alenchikov 


(Moscow) 


Graphical methods for the solution of nonlinear problems, which are encountered 
in the design of complex, reactor control systems are set forth. A graphical method is 
given for plotting the path of the operating point (POP) of a saturable reactor, taking 
into account as a first order approximation, the losses in the saturable reactor, depend- 
ing on the operating current and the voltage drop in the saturable reactor. 


Methods of plotting the output load characteristics are described in detail. For 
this it is proposed to augment the POP by functional relationships between one of the 
output quantities and the specific current in the working windings. As a result, it turns 
out to be possible to determine directly the characteristics of output quantities for fixed 
values of the SR control current. 


A method is proposed for determining initial data for the choice or the design of 
a controller in the case when the law governing the mutual variations of the saturable 
reactors output load quantities is known. 


This work is a more detailed account and a further development of the graphical methods for the calculation 
of the static characteristics of reactor control systems, proposed by V. S. Kulebakin [1). 


A magnetic amplifier, or saturable reactor (SR), is an inductance-resistance, controlled, nonlinear, addi- 
tional impedance, which is connected in series with a load to control its current and voltage. 


Nowdays, a SR is taken to be, in the calculation of characteristics, an inductive impedence, while the inphase 
component of the voltage drop across it is taken to be proportional to the current in the SR operating windings. 


At the same time, a first approximation of the losses, depending on the operating current is taken into ac- 
count. In addition to that, the resistive losses in the SR consist of steel losses APs, losses due to the operating 
current’s first harmonics AP); and higher harmonics DAPj4j, and additional losses in the SR AP. 


i=n 
i=2 


Let us consider briefly the way in which these losses vary, when the SR operates under load. The redistribu- 
tion of voltages between the load and the SR is connected with variations of the fractional cycles of the dynamic 
hysteresis loop [3]. The larger the voltage drop across the SR, the larger is the portion of the dynamic hysteresis 
loop that the fractional cycle of magnetization constitutes. This is accompanied by greater hysteresis and eddy- 
current losses in the steel. 


On the other hand, the redistrubution of voltage between the SR and the load is connected with the variation 
of the magnetic circuit's saturation angle (firing angle) [2, 3]. The larger the voltage drop across the SR, the 
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larger is the saturation angle [3]. An increase in the saturation angle is accompanied by an increase in the form 
factor of the current in the SR operating windings [4]. For equal values of output current of the SR copper losses 
increase at the expense of the higher harmonics FAP). Losses due to higher harmonics depend on the voltage 


drop across the SR and the operating current. Losses due to 
PSR. the fundamental harmonic AP)yjincrease as the operating 


jojo teats current increases. 
oped “fe Pe 2 An increase in the operating current is accompanied 
a ; by increases in the composite magnetizing force and the 
leakage flux. The latter interacts with the current in the 
|" operating windings,displacing it out to the edge of the 
0 Tov aw., conductor. This results in an increase in copper losses, 
Fig. 1. ndence of cos ¢p aw,, for dif- equal approximately to AP,. Extra losses arise as the 
obtained with operating cument ts 
a SR having a nontoroidal magnetic circuit From what has been said it follows that the resistive 
made of transformer steel andhaving small losses in a SR, and, therefore, the inphase component of the 
equivalent gaps in the magnetic circuit. voltage drop across the SR depend on the operating current 


as well as the voltage drop across the SR. The relationship 
between the listed components of the losses [1] depends on 

the material of the magnetic circuit, and the construction and parameters.of the SR. For a toroidal magnetic cir- 
cuit of high quality materials, the value of the losses is affected mainly by variations of the load current form 
factor, while leakage and iron losses have a small effect. For SR with nontoroidal shaped magnetic circuits made 
of low grade types of transformer. steel, and considerable gaps, the manner in which the losses vary is affected 
by the additional losses, while the losses brought about by changes in the current form have less effect because 
the saturation angle of the magnetic circuit is expressed weakly. 


Analysis of experimental data for different SR with 
nontoroidal magnetic circuits of transformer steel shows 
UI that on the basis of losses they can be divided into two 

' & F main groups. The first group includes SR , the magnetic 
circuits of which have small equivalent gaps (E, and O- 
type with a gap in one arm and Il-type with a widened 
it F back). For such SR the ratio of the inphase component of 
%% the voltage drop to the total voltage drop across the SR is 
me r 6,0, observed to be approximately constant (Fig. 1): 


m,-phase 


hase 
4 The constancey of cos is connected with the averag- 
| ing of losses, The true cos ¢gp fluctuates about some mean value, 


The second group includes SR with nontoroidal mag- 
netic circuits and appreciable equivalent gaps. For such 
SR greater accuracy results from averaging the inphase com~ 
ponent of the voltage drop across it for fixed values of cur- 
rent in the operating windings. 


In general, the relationship AU, = f(Isp) is nonlinear. 


u- For SR, which fall in between the first and second group, the 

c losses can be characterized in a composite manner. SR 

; Fig. 2. 1) Power SR; 2)load with output with toroidal shaped magnetic circuits of high grade mag~- 

- quantities B, and B,; 3)controller for SR netic materials are not considered by us, since power SR 
control current; 4) drive; 5) feedback loops with such magnetic circuits are not used as a rule. 

ion involving output and intermediate qan- A SR can be characterized as a nonlinear, complex, 
tities in the reactor control system. 


additional, controlled impedance by a family of volt-ampere 
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characteristics for different constant values of control current. Such a family is usually determined experimentally, 
In order that one family of experimental, volt-ampere characteristics could be used to characterize a large group 
of SR of similar construction, with magnetic circuits of the same material, differing in dimensions and winding 
specifications, it is convenient to express the volt-ampere characteristic in unit quantities. The unit voltage drop 
yu can be expressed in terms of the absolute voltage Usp by means of the relationship 


The unit current in the operating windings aw, is defined in terms of the absolute current LL: 


I. = kyl. (4) 


The unit, resultant, control current awyyp is expressed in terms of the mean values of currents in the control 
windings: 


N=n 
TonWen 


= 2 lew. ©) 


N: =] 


The following notations are introduced in equations (3)~{5): f is the frequency of the applied voltage, in 
cycles per second, S, is the active cross section of magnetic circuit, in cm*, Ww is the number of turns in the 
operating winding or the N-th control winding, 7, is the length of the middle line of force in the magnetic cir- 
cuit, incm, ky, kz, kgy are coefficients of proportionality between the unit and absolute values of, respectively, 
voltage, operating currents, and control currents of the N-th control winding, c,, c, are coefficients, depending 
on the scheme for connecting the operating windings. For opposing~parallel connection of the operating wind- 
irigs cy= 1, cy = 0.5. For opposing-series connection c, = 0.5, c,= 1. 


~ 


Fig. 3. Plotting the path of the operating point of a SR operating in~- 
to complex, nonlinear load. 
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The coefficients ky, ky, and kg are determined by the structural parameters of the SR, the network frequency 
and the connection scheme for the operating windings of the SR. 


Unit volt-ampere characteristics can be obtained by means of electrodynamic instruments. In this case, 
they will represent relationships between equivalent, sinusoidal voltages and currents. 


Figure 2 shows the block schematic of a typical, complex reactor control system, the equivalent circuit 
of one phase of its power section, and a vector diagram of equivalent sinusoids of voltages and currents in the 
equivalent circuit. For any variations in the impedances of the load and the saturable reactor, it is correct to 
express the relationship between the vectors of the equivalent sinusoidal voltages by the equation 


U =U, + Usp. (6) 
From triangle OAB (Figure 2) we find 
= (UL cos + Usp 008 + (YU, sin + Usp sin (1) 
Multiplying both sides of equations (6) and (7) by, respectively, k, and kj, we obtain 
Uc = Uy + (8) 
ug = (uw, Cos + + (uy sing, + Sin (9) 


Let the losses in the SR be characterized by the constancy of cos?cp. The load impedance is completely 
described by its volt-ampere characteristics and the dependence of cos#,, for example, on voltage. The load 
volt-ampere characteristics can be reduced to unit quantities by means of coefficients ky and ky, and can be 
plotted in quadrant IV in Figure 3. It is convenient to choose the aw,, axis to be common to the volt-ampere 
characteristics of the SR and the load. In quadrant ILI of the same graph we can plot the dependence of cosgy, 
on up. The scale for u and uy, should, in all cases, be chosen the same, while that of cose should be arbitrary and 
convenient for reading. As a result, connections exist between the unit volt-ampere characteristics of the SR 
(quadrant I in Figure 3), the load (quadrant IV), and the curve relating cosy, to the unit load voltage (quadrant 
Ill). In order to find the path of SR operating point (abbreviated POP) it is necessary to draw vector diagrams. 
For this a circle is drawn with radius cosy = 1 (dotted line in quadrant II of Figure 3) and another with radius 
Ug = kyUsg. The circle ug is the locus of the end of the network voltage vector expressed in unit quantities. 


Let it be required to find the point on the POP corresponding to point a on the curve uy. The abscissa of the 
point on the POP is determined by the position of point a, since the awa axis is common to u and u,. In order to 
find the ordinate of u point a is projected onto the curve cosg,, (point a*) and a vertical line is drawn to its inter 
section with the circle cosw = 1 at the point a". Then a straight line is drawn through the origin and point a”, 
and a distance Ob, = ab = uy is measured from the origin. From point b, a straight line is drawn at an angle 
sR to the horizontal axis, up to its intersection with the circle u, at the point by. Then the segment byby is 
measured and laid off in quadrant I (segment bb}). 


In order to simplify plotting, it is convenient to lay off the angle gsr once, making use of the auxiliary 
circle cos® = 1, and then to move this straight line parallel to itself. If point b, is connected to the origin , then 
it is easy to see that triangle Ob,b, in Figure 3 corresponds to triangle OCA in Figure 2. Segment Ob, corresponds 
to the vector of the voltage across the laod, segment byb, corresponds tothe vector of the voltage drop across the 
SR, and Ob, corresponds to the network voltage. The only difference is that in Figure 3 voltage components are 
expressed in unit quantities, which is achieved by multiplying all the components of the vector diagram by a 
constantscalar k,. Moving point a along the curve uy, and correspondingly point a’ along the curve cos; and 
repeating analogous plots, we can determine graphically the whole POP (solid curve in quadrant I of Figure 3). 
Up until now the POP has been determined analytically [2] by the solution of an equation analogous to (9), but 
under the assumption that cos?gp = 0. 


No restrictions were imposed on the nature of the load's volt-ampere characteristics, and, consequently, 
the proposed graphical method of plotting POP is applicable to nonlinear loads in which the ratio of the load 
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Pig. 4. Plotting load characteristic of SR operating into a resistive, constant load. 


voltage to the current does not remain constant. Analytical methods of calculating points on the POP are quite 
. complicated, when the path of the operating point is plotted for the case of a nonlinear load impedance [2]. 


In the investigation of a SR one often determines its load characteristic which represents the dependence 
of the SR output current on the control current, when the SR operates into a given load. 


, 
“Os ka 


Fig. 5. Load characteristics of SR operating into 

a symmetrical, wye-connection or resistances, Ry = 
= 2.3 ohm, RL, = 4.65 ohm, Usi = 220 v. Dotted 
curves are theoretical ones with Ipj = 0.485 a, 
solid curves connect experimental points with 
Ipi = 0.4, Roy = 340 ohm, Rp ir = 360 ohm. 
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Let us examine the determination of the load 
characteristic for the case of the three-phase con- 
nection of a SR with a feedback loop, (with self- 
saturation) when it operates into a resistance Ry 
(Figure 4). For the given circuit of connection of 
the operating windings c, s 1, c, s 0.71, and the 
unit volt-ampere characteristics of the SR are of 
the form shown in quadrant I of Figure 4. Theunit 
volt~ampere characteristic of the resistance R, is 
a straight line passing through the origin and a point 
corresponding to an arbitrary L, 0. The abscissa 
of the point is determined from equation (4) and the 
ordinate from equation (3). It is convenient to write 
the latter, in this case, in the form 

uy =k, RT. (10) 

Figure 4 shows the volt-ampere characteristics 
for Ry; < RLz. The construction of the POP is sim- 
plified to some extent because cosq =1 and the 
vector uy, coincide in the direction of the horizontal 
axis. In order to describe a better point a on the curve 
Ria (Figure 4), the plotting of the path of point b¥, is 
shown by menas of arrows. 


ad 


Fig. 6. General view of a three-phase SR, 
having a nominal power rating under load of 
12.5 kw. 


The point of intersection of the POP with the SR 
volt-ampere characteristic determined the unit, resultant 
ampere-turns of the control and operating windings. In 
order to determine the actual values of the control and 
load currents, it suffices to divide the obtained values of 
and by, respectively, ks and ky. By determin- 
ing the points of intersection of the POP with other SR 
volt-ampere characteristics one obtains the load character- 
istic of the SR (Figure 5). The plotting carried out in 
Figure 4, and the characteristics shown in Figure 5 relate 
to a practical SR (Figure 6), connected in the circuit shown 
in Figure 4. 


The magnetic circuits of three-phase SR (six mag~ 
netic circuits) are made of pi-type laminations with a 
widened back edge. The width of the lamina bar is 40 mm, 
the lamina material is type £-310 heat treated steel. The 
magnetic circuit's cross-section S, = 23.4 The load 
winding has 162 turns, the control winding has 540 turns, 
and the bias winding has 200 turns. 


For the given SR,cosgsR is constant (see Figure 1, 
upper family of points). 


& & 
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Fig. 7. Method of plotting the family of characteristics B, = f (B,) of a 
closed-loop control system for the case when for fixedaw., uj * Const. 


Figure 5 shows that the coincidence between calculated and experimental curves is quite satisfactory. We 
note that similar accuracy can be achieved in the calculations of curves without taking into account all the losses 
in the SR, but, in this case, one is obliged to use empirical correction formulas ([2], page 293). A shortcoming 
of empirical formulas is the fact that they are related to the SR load characteristic, to a particular circuit of 
connections for the SR operating windings, and to the type of load. 
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If the SR output (load current) uniquely determines the load's output quantities (thermal, luminous, etc.), 
then the load characteristic can serve as a basis for determining the dependence of output quantities on the SR 
control current which are of interest. 


When a SR operates into a motorload, the output current of the SR uniquely determines, at best, only one 
output quantity cf the load (magnetic moment of a direct current motor with semiconductor rectifiers), When the 
SR operates with an asynchronous motor, the output current of the SR does not determine uniquely a single output 
quantity (slipping torque). At the same time, for a motor load it is necessary to know the functional dependence 
of one quantity (B,) on another (B,) in the presence of various methods of variation of the SR control current. 


Let us examine only two of the more prevalent cases of plotting the output characteristics of a reactor con- 
trol system. In the first case it is required to plot a family of curves, establishing the dependence of one output 
quantity (By) on another (B,) for different, constant values of the SR control current (the, so-called, open-loop 
characteristics of a control system). 


A motor load is characterized by the fact that its equivalent impedance (ratio of the voltage at the output 
of the SR to the current) does not remain constant but depends not only on the current, but also on the rpm of the 
motor. In this case, the equivalent impedance of the load can be called nonlinear and variable. Such an im- 
pedance can be characterized by a field (or family) of volt-ampere characteristics. The dependence of the volt- 
ampere characteristics of the load and cos@ on output quantities can be determined analytically or experimentally 
[5]. In order to plot a family of characteristics for an open-loop control system, it is more convenient to plot in 
quadrant IV a family of unit, volt-ampere, load characteristics, and in quadrant 
Il, a family of curves relating cosy, to u, for several fixed values of one of the 
mmo 7? 3 4 output quantities, for example, B,. A family of curves establishing the de- 
pendence of the other output quantity (B,) on the unit current in the operating 
windings (Figure 7), for the same values of B,, is plotted in quadrant IV. The 
scale of B, can be arbitrary, while the scale of aw, must be the same as for the 
unit volt-ampere characteristics of the SR. The POP is plotted analogously 
from the volt~ampere characteristics and the curves relating cos¢1, to u,. 
0 Since earlier on (Figures 3, 4) the method of plotting the POP was illustrated 
for the case when cos¢.,, » const, in Figure 7 curves are plotted for the case 
= SPs when for every value of conch the inphase component of the voltage drop is 
ead aa approximately constant. In this case, quadrant I (Figure 7) should contain 
p control system, plot- 
eal ihieik tie Matnater~ curves of the quadrature component of the unit voltage drop across the SR as 
istics in Figure 7 a function of aw~. In connection with this characteristic of the SR, the in- 
; phase component of the voltage drop across the SR can be included in the 
load yolt-ampere characteristics, which is achieved by adding vectorially the 
unit, inphase component of the yoltage drop across the SRand the unit load voltage: 


= Uj. (11) 


Here u,’ is the vector of the unit load voltage with the inphase component of the voltage drop across the 
SR dy taken into account. 


Cosg, is converted at the same time. In quadrants III and IV of Figure 7, corresponding quantities with the 
inphase component taken into account are shown as dash lines. Obviously segment byb, (Figure 7) is perpendicular 
to the horizontal axis, since in this case the quadrature component of the voltage drop across the SR is being 
determined. 


Near the POP curves in Figure 7 are marked fixed values of the quantity B,, for each one of which the cor- 
responding unit volt~ampere characteristic u;, is plotted. The point of intersection of the POP with any unit 
volt-ampere characteristic of theSR uniquely defines the unit value of the resultant SR control current awyp, the unit 
current in the SR operating windings (load current) aw _— and the output quantities By and B,, because the cur~- 
rents through the SR and load are common and the family of POP is plotted for the same values of B, as the family 
of curves relating to B, to aw,,. If one passes through the points of intersection of the SR volt~ampere character~ 
istic, for a given value of awyp, with different POP and determines the quantities By, and B, corresponding to 
these points, one can plot one characteristic B, = f (B,) for the chosen values of awyyp. Choosing other values of 
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he 


awyp one can plot a family of characteristics for an open-loop control system (Figure 8). If the scales of B, are 
the same in Figures 7 and 8, then the construction of the characteristics in Figure 8 reduces to measuring values 
of B, in Figure 7, and transferring them to Figure 8 for the values of B;, corresponding to the point of intersection. 


Fig. 9, Plot of the control characteristics awysp = f(B,). 


Fig. 10. Given character- 
istics of By= f (B,) and calcu- 
late ‘, from Fig. 9, control 
characteristics awtyp™= f (By) 


Let us examine the second case where the law governing the mutual 
variation of the output quantities is known, and it is required to find the 
limits within which the given law is obeyed, and to determine the initial 
data for the choice or design of a controller. In this case it is necessary 
to determine how the SR control current depends on one of the output 
quantities. 


In contrast tothe case just discussed, the unit, volt-ampere character~ 
istics of the load and the curves of B, as a function of aw~ are plotted ac- 
cording to the law governing variations of the output quantities, making 
use of the field of load characteristics. Figure 9 shows curves of uj, Ccos?,, 
and B, for three given laws governing the variation of B, as a function of 


By (P Pa» 


The point of intersection of the POP with the unit volt-ampere 
characteristic of the SR determines aw~, awyp, B, and p, which makes 
it possible to plot aw;yp as a function of B, for a fixed value of p. 
Obviously, for this it is necessary to find the points of intersection of the 
given POP with different SR volt-ampere characteristics. Calculations 
yielded the result that all the POP ps lies in the region of residual voltage 
drops across the SR, and the required characteristic of p, (Figure 10) 
cannot be ensured for any variations of the control current. 


The characteristics of a reactor control system were plotted without use of the SR load characteristic. 


The POP determines the power efficiency coefficients and the power factor of a system only in particular 
cases when cos, = 1[2]. In general, when cosy, # 1, the system's power factor is determined graphically with 
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the help of the POP, the load volt-ampere characteristic, and the relationship between cos, and uy (or cos? 
” ti). For this a complete vector diagram is plotted in practice. In order to lay off the cosine of the angle 

¢ formed by the vector of the unit line voltage and the horizontal axis (Figure 9), one can make use of the 
aes circle of radius cos? = 1. The graphical construction of cosy, is shown in Figure 9 by means of arrows 
for point A. Knowing cos?, one can find the efficiency of the system if the output quantities determine the 
useful power Poy. 


The efficiency of the power circuit is 


Pouk, 
mUsf aw_ cos@,” 


(12) 


where my, is the number of phases, Ugg is the phase value of the m,-phase eas sie (when m, = 1, it is the 
supply voltage). 


Taking into account the power required for operation Pro 


Pou 


mille % + Fro (13) 


In the case when one part of the SR losses is characterized by the inphase component of the voltage drop 
across the SR being constant for fixed values of aw,,, while the other is characterized by cos¢_» being constant, 
the constant component u; is included in the load volt-ampere characteristics while the variable one is taken into 
account by a corresponding inclination of the line byb, (Figure 7). This case is actually an elementary combination 
of the two cases considered above and is not considered here. 


CONCLUSION 


1. The proposed graphical method for the calculation of the static characteristics of SR enables one to plot 
the path of the operating point of the SR taking into account, as a first approximation, the resistive losses in the 
SR, which depend on the operating current and the voltage drop across the SR. The method’s accuracy is sufficient 
for practical purposes. 

2. The method of calculation does not depend on the nature of the nonlinearity of the load resistances and 
can be used in the case when the load's volt-ampere characteristic cannot be approximated exactly mathematic~ 

ally. 


3. The proposed method for plotting directly the characteristics of the load*s output quantities by adding 
to the load yolt-ampere characteristics the functional relationship between one of the load's output quantities 
and the unit current, increases the potentialities and the significance of the path of the operating point, and de~ 
creases the difficulties of computing the characteristics of a reactor control system for open-loop control. 


4. The proposed method for determining directly the initial data for the design (or selection) of a regulator 
for the SR control current is also based on the use of the POP and the functional relationship between one of the 
output quantities and the unit control current for a given law governing the mutual variation of the output quan~ 
tities. 
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CALCULATION OF AC MAGNETIC REACTOR AMPLIFIERS 
WITH INTERNAL FEEDBACK 


N. A. Kaluzhnikov 


( Kharkoy) 


The most typical circuits of magnetic reactor amplifiers with internal feedback 
are considered. The circuits are analyzed in two cases which are extreme from the point 
of view of nonsinusoidal distortion. Possible calculation errors are estimated. The re- 
sults of the analysis are given in the form of design curves. 


General Assumption 


An important part of the problem of the design of magnetic amplifiers (MA) is the calculation of the ac 
circuit. The essence of this consists of the determination from given load parametes (resitance Ry, inductance 
L,» maximum current I; 4x" maximum voltage Uy max) Of 4 series of unknowns; U, (supply voltage), U, may 
(maximum voltage on the ac windings of the MA), ug.max (amplitude of the reverse voltage on the rectifier). 


The quantities, which are determined from the calculations are in turn used in the design of the supply trans- 
former, rectifiers, and, finally, the magnetic amplifier itself. 


The quantity ud.max represents the maximum, instantaneous value of the diode reverse voltage, which it 
is necessary to know to select the appropriate rectifier. The quantities U, and Uy may can be expressed in either 
effective or half-period mean values. In the discussion that follows, mean values are used. 


The volume of the core steel Vs in the magnetic amplifier can be calculated if the, so-called, size power 
of the MA is known: Py max = Ur.max ‘ly max* 


P r.max 


1 
st = mas ()) 
Here f is the frequency of the supply voltage, Bmax 4Md Hma, are the maximum values of the induction 
and the intensity of the alternating field (see below). 


The number of turns in the ac winding of the magnetic amplifier is determined from the well-known 
formula 


Ur, max 
~~ 4fBmax5 st (2) 
where S,, is the core package cross-section. é 
The quantity U, is used in an analogous manner in the design of the supply transformer. 


In this way Us, Uy max 20d Ug max ate necessary intermediate parameters in the design of magnetic am- 
plifiers, and their determination from a given value of Uy max is, in fact, the main problem in the calculation 
of the ac circuit of a MA. 


u 
The ratios k, = Us/ Uy, = Ur.max/ UL. max: 294 kg = Ug max/"s max > ¢ where Us max 
s 


is the amplitude of the supply voltage) are in general, variable quantities depending on the chosen mode of 
operation of the MA magnetic circuit, its circuit details and the type of load. Their determination is the object 
of this article. Moreover, in certain cases the relationship between the currents I, , flowing through the MA 
windings, and the current through its load I, is nonlinear. Among the problems considered further on is that of 
determining the relationships kj = I; / . 


Much has been written about the design of magnetic circuits [1, 6, 4]. Without stopping to deal with this 
problem in detail, let us limit ourselves to a few general comments, 


In most practical design cases it is not necessary to calculate the whole control characteristic of the mag- 
netic amplifier 1, = f | lin], and it suffices to determine only its limiting values. Two points on the magneti- 
zation curve of the MA B= f | H.. H| correspond to these limiting values. One of them is usually chosen near 
the bend in the main magnetization curve B= f | H,,| where H= = 0, and is described by the coordinates Bmax 
Hwmin: Hz = 0. The other point corresponds to the maximum value of the magnetizing field and is determined 
by the coordinates Brin» Homax: H=max- 


The relative position of these two points can be expressed by the following dimensionless parameter 


Kp—_max _ y — 4~max _ 


Bmin Almin min” 


The quantities k,, kp, kg and kj are functions of kg and k,,. However, for large values of kj; (of the order 
of 7-10 or higher) the dependence of the parameter on ky can be practically ignored. In the subsequent dis~ 
cussion ky is taken to be oo, and the obtained results are applicable only for high values of kj;. Since magnetic 
amplifiers with cores of iron-nickel alloys and the best types of electrotechnical steels operate usually with 
high values of ky;, this restriction on the application of the obtained results is not particularly severe. With 

ky = ©, Uy max = Use Ke = ks, and kp = Us/U, 144+ Which is what will be assumed in the subsequent discussion. 


Besides this assumption, we also assume that all the diodes are ideal (R¢.= 0, Rpey = 0), the dissipation and 
resistance of the MA ac windings are zero, and there are nocore losses. 


However, even with all these assumptions, the calculation difficulties created in taking into account a 
practical magnetization curve would be practically insurmountable in the analysis of even relatively simple MA 
circuits, while the value of the obtained results would nevertheless be quite conditional in view of the appreciable 
scattering of characteristics, which is inherent in magnetic materials. Consequently, the more rational methods 
of calculation seem to us to be those which enable one to obtain the simplest solutions, but which also enable 
one to determine the value of the error, which is made in the calculation. 


In the light of what has been said, two extreme cases are considered in this work. 
1. The case of an *ide71" MA having a core with a rectangular magnetization curve. 
2. The case ofa “linear"MAin whichthe function B = f (H) is linear. 


The theory of the “ideal” magnetic amplifier is used widely in the literature [1, 2] in describing the physics 
of processes in MA. It is used much less frequently in calculations. 


The assumption of a “linear” magnetization curve at first glance contradicts the principle of operation of 
a MA, which is a typically nonlinear device. However, it is used widely in calculations (1, 4], when the currents 
and voltages of the magnetic amplifier are assumed to be sinusoidal. Moreover, when the MA operates with a 
strongly saturated magnetic circuit, the so called, residual inductances, which by their nature are close to being 
linear, are of primary significance. Exactly such regimes are required when it is desired to design a MA with a 
small volume of steel (see, for example, [4]) and, consequently, the concept of a “linear*MA receiv.s a definite 
physical meaning. 


The "ideal* and “linear* cases are extreme from the point of view of possible nonsinusoidal distortion in 
the ac circuit of the MA. Hence, having compared the results for one and the other case, one can choose the 
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Fig. 1. Circuits of magnetic amplifiers with internal feedback 
loops. 


worse one, and at the same time, evaluate the maximum possible error. Analyses of thistype were carried out 
in [8] for magnetic amplifiers operating with a single phase bridge rectifier. 


In this paper the same means are used to analyze a row of typical circuits of MA with internal feedback 
loops (Figure 1). Existing papers have actually not touched on the question of the characteristics of the design 
of such circuits, while the results of the analysis of MA with an external feedback loop are mechanically trans- 
ferred to MA with an internal feedback loop. Such an assumption is based on the well known similarity of the 
indicated circuits, but a complete analogy can be established only for the case of the operation of an “ideal” 
MA into a resistive load. Considerable divergence in the properties of various modifications of circuits with a 
feedback loop exists in practical MA [5]. This situation forced the author to devote a special investigation to 
magnetic amplifiers with an internal feedback loop, with only two types of load on the MA being considered: 
resistive, and inductive -resistive, with a large value of 


(= 


The physics of processes in ideal magnetic amplifiers with an internal feedback loop are described in de~ 
tail in the literature [1, 2]. The operation of rectifying circuits with current-limiting, linear inductances is also 
known sufficiently well [3] Hence in the subsequent discussion questions of the physics of processes are analyzed 
from that point of view which is necessary for determining the design coefficients ks, kg, and k; 48 functions of 
kp- 

Half-Wave Circuit (Figure 1, a) 

a. Resistive Load 


Ideal case. The physics of processes in such a circuit are examined in detail, for example in[1]. The 
main conclusions of the analysis carried out there can be reduced to the following. If the applied voltage obeys 


the law u, = Us max Coswt, then during the conducting half-period in the interval ut=(--£) (p is the 
saturation angle) all the input voltage is applied to the ac winding of the MA, at the same time no further cur- 
rent flows and uw, = 0. Further when w= y= the core becomes completely saturated and the whole input 
voltage is applied across the load. During the nonconducting half-period in the interval wt = —-(n + ¥) 
the reactor again takes up the whole input voltage, at the same time returning to the source the energy stored 


ty during the first half-period, while when wt = (4 + #2 " 
Be the diode blocks the remainder of the input voltage. 
To sum up, we can write relationships for the mean 
values of U, = U,+ Uy + Ug, where Ug is the diode reverse 
a voltage. Moreover U; = Ug = Uy, Hence, Us= 
06 4 = Ur+ 2U,. Since U, = U,/ kg, it follows from here that 
k 
Fig. 2. Characteristics kg = f (kp) for an 
*ideal” magnetic amplifier. = 
Insofar as the current passing through the MA and the load a 
is one and same, then for the circuit inFigure 1, a, we obtain I. 
in all cases d 
I q 
r, max, 
P Let us now find the quantity kg. Obviously, the maximum 


reverse voltage, which is possible in the circuit in Figure 1,2, q 
is equal to the amplitude of the supply voltage us max. How- 7 


? ever, the amplitude of the input voltage is applied to the diode i 
\ only when the saturation angle is — 0. In the re~ 
maining case (when » = 0 —~ w/ 2) the diode reverse voltage q 
\ \ will be less and will be determined by the expression : 
Ud.max = Us,max©0s (5) q 
2 
Ne To plot the curve ky = Ug. max/ Ys, max 2Xug,max/U,)= 
SS = f (kp) it is necessary first of all to express U, as a function 1 
= — of the angle p: a 
Fig. 3. Characteristics ky = f (kg) for a 

“linear” MA with a resistive load. 0, n4f Us max 008 wt dwt = U, a » (8) 

1) For the circuit in Figure 1, a; 2) for —th 

the circuit in Figure 1, b; 3) for the cir- g 
cuit in Figure 1, c. whence : 
ke 

Be T¥aing 
From equation (5) we obtain 7 

“s.max “s.max (8) 


Solving equations (‘7) and (8) together we determine the sought-after relationship k, = f (kg) in the interval 
kg = 1——2. When kge 2 kg = 1 (Figure 2). 7 


Linear case. This problem is described in detail in[3}. Only the main results of the given analysis are 
presented below. 


The instantaneous value of the current is determined by solving the differential equation 


di 
Us mak 008 wt = + x, Toi (9) 


| 
| 
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with the initial condition i;(—#*/2)= 0 


(ot + singe * |, (10) 


where T = arc tg y = X,/R,- 
Mean values of the current I, and voltage Ur will be the solutions of the integrals: 


8 
i d 
iydot, U, =+ x, dot. 
After a series of transformations we obtain 


+sing 

OU, (12) 
(14) 


(11) 


8 T+sing 
2 

(15) 


The angle ¥ determines the moment at which the diode cuts off, and is found from the condition i, (¥)= 


The angle § correspondstothe moment when the instantaneous value of the voltage across the reactor u, 
passes through zero, and can be found from the expression di, / dwt = 0. 


Finally we obtain the following system of differential equations determining 


cos(>—)+singe * =0, (16) 
Solving (14)-(17) jointly, we determine k, = f (kg) (Fig. 3). 
In order to find ug max it is necessary to analyze the curve of the input voltage us max coswt for an 
extremum in the interval between » and (3/2). This permits one to calculated kg = F (kp) (Fig. 4). 
b. Inductive-resistive Load 


The case of an inductive-resistive load with a large value of T,, can be reduced without appreciable loss of 
accuracy in the calculations to a purely inductive load. 


Ideal Case. The nature of the process under these conditions can be obtained on the basis of the theory 
of the “ideal* MA and is illustrated in Fig. 5. 


Let usgive a few explanations; in the interval wt=2/2- magnetization of the core takes place 
and the current is equal to zero; in the interval wt= %¢—~ #/2, the core is saturated and the magnetiza~ 
tion of the load by the current i, takes place; in the interval wt= 4 /2—'(x — p), the load returnsthe stored en- 
ergy to the source,while the core is saturated; in the interval wt = (x ~~) —~(3/2)a, the demagnetization of the 
core takes place and the current is equal to zero. 


Finally we obtain Us = Up + U, and 
k 
= 2 = (18) 
“Pp ae In the ideal case, there 1s no voltage across the diode 
o . (kg = 0). Actually there is a small reverse voltage resulting 
/ from the existence of resistances in the circuit, which were 
taken to be zero to simplify the analysis. 
/ To determine the instamttaneous value of the current, it 
7 is necessary to solve the equation U, may Coswt = wl; (diy, Mwt 
od with the initial condition i,(p)= 0. This gives 
‘ d 
Fig. 4. Characteristics k, = f (kg) for toa (sin wt — sin $) (19) 
a"linear” MA with a resistive load. 
1) For circuit in Fig. 1, a; 2) for the and 
circuit in Fig. 1, b (diodes 1, I); i, dot = 
3) for the circuit in Fig. 1, b (diodes lL" 
11, IV); 4) for the circuit in Fig.1; c. 
cos — (w/2— 9) sing 
U, | 
For the mean values of the voltages across the load and 
reactor we get 
x-2 
? 
\ 
4 (Comparing equations (20) and (21) we see that Xeq: * 
= U_/ |, will be a variable quantity depending on #. Let us 
introduce the concept of a coefficient of nonsinusoidal dis- 
tortion 
Fig. 5. Processes in an “ideal” MA OL i—sin » 
with an inductive load connected ac~ "Fol, 
cording to Fig. 1, a. 
From equation (22) we obtain 
2 
4) 


Solving simultaneously (23) and (24) we determine Ax = f (kg) (Fig. 6). Now we will have 
U, = (25) 


We can also obtain directly the curve of the function 


a 

q 

ii 


ia 
4 
| 
| 
a 
q 
of q 
y 
the 
q 
¥ 
1 


shown in Fig. 17. 


Linear Case. The circuit consisting of a diode and two series-connected, linear inductances, to which the 


given case leads, is an elementary one. For it U, = U;+ Uj, and k, = kp kg 1) while the load current is defined 
by the expressions 


This gives Ay = Uy /IjjwL, a constant value equal to 2/m. In other words, because of the presence of a con- 
stant component of current the inductive reactance decreases. For kg we obtain the expression 


Us 2 kp a 
(29) 


ug here is also equal to zero. 


— 


§ 
== 


\ 
Fig. 6. Calculated coefficient 2, = f (kp) 
for an *ideal® MA, having the circuit 2 ERs mon 
shown in Fig. 1, a, with 2, = why, 
UL Fig. 7. Characteristics kj= f (kg) for an 

1, = kp= *ideal* MA with an inductive load. 


1)For the circuit in Fig. 1, a; 2) for the 
circuit in Fig. 1, d. 


Bridge Circuit (Fig. 1, b) 
a. Resistive Load 


Ideal Case. A detailed analysis of this is given in[1, 2]. The main results are the following: 


U,=U, +U,, U4g=U,, 2h. (30) 
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and 
#3 20) 


As a result we find 


U, kp 
ky = 0.5. (31) 4 
The curve kg = f (kp) for the circuit in Fig. 1, b, has the same form as for the analogous case of the cir- di 
cuit in Fig. 1, a. 4 
Linear Case. A feature of the given circuit is the fact that flowing through the load is the maximum of the @ 
currents flowing through the cores [5]. The nature of the processes is shown in Fig. 8. We will give a few ex- A 
planations. 
In the interval wt = (~"/2)— » (switching) diodes 
Yr I, 0, and Il are open. The current flows through diodes I 4 
and Il, thus bypassing the load. Its value is determined from { 
s.maxC0S = Ar 
with the initial condition if,(~=/2)= 0. Current if, passes 
through the load and diodes II and II, so that flowing through 
4 ty II is difference of the currents — . 
In the interval wt= + #) (work interval for core 
oe ’ ‘wt I) at the moment wt = », the current iz, becomes equal to 
Ver and diode II closes. Now current flows through the 
4, load (if; = {,,) and diodes I and IV. The value of if, is de- | 
m Yen termined from the equation us, max Coswt = + 
Wh — / dwt), At the same time current if; falls to zero. 
In the interval wt = ( + (switching) when 
wt = © + » current if,, which began to flow when wt = «/2, 
becomes equal to if;, after which diode Il opens. The current 
Fig. 8. Processes in a “linear* MA with fF, is defined by equation us, max Coswt = Xp(dir Alwt) and 
flows through diodes I and I, bypassing the load, When ' 
0 eee OP Pi wt = @both the current it becomes equal to zero and diode I 
closes. 


In the interval wt = @~%/, # (working interval for core II) the current iy, is equal to zero. The current 
flowing through the load is ify. 


An analogous process, shifted by 180°, takes place in core II. 
Let us determine expressions for the current i,; successively for all the listed intervals: 


i, = sinot + A = (sin ot + 4), (32) 


The constant of integration A is found from the condition f4(—# /2)= 0: 


wt 
it =i, (wt — 9) + Be 
_ sine 


0) = [con ot ot —9) 


i 
| 
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Here = tg y = X;/R,, while Bis determined from the condition = + 
The equality = enables one to find the angle 


+ sing) = cos (wt — 9) — | 


4 sin = — sin (p— 


(34) 
= sin ot + C = (sin ot + 1 + 2sing), as) 
where C is found from the relationship #%(" + ¥) = %¢¥). 
The angle @ where diode I cuts off is determined from the condition #%,(@) = 0: 
sin § = — (i + 2sin9). (36) 


Now we can determine the mean values of U;, |, U;, and I, : 


on. dwt = 


The second term in this expression becomes zero because iid +p) = Mr,(W). Finally we will have 


U, = —U, sin 9, (37) 
Us 

— (39) 


In determining U,, it should be noted that its mean value over a period must be equal to zero. Hence, 
Uy can be found in the following manner: 


8 

4 d i di 

—n/2 


+i O—in WI = 


Here, the fact that = while if,( ~ 1/2) = 0 was taken into account. 
The angle & is calculated from the expression di,,/dwt = 0, which yields 


sin — 9) = 


Making use of the obtained equation we finally get 
Ur 


(41) 


(42) 


hz = 
— cos wt + + — cos ot +4 + = 
= E+ + 0 — 0 — 2008 4. 


In contrast to the ideal case, here inthe intervals 3/2 and correspondingly (9 ~ — x/2 only diodes - 
I and II have the input voltage across them, while the reverse voltage on diodes II and IV in the intervals , 
¥ —( + p) and correspondingly (x + ~)—~(2e + #), is equal to the voltage across the load. The maximum q 
reverse voltage will occur at § and 6 + w and is equal to uj(5)= f4(5)R,. i 


The results of this analysis make it possible to determine ky = f (Ig) (Fig. 3), kg = f (kg) (Fig. 4), and ky = j 
= f (kg) (Fig. 9) 


b. Inductive-resistive Load 


Ideal case. The results of the analysis carried out by Storm [2] can be reduced to the following relation- 
ships: 


2 
kg = (44) 
(45) 


The simultaneous solution of equations (44) and (45) permits one to obtain ky = f (kg) (Fig. 9) 
For kg we obtain the following expression: 


k 
(46) 


. The diodes in the circuit are under the action of the 

reverse voltage, equal to the input voltage in the intervals 
wt = — and correspondingly wt = (« + 
The nature of the relationship ky = f (kg) is the same as that 
for a resistive load in the case of an “ideal” MA. 


= Linear Case. The process in this case, illustrated in 
Fig. 10, follows the following course of events. 

In the interval wt = (~2x/2)—~ » (switching) core I is 

. % magnetized through diodes I and II, the current being defined 

Fig. 9. Characteristics kj = f (kp) for by the equation us maxCoswt = X_(diry /dwt) with the initial 

the circuit in Fig. 1, b. condition 0: (Us max/ Xr) (1+ sinwt). The 

1) For an *ideal* MA with an induc~- load current is shorted through diodes II and IV. 

tive load; 2) for a “linear® MA with a 

a teed. itbecomes equal to and diode 
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and 
q 
/ 
1) 
8) 
| 
— 
40) 4 
_| 


In the interval wt = x/2 (working interval for core 1) f= I, and the whole input voltage {s applied to the 


load (since diy/dwt = 0). 


Fig. 10. Processes in a “linear*® mag- 
netic amplifier with an inductive load, 
connected as shown in Fig. 1, b. 


In the interval wt = 1/2 © (switching) core I is de- 
magnetized, with current ff, defined by equation ug maxcoswt = 
= (digg /duwt) with the initial condition 4/2) = 1, = fy 
1 ™ (Us max / Xz) (sinwt+ sin). The load is once again 
shorted through diodes II and IV. 


In the interval wt = © 3/2 x current i, = 0, and the 
input voltage is blocked by diode I. The angle @ is determined 
from the condition f,(9)= 0, which corresponds to 9 = © + yp. 
The same process, but 180° out of phase, takes place in core II. 


As aresult kg = kp /\kp ~ 1J while the relationships have 
the same form as in the case of the “ideal*® MA. 


Let us determine the mean value of the current Jy: 


—n/2 


Since (ls max/ Xr)(1+ sin~) = we finally get I, = 
= I, /2 and kj = 0.5. 


Circuit with Zero Ouput (Fig. 1, c) 
a, Resistive Load 
Ideal Case. The results of the analysis, which is well 


known from the literature [1,2] reduce to the fact that the processes in the circuits in Fig. 1, b and c coincide in 
this case, with the exception of the fact that in the second case ug max Will be twice as large. Therefore, if we 
denote kd = ug max/2Us.max, then to determine kg as a function of k, we can make use of the curve in Fig. 2. 
ky for the circuit in Fig. 1, c is equal to 0.5 in all cases. 


Linear Case. The process in such a circuit for a linear inductance has been examined in detail in[3]. Hence 


below are presented only the final results. 


The quantity kg can be determined from the expression 


2 
+ tc0s (48) 


where angle » is determined from the simultaneous solution of the two equations: 


cos = — 


= sing + 1+ 


2 (8—9%)> 
|. (50) 


= (61) 


tN, — 


Here the angle & is calculated from the following expression: q 


=" 


The quantity kg is defined by the equation , 


“d.max 
= cose — sin 


ka = (53) 


$.max 


To determine the angle ¢€ we must solve the equation 
sins = cos sin (p —e) + cos (0—g)e — (54) 


The results of the analysis are presented as curves of k, = f (kg) (Fig. 8) and kg = f (kg) (Fig. 4). 


b._Inductiye-resistive Load. 


Ideal Case. The resultsof the analysis conducted by Storm [2] can be reduced to the following: the relation- 
ships ks = f (kp) and kg = f (kg) and have the same form as for an “ideal” MA with a resistive load, 


Linear Case. The process in such a circuit in the case of a linear inductance is examined in detail in [3]. 
The results of the analysis show a complete analogy in the form of the relationships ks = f (kg) and kq = (kp) 
for the case of the “ideal® and “linear” MA. 


Full-Wave Circuit with AC Output (Fig. 1, d) 
a. Resistive Load 
Ideal Case. The nature of the process is well known from the literature [1]. 


, For an ideal MA connected in a circuit, as shown in Fig. 1, d, with a resistive load, we can write k, = 
= kp/ikp — = 0.5. 


The value of kg for the circuit in Fig. 1, d is equal to zero in all the cases being considered. 
Linear Case. Here the process can be described by the following system of differential equations: 


di 

ug = + Xr jo 
ing tire = iy (57) 


Adding (55) and (56), taking (5'7) into account, we get 
xX, di 
us = : (58) 
Substracting equation (56) from (55) we obtain 


1) 


dig 
(59) 
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°° °° 
= 


This is the equation of the mesh current which flows 


through the ac windings of the MA and, which has a constant 
value if the resistances in the mesh are neglected. 
ANN Finally we obtain 
Nu 
i, cospcos(wt— 4), (60) 
Fig. 11. Relationship Igy = f (kp) for an 
*"ideal*® MA with a resistive 1 and an in- + cos (wt— ¢)]}. (61) 
ductive 2 load. 
Here = arc tg(X;/ 
The mean valuesofl, and ly take on the following form: 
2 
(62) 
Tr = (63) 
Now we can calculate kj: 
n 
(64) 


To determine k, we must solve the equation U3 = U?2+.U*,, which is analogous to the case of a MA with 
‘no feedback loop. As a result we find 


65 
Vani (65) 

b. Inductive Load 
Ideal Case. The nature of the process is analogous to the case of MA circuit with no feedback loop, dis~- 


cussed in [7], with the one difference tha here kj= 0.5. The curve of k= Us/wLil, = f (kp) is shown in Fig. 7. 


Linear Case. The operation of the circuit can be described by the system of equations (55)-(57), if we re= 
place in them i,R, by X,(di,/dwt). As a result of the analysis one can come to the conclusion that the values 
of kg as expressed for MA without a feedback loop and with an internal feedback loop also coincide in this case, 
and k, is expressed by the equation k, = kp /fics~ 1) white ky = 0/4. 


In concluding, let us present curves of 8 = ty ef /l, ay=f (kg) for an "ideal" MA since in the case of an ac 
output, it is the effective value which is of interest (Fig. 11). 


CONCLUSION 


In order to simplify the use of the results of this paper all the main conclusions are tabulated. 


The values of ks, kg, and k; calculated for *ideal* and "linear"MA often diverge quite appreciably. Hence 
we should further clarify the question of the final choice of these design quantities. 


*Ideal® MA as a rule yield large values of ks. Naturally, in this case Py max 20d Vs_ will be larger, for the 
same load power, than in a linear MA; moreover, the design value of the control power will be somewhat higher 
than the actual power necessary. - In other words, the MA turns out to be slightly too big, but reliable in control. 
If the dimensions are not too important, then it is recommended that the values of kg obtained for an “ideal” MA be 
used. 


5) 


Formula or curve { 
Type of a 
Circuit Type of load approxima- g 
tion ka ky 
q 
Fig. la | Resistive Ideal 2 kB/(kp—1) Fig. 2 i q 
Linear Fig. 3 Pig. 4 i gq 
Inductive Ideal Fig. 
Linear kp/(kp—1) 0 
Fig. 1b . | Resistive Ideal kp/(kp—A) Fig. 2| 0.5 
Linear Fig, 3 Fig. 4| Fig. 9 q 
Inductive-resist- | Ideal Pig, 2 Fig. 9 
ive kp/(kp—1 Fig. 2 5 
Fig. Resisti e kn/(kp—4 Fig. 2| 0.5 
ig. 1¢ ve Linear 4| 0.5 
Inductive-re- Ideal kp/(kp—1) Fig. 2 0.5 q 
sistive Linear kp/(kp—1 Fig, 2 0.5 
Pig. 1d Resistive Ideal kp/(kp—1i 0.5 
Linear kp/V 0 n/4 
Inductive ideal Fig. 7 (ks’) 0 0.5 q 
Linear kB/(kp—1) 0 n/4 


The design based on the “linear” case yields slightly low values of steel and control power, which makes 7 
it necessary to increase by about 10-20% the control ampere~turns obtained from the design [1). q 


Other recommendations can be given on the use of the obtained results. 


Thus, in the design of cores with a low volume of steel [4] one can make use of the value of k, for a *linear® 
MA. As a matter of fact, in this type of design high values of kp are usually chosen when the deviation between 
the “ideal” and “linear cases” are small. For reliability it is useful to increase kg by 10-15% above the values of 
the range of deviations Ak, = kiq ~ ks. lin. 


For high quality steels and iron-nickel alloys of medium quality one can usually use approximately mean 
values of kg, while for the best materials with rectangular hystersis curves, one must use values of k, for the 
“ideal” MA, allowing a reduction in them of not more than 15-25% below Ak,. 


It is quite obvious, of course, that these recommendations are quite approximate, since beside the quality 
of the material, the nature of the nonsinusoidal distortion is strongly affected by the chosen mode of operation 
of the magnetic circuit. Thus, for small multiplicities of the field strength or high multiplicities of the induction, 
the MA approaches, in its properties,a “linear* one. Approximately the same reasoning can be applied to the 
choice of k;. 
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CONTROL CHARACTERISTICS OF THREE -PHASE 
MAGNETIC AMPLIFIERS® 


A. L. Pisarey 


(Moscow) 


The paper deals with two three-phase magnetic amplifiers circuits, which differ 
by the method of connecting the phase control windings. The static operation of the 
circuits is described. Analytical expressions for the control characteristic of the ampli- 
fiers of both types are obtained, On the basis of this analysis the conclusion is arrived 
at, that it is inadvisable to use a circuit with series connected phase control windings. 


Magnetic amplifiers used for controlling three-phase loads can be connected in various ways. In some of 
the circuits (for instance when the magnetic amplifiers are connected in series with delta connected phase loads) 
their operation is similar to single-phase working and does not require a special analysis. 


The present paper deals with magnetic amplifier circuits 
without a positive feedback whose operation is very different 
from single-phase working. These circuits with a star conmnec~ 
tion of the magnetic amplifiers are often used, but to date the 
operation of these circuits has not been properly investigated. 


Each of the above circuits consits of six reactance coils 
connected a pair per phase in series with the resistive load. The 
load windings of each pair of coils are connected in series . 
Each pair of coils (phase) can have either a common control 
windings or separate control windings connected in series op- 
posing. In the first of the two circuits (Fig. 1), the phase con- 
trol windings are connected in parallel and fed from a supply 
with a small internal resistance. If only a source with a high 
internal resistance is available, it is sufficient to shunt the in- 


Fig. 1. A magnetic amplifier circuit 


with its phase 1 windings con- put with a capacitor, 

nected in parallel. In the second circuit (Fig. 2), the phase control windings 
are connected in series and are fed from a source with a large 
internal resistance. 


The first circuit has free evenharmonicsof the current in the control circuit, in the second circuit the 
harmonics are suppressed. 


The above two operating conditions exist in their pure form when the afore-mentioned stipulation regard- 
ing the internal resistance of the source is observed. With a finite resistance of the supply source, the operatihg 


*The paper was read at the seminar on magnetic amplifiers and contactless magnetic elements at the Institute 
of Automation and Telemechanics of the Academy of Sciences, USSRon March 26, 1958. 
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conditions of the circuit in Fig. 1 approximates that with 
free even harmonics of the current in the control circuit 
and the operating conditions of the circuit in Fig. 2, ap= 
proximates that with suppressed harmonics. 


The following assumptions were made in the analysis; 


1. The reactance coils have an ideal magnetization 
characteristic (Fig. 3). 


2. The resistance ¢ the reactance coils is negligible. 


3. The supply (phase) voltage conforms to equa~ 
tion: 


4 
Uy, = ww: 10°8. (1) 


Fig. 2. A magnetic amplifier circuit with its Reasons for selecting this value for Um are given 
control windings connected in series. below. 


Let us now examine the two afcre- mentioned 
magnetic amplifier circuits. 


Il. A Three-Phase Circuit with the Phase Control Windings Connected in Parallel 


The circuit is shown in Fig. 1. Let us adopt the following notations: ug, Up, and Ug are the phase voltages 
(instantaneous values), Up) and U are the phase voltages (the amplitude and mean values), ugh, Ube and Ugg are 
the line voltages, ¢aq is the emf in the reactance coil a of phase a (and similarly for other phases), i,, 1, and ig 
are phase currents (instantaneous values) iy, is the current in one control phase winding, I}; and I are the mean 
values of the load and control currents, R is the phase load resistance, wy and wy are the number of turns in the 
load and control windings. 


Ug = U,, sin 8, (2) 
uy = Umsin(@—=*), 


Uc = Umsin( —*), 


Us = V3U,, sin (o—*), 
Ure = V3U sin 


Uca = V 3U m sin 


For circuits with parallel phase control windings the following propositions, known from the theory of single- 
phase magnetic amplifiers with free current harmonics in the control winding will hold. 


1, When one of the coils of a phase is saturated,the flux of the other does not change. 
2. During the change of the coil fluxes of the phase: 
d®, dD, 


dt (8) 


_8. In the coils of each phase the ampere~turn law is preserved (providing both coils are not saturated simul- 
taneously): 


| a 6 c 
| 
(5) 
+ lypwy = 0. (9) 
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ia + iy + ie = 0. tio 


In above circuits the following magnetic conditions can exist 
in the coils: 


| | 1. The coil of one of the phases is saturated. Its paired coil 
% has an unchanged flux. The fluxes of the remaining two phases 


change. 
hy 2. In each of the two phases one coil is saturated. The fluxes 


of their paired coils remain unchanged. The fluxes of the third phase 
coils change. 


3. In each of the three phases one coil {s saturated. In the re~- 
maining phases the fluxes are unchanged. 


4. All the coils of the circuit are saturated. 
Combinations of the above magnetic states of coils produce 


three operating conditions for a magnetic amplifier and the boundaries of these conditions are determined by the 
angle of saturation ( 6s): 


Fig. 3. Approximate magnetization 
characteristic of a reactance coil. 


| K 


A 


Fig. 4. Voltage and current diagrams for a magnetic amplifier with a Fig. 1cir- 
cuit. 


also holds. 
“hy 
=# 
(12) 
(33) 
| 
c 
(9) 


Angle 9, cannot exceed 51/6 since otherwise condition (10) cannot be fulfilled. 
Let us examine the operating conditions of the magnetic amplifier (Fig. 4). 


The First Mode. In this mode, coil conditions one and two alternate. 


In the interval 0 <9 < 1/6 (Fig. 4, A) coil aB is saturated, and the flux of coil aa has a changed value. 
Coils ca and cB are in a similar condition. 


The flux of coil b8 is increasing and that of coil ba is decreasing. The following equations hold: 


= — AM, = 


1.5 up, d6, (15) 


(16) 
(17) 


2R 


ip = 0. 
The validity of equation (15) follows from Fig. 5. It is obvious that in this interval the potential of the 


zero point in Fig. 1 circuit lies in point O’ (Fig. 5). Coils ba and bB have voltages which are represented by 
vector O°B applied to them. 

In the interval ¢/6 < 9 < 4g, the flux in coils ac and b§ increases and in coils aB and ba decreases. 
Coil ca is saturated and the flux of coil c8 is unchanged. The emf ecq= 0, and ecg = 0. 

Coil ca remains saturated despite the load current i, having stopped flowing through it. This is explained 


by the fact that a decrease in the flux éca will only be possible when voltage u, becomes negative. In this 
interval the following equations hold: 


(18) 
Ure + op + €pa = 0, (19) 


sb 
108 
= — = = Ure 


ig=h =i, = 0. 


In the interval 9,43 <9 < 2/2 coils b8 and ca are saturated, the flux in coils ba andc8 is unchanged, 
the flux in coil aa increases and that of coil aB decreases. The following’ equations hold: 


(24) 


Ure — if 4+- O, 


(25) 


= — AD = 0, (22) 
(23) 
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25) 


26) 


In the interval w/2 < © < Osaq the coil b6 is saturated, the flux in coil ba is unchanged, the flux in coils 
‘a6 and ca rises and the flux of coils a6 and ca falls. The following equations hold; 


(27) 
(28) 


(29) 


(30) 


= — = 
i, = i, = i, = 0. 
Then coil aa saturates. 
The processes and the equations., which represent them,repeat themselves with different subscripts. 


B The end of the first mode will correspond to 9, = #/2. Moreover 
the coils will not be in their first magnetic condition. 


It will be seen that in the first mode there are no instances when 
both coils of one phase are saturated at the same time. Hence the law 
of ampere-turns is preserved (9). 


The control current in each phase follows the equation: 
iyg= + 
Fig. 5. Vector diagram for the 
voltage in the interval 


<O< 


where the plus sign corresponds to the positive half pericd of the load 
current and the minus sign to the negative one. The control current 
control currents: 


ty = bya + typ + bye. (34) 


Let us examine the mode of operation of a magnetic amplifier without superposed magnetization. Let us 
show that if the supply voltage is chosen according to (1) the flux of each coil will vary from @; to —@,, more~ 
over, coils of different phases will not be saturated simultaneously. 


This mode is the limiting case of the first mode of operation and corresponds to saturation angle values of 
6, = 54/6. Moreover at any instant iz; = 0. Equations of the first mode, which have a load current equal to zero 
simultaneously in all the phases,(i.e., the families of equations (18)-(23) and (27)-(32))will satisfy this mode. 
Moreover, 9shg = 4/2 and will hold. 


according to (29). 


Substituting (77) in (20) and (5) in (29) and integrating we obtain: 


®saa 
. 108 
= — = — \ Und, 
© A® 
= — tow \ = 
) LE (31) 
) (32) 
3) 
2) 
0) 
1) 
2) 
= at + <6 <> 
Ao’ V3U,, 108 
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on 


° V3U,,-108 
= 


The total change of flux in the interval 1/6 <@ < 51/6 is determined by expression: 


V3u,,- 108 


Choosing Um according to (1), we obtain from (37): 


=> 29,. (38) 
This value of Um is limiting where there is no simultaneous saturation of coils in two phases,and hence 
in the absence of a saturation current conditions ij; = 0 always holds. 


From the above the conclusion can be drawn, that, if the coils are connected according to Fig. 1 and with 
a given supply voltage Um, the design voltage for a coil will be Up, /V3 and not U,,,/2 as it would appear at 
first sight. 

It can be shown that this conclusion also holds for circuit of Fig. 2. 


Second Mode. The second mode of operation of the coils (Fig. 4, B) starts at 6,= 4/2. It differs from the 
first by the existence of only the second magnetic condition of the coils. The duration of the load current flow 
in this mode is constant and equals to 21/3. 


In the interval 0<@ < &),g the magnetic condition is the same as for the first mode at 0< 0 < 4/6. The 
corresponding equations hold with the exception of equation (15) in which the upper integration limit will be 
Osbp instead of x/6. At & = 6.48 coil bB is saturated, and at this instant one coil in each of the phases also 
_ becomes saturated. The circuit becomes noninductive if the assumed conditions hold. At this instant, however, 
the current in phase a is in antiphase with respect to its phase voltage and therefore it falls instantly to zero. 


Any change in the phase of the current i, is prevented by the changing fluxes @,, and ag. 


In the interval 6,p8 < 9 < Ssacx the magnetic condition of the coils is similar two that of the first mode at 
Osbp < @ <#/2, Equations (24) and (26) hold. In equation (25) the upper integrating limit is saa instead of 


After this the processes and the equations which represent them repeat themselves, with different subscripts. 
The second mode ends at 9s= 4/3. A further decrease of 9, produces the third mode in addition to the second. 
The control current in the second mode similarly to the first mode is determined by equations (33) and (34). 


Third Mode. The third mode of operation of the coils (Fig. 4, C) starts at 0,= w/3. It differs from the 
second mode by the alternation of the second and third magnetic condition of coils. The duration of the load 
current flow in the third mode is variable and always larger than 21/3. 


Let us assume that at © = 0 the coilsb8 and ca are saturated (Fig. 4, C) the fluxes of coils ba and cB are 
unchanged, the flux of coil ac is increasing and that of a8 is decreasing. In the interval 0 < @< 9,,,, the fol- 
lowing equations hold: 


Ud. 

ic = —ih= (40) 

= — = 1.5uad0. 

0 


— 
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In the interval @saq < © < 4/3 coils ac, bB and ca are saturated. The flux of coils a8 , be and cB re- 
mains unchanged. The following equations hold; 


Uca — igR + = 0, (42) 
Ute — i + i,R = 0. (43) 


By solving them in conjunction with (10) we obtain: 


(44) 
(4) 
(46) 


In subse quent intervals the processes and the equations, which represent them repeat themselves with dif- 
ferent subscripts. 


The third mode ends at 0, = 0. Moreover, 
a 4 4 the load current is sinusoidal and in phase with its 
phase voltage. 
we Le. In the third mode as in the first two modes 
NE the control current is determined by equations (33) 


and (34). 


Fig. 6. Oscillograms of voltages and currents for a 


Fig. 6, a, b, c shows oscillograms of phase 
magnetic amplifier with a Fig. 1 circuit. voltages, spoiding phase a ail , 


currents for the three modes of operation discussed 
above. It will be seen from Fig. 6, saat, the 4. 


Let us introduce the notation 


(48) 


In (48) L, is taken to be the mean value of the control current of one phase. Hence the value of uy from 0 
to 1 will satisfy according to (33) the equation: 


Uy = Wy. (49) 


The design control characteristic of a magnetic amplifier uy = f (uy) is shown in Fig. 1. The dots on the 
same graph show the same characteristic obtained experimentally. 


The gain of the amplifier based on ampere~-turns for a working interval of the characteristic (0 < uy < 1) 
is equal to: 


kaw = (50) 


The relationship kaw = f (uy) is given in Fig. 7, 


The current gain is: 
dl, 


a, (51) 


Uy 
b= 
u. 
8) i. = Ez 
ith 
the 
4 
er, 
wy 
cond. 
e 
b 
40) 
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In conjunction with (50) we obtain: 
ky = (52) 
The power gain,according to [1},can be represented in the form: 
ky = kaw kj. (53) 
It is known that a coefficient of the form kf is represented by the expression: 


, 
Jeff (1/*f Geo) 


The design relationship between kf and wy is given for this case in Fig, 7. 


Il. A Three-Phase Circuit with the Phase Control Winding Connected in Series. 


This circuit is shown in Fig. 2. The notations for voltages,amfs and currents is the same as for the circuit 
in Fig. 1 


. The alternating magnetic conditions of coils produce 
Henke Uy in this instance four modes of operation of the magnetic 
amplifier. 


First Mode. In this instance, there can be the following 
magnetic conditions of coils in one phase: a) both coils are 
saturated; b) one coild is saturated and the other is not. 


In case b), equation (33) holds. Hence it can be as- 
sumed that, when equation (33) does not hold,both coils are 
saturated. 


Since during the operation of the circuit in each phase 


My one of the coils is saturated and also because the coils, as it 
Fig. 7. Control characteristic of a mag- will be shown below, are acted upon by line voltages with an 
netic amplifier with a Fig. 1 circuit. amplitude of ¥3Um, there is a complete reversal of mag- 


netization of the unsaturated coil. Moreover, its flux changes 
from + @s to —@,, remain at that level during a time, corresponding to angle y and then grow from ~¢s to 
+. 


Let us denote the phase differences between phase voltage and the beginning of the current flow in that 
phase by a. The graph of voltages and currents for the first mode is shown in Fig. 8, A. 


Let us write the equations of electrical balance,interval by interval,in order to determine the emf of one of 
the coils (for instance, that of eg, ) during all the portions of the interval when the emf does not equal zero. 


1, a< © <a+y (Fig. 8, A). Coils ba, bB , act and ca are saturated. In coils a8 and cB the flux is 
changing. The following equations hold: 
ip = — (ia + ic) = — 21 y, (55) 
Uar — inh + igR + Cap = 0. 


(56) 


B Cc 
Fig. 8. Diagram of voltages and currents for various phases of a magnetic am- 
plifier with a Fig. 2 circuit. 

2. a+ y< @<a+n/3. Coils ca, cB ax and b6 are saturated. In coils aB , and ba the flux is changing. 

The following equations hold: 


(57) 


(58) 


3. a+"/3<0<a+n/3+y. Coils aa, a6 , b8 and cB are saturated. The flux in coils ba and ca is 
changing. The following equations hold: 


Uca — Cag — igh = 0. 


eap = 0. (60) 


4. y< < a+ 2n/3. Coils ba, , ax andcB are saturated. The flux in coils a6 and ca is 
changing. The following equations hold. 


a= (61) 
Uar + +- Cag = 0. (62) 


5. a+2n/3< 6<a+2n/3+y Coils ca, cB, ax and ba are saturated, in coils a8 and b8 the flux is 
changing. The following equations hold: 


ie = — (63) 
Uca — — + = 0. (64) 


Solving the equations (56), (58), (60) and (64) with respect toeag and changing the currents in the load 
circuit by control currents, we shall obtain for the above five intervals the following relations. 


= — — (85) 


(66) 


le \" 
7 
| 
Cap = Uca — TyR = 
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Cag = 0, 
Cap = — Ug — Ly —t, 
H 


Cap = tga — (69) 


In the interval a<@<a+ 4/3, the flux in coil a8 changes from + @, to —&,. In the interval a+ 2/34 
+ y< <a+2n/3+y it changes from to + 


The increment of flux in the core is related to the emf at the coil terminals by the equation: 


(70) 


From ('70) on the basis of (65), (66), (68) and (69) and in conjunction with (1), (5) and (7), it is possible to 
write: 


E V 3U » sin (0 —3lyR 4 


v3 
rx w. 


a+ 


After integration between the limits expressions (71) and (72) assume the form: 


V3 U moos (a + (4+ (74) 


Let us determine the boundaries of the first mode. The mode will start with a very small value of the 
control current ly 0. 


Fig. 9. Oscillogram of voltages andcurrents ina magnetic amplifier with 
a Fig. 2 circuit. 


(68) 
= — (edd. 
a+ = 
| | (0—F)— 1 
ty 
a+ +y 
[v3u., sin (o— “2 a0 —¥3 
+ \ m sin ( y (72) 
V5Umsin +1) — UyR =P Um 
7 
ly 
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L) 


2) 


3) 


4) 


Assuming that in equations ('73) and (74) ly = 0, we find that a = 90° and y = 60° at that point. 
Moreover the load current curve has no discontinuity. 


Angle y decreases with a rising control current, since the voltage applied to each coil is decreased by 
quantity TyR( Wy / Wy) and the time of negative saturation also decreases,thus reducing angle y. 


Let us find the other extreme point of the first mode, a point determined by the conditions y = 0. We find 
from (73) and (74) that for this point a = 66° and Uy = 1.06. The oscillogram for this mode is given in Fig. 9, a. 


In conjunction with (48) equations (73) and (74) can be reduced to 


sin (« + 7) — 0.736 wy(y + %) = 0.5, (75) 


cos (a + — 0.736 + 7) =—05. (76) 


In equations (75) and (76) Uy is the independent variable and a and y are functions. By taking various 
values of Uy in the interval of 0< uy < 1.06, 1.e.,in the limits of the first mode it is possible to obtain correspond 


ing values of a and y. The solution of these equations, however, is very difficult. In practice it is easier to take 
y as the independent variable and thus determine the values of a and uy. 


kenty D On the basis of the above and Fig. 8, A it is pos- 
20y 10 al Bs sible to obtain a characteristic control equation for the 
6 Fate : amplifier in the first mode of operation in the form: 


Uy = (5 +2). (77) 


04+ |. = if Fig. 10 shows the calculated control characteristic 
P | ie \ of an ideal amplifier. On the same diagram the experi- 
wo iby, mentally obtained characteristic is given in dots. Length 


Fig. 10. Control characteristic of a magnetic OA corresponds to the first mode. 


amplifier with a Fig. 2 circuit. Second Mode. With ty = 1.06, the second mode of 
operation of the magnetic amplifier comes into force 

(Fig. 8, B),a mode characterized by condition y = 0. Physically this means that in the coil whose flux is chang- 

ing, there is not a complete reversal of magnetization (from + @, to ~@s). The load current has the shape of 

rectangles with ordinates of ly(wy/ w};) and abscissas of 24/3. The control characteristic of the magnetic ampli- 

fier is represented by a straight line inclined to the X-axis at an angle of =tan-'2/3= 33°40". The beginning 

of this mode will occur when Uy = 1.06 (point A in Fig. 10). 


The second mode ends when coils of two phases become saturated simultaneously. This instant will cor- 
respond to the beginning of the third mode and will be determined from equations of the third mode. For de~- 
termining u,, in the second mode the following equation is used; 


= (78) 


The second mode corresponds to length AB in Fig. 10. The current oscillogram for the second mode is given 
in Fig. 9, b. 


Third Mode. The third mode is characterized by the coils of two phases becoming saturated over a part of 
the period B (Fig. 8, C), one of the third phase coils being also saturated and the other having a changing flux. 
In order to determine the relation between quantities a, 6 and uy let us find,similarly to the first mode, the ex- 
pression for the emf of one of the coils (for instance €qg) from the equation of the electrical balance of the cir- 
cuit in intervals. 


The sum of the integrals of these interval emfs will provide a measure of the total flux in the coil, which 
is equal to zero, since the flux starts changing from a saturated condition,and again returns to a saturated condition. 
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The flux a6 changes ir the interval a+8 <0 < a+ 20/3 (Fig. 8, C). The equations of the electrical 
balance in intervals will be the following: 


1) coils ca, cB, axand are saturated, 


Uca — Cag — ight = 0. : (79) 
2) a+ coils ca, cB, ba, bBandaw are saturated, 
Uca — Cap — taht + i-R = 0. (80) 


3) a+ 5 +B<0<a+ coils ba, bB, cBanda« are saturated, 
Ug» + igh + €ag = 9. 
The emfs for the corresponding intervals will be : 


Cap = Uca — ial, 


It is not difficult to find the current in the second interval: 


2 (86) 


From (85), (86), (5) and ('7) and equating to zero the sum of the integrals of the corresponding emfs in in- 
tervals, we shall obtain after transformations; 


VE [Sain $) e494 


The second equation for determining « and B by given u,,we shall obtain from the condition 


w. 
= (88) 
From the equation of the electrical balance in the interval a < 9 < a+ 8 (Fig. 8, C) we have 


Uca — ight = 0, 


In conjunction with (10) we obtain: 
Wy 


Uca I, 


By substituting (90) in (88) and considering (7) and (48) we have: 
uy = V3 sin (a+p—). 


(82) 
Cap = — Ue — if. (84) 
In the interval under consideration 

ig= TI, (85) 

(89) 
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5) 


(89) 


(90) 


(91) 


Taking for the beginning of the third mode 8 = 0, we find from (87) and (91) c= 55°20" and uy™ 1.18. 
These values of a and uy also correspond to the end of the second mode. 
For the end of the third mode assuming 8 = 4/3 we find a = 3°24" and uy = 1.49. 


By assuming values for uy it is possible to obtain from (87) and (91) a and 6 for the intermediate +4" 
and then find ujy. The calculated control characteristic for the third mode is given in Fig. 10 (length BC). 
oscillogram of the current in the third mode is given in Fig. 9, c .. 


Fourth Mode. When uy grows from 1.49 to 1.57 a new mode of operation of the magnetic amplifier is pro- 
duced ,. which differs from the third mode by all the six coils. being saturated in certain intervals. In these inter- 
vals the load current is in phase with the corresponding phase voltages. With the rise in uy the duration of these 
intervals increases and at uy = 1.57, all the coils remain saturated during a period. Thus, the load current becomes 
sinusoidal and in phase oe the corresponding phase voltage. Length CD of the amplifier control characteristic 


corresponds to the fourth mode of operation. The oscillogram for the beginning of the fourth mode is given in 
Fig. 9, d. 


The gain,according to the ampere~turns character- 
istic, varies with the parts. 


The relation between k, y and Uy is shown in Fig. 10. The current gain changes accordingly. The power 
gain similarly to the previous case is determined from formula (53). 


The relation of the current form factor (54) to uy is shown in Fig. 10. 
On the basis of the above, the following conclusions can be drawn: 
1. The value of the supply voltage for both of the circuis considered must be chosen according to equation(1). 


2. The control characteristic of a circuit with its phase control windings connected in parallel has a slope 
of 45° constant over the whole operating range. 


slope. 


For complete saturation of the magnetic amplifier in this circuit 1.57 times greater ampere~turns are required 
than for the first circuit. 


4. Tests have shown that in both circuits the control characteristic with a complex load is almost the same 
as with a similar resistive load. . 


5. When a three-phase amplifier circuit is made up of three single-phase amplifiers, their control windings 
provide savings in copper, but have a lower gain as comparedwith a groups of single -phase amplifiers, connected 
according to the circuit in Fig. 1 with a nonlinear control characteristic. 
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DETERMINATION OF THE BASIC PARAMETERS OF 
ELECTROMAGNETIC RELAYS 


M. I. Vitenberg 


(Leningrad) 


Most important relations between nominal work, consumed power, overheating 
temperature, dimensions and weights of electromagnetic relays are considered. 


Experimental curves are appended and empirical formulas deduced for determin- 
ing the operating conditions of a relay in terms of the cross-section of its core, power 
consumed, overheating temperature and its weight. 


The design of an electromagnetic relay is a difficult and complex question which requires the solution of 
very complicated and diverse problems. One of the most important problems in designing a relay is the correct 
determination of the main dimensions of its magnetic system in terms of loading and the power consumed. 


The general design technique of the basic magnetic system parameters of relays has been as yet insuf 
ficiently developed and does not provide the possibility of determining with an accuracy sufficient for practical 
purposes, the main dimensions, consumed power,and weight of the relay in relation to its nominal work. 


Nominal Work of an Electromagnetic Relay 


The force of attraction of an electromagnetic relay armature can be expressed by the following formula [1]: 


(1) 


where Gy is the full permeance of the relay, Gg is the total permeance of the effective air gap related to the 
armature movement, AW ~ the effective ampere~turns and § ™ the distance between the armature and core (ef- 
fective gap) of the relay. 


Force F in formula (1) is expressed in newtons. 
If the stray fluxes are neglected, the total permeance of the relay becomes: 


i 


where G,, is permeance of parts of the magnetic circuit, independent of the armature movement (core, body, 
base, armature, butts and nonmagnetic distance pieces). 


The nominal work of the relay is at a maximum, when the two permeances Gg and Gp are approximately 
equal (Gg~wG,,) [1]. In this case: 


AW*G2 dG. 
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(2) 


y, 


ately 


Hence in case of an optimum relationship between the reluctance of the variousrelay components the for- 
mula for the force of attraction will assume the form: 


R Aw? 4G, : 
(3) 


If the length of the effective gap & is small in comparison with the core cross section area S, it is possible 
to assume that: 


where pig™ - 10°" ts the reluctance of vacuum. 


In this case we obtain for the optimum force of attraction of the relay expression: 


AWG? 


The maximum value of the nominal work of the relay will, obviously, be: 


(5) 


In the last formula C is the meanresistance ofa turn of the winding, and P is tha gouge coummed by Se 
relay winding. If the dimensions are expressed in cm, the power in watts, and jig™ 4x- 10~*, the optimunr value 
of the nominal work will be: 


A nr AWG, 


It follows from this formula that if there is no saturation of steel, the flux is evenly distributed in the ef- 
fective air gap and the effect of stray fluxes is neglected, the relation between the optimum nominal work in a 
relay of the same type and the power applied is linear. 


The thermal conductivity of the surrounding medium to the heat dissipated by the winding is: 
G = 
The electrical conductance of the winding port is: 
G =+. 


Substituting in equation (5 a) for P and C their values obtained from the latter expressions we find for the 
nominal work of the relays a relationship similar to the one obtained by B. S. Sot-skov [2] 


Aop = 1,6 Ls. (5 b) 


«= 
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Three different types of relays with a swiveling armature (of the drop type) and four relay models with cores 
of different diameters in the range of 2-18 mm were used for testing. The main data of all the relays used for test~ 
ing are given in the table attached. 


Type 


Armature 

ppovements 
, mm 

Distan 
eces 
m 


Core 
diameter, 
mm 
me 
essss A 
if 


Relay model No, 1 
Relay type RKN 


Omoooso 


Figure1 shows the experimentally obtained relations between the nominal work and the power consumed for 
seven different types of relays with a swiveling armature. These curves were taken with a normal armature moye-. 
ment and distance pieces for each relay. 


a G*cm 


700 


a aa7i 23 57H BHP w 


Fig. 1. Relation between the nominal work of the relay and power 
co.sumed. 

1) Relay model No, 1(d=18 mm); 2) relay type RKN; 3) relay type 
RKM; 4) relay type RSM; 5) relay model No. 2(d= 3 mm); 6) re- 
lay model No. 3 (d= 2.4 mm); 7) relay model No. 4(d = 2 mm). 


Two versions of the relay type RKN and the model No. 1 were used, the normal (with pole pieces) and the 
special (without pole pieces). The curves for the first two types of relays in the table without their pole pieces 
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are shown in Fig. 1 by dot-dash lines. The ratio of the core length to its diameter in the first three types of the 
table is about 7.8 and in the four last types, about 3.2. Therefore, the value of coefficient C (the characteristic 
resistance of one turn of the value of winding) is smaller in the first three than in the last four relays. 


Hence in order to eliminate the effect of the ratio 2/d on the power consumed, the curves showing the re- — 


lation between the nominal work and the power consumed for the last four types of relays were reduced tothe values 
of 1/d= 7.8 and for the first three types were reduced to 1/d= 3.2. 


The corrected curves of the nominal work are shown in Fig. 1 by a dotted line and marked by corresponding 
numbers with a dash. From these curves it follows that the relation between the nominal work of the relays and the 
power consumed {s not linear, especially so for small relays. (In logarithmic scale of Fig. 1, a linear relation {s re- 
presented by a straight line drawn at an angle of 45° to the X-axis). 


Relays type RKN, RKM and model No. 1 have curves which at their beginning are close to a linear relation- 
ship, but they bend in their higher portion due to the saturation of the core. 


Relation of the Operating Conditions to the Cross Section (Diameter) of the Relay 
Core. 


The values of quantities Gg and C, which appear in formulas (4) and (5), depend on the construction and 
dimensions of the relays 


inp (Dy + h) 
= = and C= 


In these formulas d, is the diameter of the pole piece, Dg is the internal diameter of the winding,J, is 


the length of the winding, h is the height of the winding, p {s the resistivity of the wire, K, is the stacking factor 
of the winding 


Quantities Ds, h and?) can be considered proportional to the diameterd of the windings 
= h = adandl, = a,d. (8) 


For the sake of simplicity let us assume that the value of the effective gap and the diameter of the pole 
plece are also proportional to the core diameter: 


= andd, = ad. 


Substituting in the expression for Gg and C the quantities 8, Dp, h and 2\,by their values we obtain: 


nad mp (1 + a) 
= OgdandC =: (1) 


Hence the optimum nominal work of the relay is 


(8b) 
where S is the cross-sectional area of the core. 


Thus, with the assumption made above and constant ampere~turns, the nominal work of the relay is pro~ 
portional to the diameter of the core, and with a constant supplied power, it is proportional to the square of the 
diameter, i.e., the cross-sectional area of the core. 
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Figure 2 shows the relations between the nominal work and the core cross~sectional area at various values of 
the supplied power, plotted on the basis of the experimental data obtained in Fig. 1. The curves have been plot- 
ted separately for the two values of the ratio of the core length to its diameter i.e., 7.8 and 3.2. The dotted line 
in Fig. 2 correspond to the dotted portions of the curves in Fig. 1. 
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Fig. 2 


From the curves in Fig. 2 it is clear that for each value of J/d, the relation of the nominal work and the cross- 
sectional area of the core is practically linear only at a power of about 0.2 w,and cross sections larger than 


0.07 cm”. (On the logarithmic scale of Fig. 2,linear relationst.'p corresponds to a straight line at angle a = 26.5° 
to the X-axis). 


% 
we 
> 


Scm? 


At powers smaller or larger than 0.2 w the relation of the nominalwork of the relay to the cross section 
of the core (at S = 0.07 cm?) is not linear. At cross sections smaller than 0.07 cm’, the nominal work of the re- 
lay drops‘rapidly. With a rising power supplied to the relay, its nominal work increases. The maximum nominal 
work of the relay is determined by the maximum permissible temperature of its winding. Hence in the formula 


for nominal work it is advisable to substitute power by the corresponding excess temperature of the winding above 
the ambient air temperature. 


It is shown in [3] that the relation of the excess temperature of the winding above the ambient air, to the 
power supplied and the weight of the relay,can be expressed by the formula: 


(9) 


P= 


For the relay of the drop type weighing between 1 and 2,000 g, the value of coefficient 9), is equal approxi 
mately to 176. By substituting in formula (8 b) quantity P by its value obtained from (9), we have: 


Aop = d?, (10) 


| 
whence 
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Relation of the Relay Weight to Diameter (cross section) of the Core 
The total relay weight consists of the weight of the core, wire, system of contacts and other components: 


= Ou + 


The weight of the contacts, base, cover, and other constructional details can be accounted for by a coef- 
ficient Ky. Thus, the total weight of the relay can be expressed by the formula; 
(11) 


Q =K, (Qo+ Qu): 


The weight of the core can be found from the expression: 
Qc = LS 


where L is the mean length of the core, S is its cross-sectional area, y, is the specific gravity of steel and K, is 
the coefficient accounting for the irregularities of the core cross section along its length and for protruding por- 


tions of the armature and body of the relay. 
The mean length of the magnetic circuit is equal approximately to: 


where J is the length of the core, J, is the length of the effective part of the armature. 
For the sake of simplicity let us assume that the length of the effective portion of the armature is proportional 


to the diameter of the core: 1, = agd. Then: 


L = 2agd + 2agd = 
cCross- 
Qo = KydS = 


The weight of copper in the sdnding to known formula: 

tna Qu = (Do + h) K "hw 

“7 where yp, is the specific weight of copper which is equal to 8.9. 

Substituting 2,, h and Dg'by their values in terms of the core diameter we obtain: 
Que = (1 + a,) K ty = 


) Hence the total weight of the relay can be expressed by the formula: 
Q = Ky + Kya?) = (14) 


Thus, the relay weight with the assumptions made is proportional to the cube of the core diameter. 
With a ratio of the core length to its diameter of 7.8, the yalue of coefficient Ky for telephone relays lies 
approximately in the limits of 325-400. Hence with a ratio J/d= 7.8 the weight of a relay can be approxi~- 


mately determined from the following formula; 


= (325 — 400) fg}, (14a) 


LO) where d is the diameter of the core in cm. 


ues of 
line 
| | 
615 


With a ratio of the core length to its diameter of 3.2, the value of coefficient K, for small relays of the 
drop type is approximately in the limits of 180-200. similarly with a ratio Z/d = 3.2 the tentative estimation of 
the relay weight will be: 

Qs,2 (180 — 220) a’. 


(14b) 
Relations of the Nominal Work to the Weight of the Relay 


Substituting ia formula (10) quantities Q and d by their values obtained from (14), we have: 


A 


079 
op 


on d? => = (15a) 


= (/ z) = = 


(15 b) 
Thus, with the assumptions made, the nominal work of the relay is practically proportional to its weight and 
the excess of the temperature of its windings above that of the ambient air. 
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Figure 3 shows the experimentally obtained curves for the relation between the nominal work and the cross 
section of the core at yarious values by which the winding temperature exceeded that of the surrounding air for 
two values of ratio 1/d. The dotted portions of the curves correspond to the dotted portions in Fig. 1. 

It follows from Fig. 3 that the curves of the relation between the nominal work of the relay and the core 
cross section at S > 0.07 cm? are represented on the logarithmic scale by practically straight lines, but these 
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‘lines for different excesses of temperature and different values of /d do not coincide and are not parallel to each 
other. At cross sections below 0.07 cm? the curves of the nominal work bend sharply downwards. 


Figure 4 shows the experimentally obtained curves of the relation between the nominal work and the weight 
of the relay, at different values of excess temperature of the winding over that of the surrounding air, for two values 
of ratio Z/d. Dotted portions of these curves correspond to the weighted (dotted) portions of the curves in Fig. 1. 


a,G-* cm 


30 50 70 100 200 300 500 7001000 2000 20006, 
Fig. 4 


It follows from Fig. 4 that,at constant values of excess temperature of the winding,the relation between the 
nominal work of the relay and its weight within the limits of 5-3,000 g is expressed on the logarithmic scale 
by practically straight lines. The variation of the ratio of the core length to its diameter in the limits of ap- 
proximately 3.2-7.8 hardly affects this relationship. The deviation of the curves in Fig. 4 from a straight line at 
small weights (2 to 5 g)is due to the relatively large movement of the armature (0.4 mm) in relay mc¢wels No.5, 
6, and 7. It should be also noted that the curves of linear relations of the nominal work of the relay to its weight 
for various excess temperatures are not parallel to each other and are inclined to the X-axis at different angles. 


Thus, the relation of the nominal work of the relay to its weight in the limits of 5 to 3,000 g,at a con- 


stant excess temperature of the winding above that of the surrounding air,can be approximated by the formula 
similar to (15 b): 


A = KySQ? = Gen], 


where Ag is the coefficient, whose value is equal to the nominal work of the relay weighing 1 g at the given ex- 
cess of temperature when the relation of A to Q is also linear in the limits of 5-1 g,and b is the tangent of the 
angle of slope of the nominal work curve to the X~axis. 


According to the assumption made in formula (15 b), the values of coefficients Ky and b do not depend 
on the excess temperature(b= 1.17). In practice, however, coefficients Ky, and b depend on the excess tem~ 


perature of the winding in fairly wide limits. With an excess temperature of the winding = 70°C, it is possible 


} 
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to express the relation of the nominal work of the relay to its weight in the limits of 5 to 3,000 g by the formula: 
Ary = 0.55 95 = 0,41 Q, 16) 


where Q ts the weight of the relay, 


With a temperature excess of the winding of 50, 30 and 10°C, the relation of the nominal work of the relay 
to its weight within the same limits can be expressed repectively by the following formulas: 


Ago = 0.4199 0.36 Q, (16 a) 
Ay = 0.28 0,30 Q, (16) 
and 
Ai = 0.1291 0,199. (16 c) 


Relation of the Power Consumed by the Relay to its Nominal Work and its Weight 


When designing electromagnetic relays, it is important to know the relation between its sensitivity (the 
power required to work the relay at a given load) and weight. For determining the relation between the power 
consumed by the relay and its weight let us use the formula (8b): 


A 


Substituting in the last expression d by its value obtained from (14), we find the formula for the power 
consumed by the relay in terms of the nominal work done and its weight: 


From the last expression we obtain for the nominal work of the relay the following formula: 


(17) 


(17a) 
A En 
Thus, the nominal work of the relay is proportional to the power consumed and the cube root of the square 
of its weight. 


Figure 5 shows the experimentally obtained curves of the relation between the power consumed by the relay 
and its weight at various values of the nominal work for ratios of the core length to its diameter of 7.8 and 3.2. 
Dotted lines correspond to weighted (dotted) curves of Fig. 1. 


It follows from Fig. 5 that the relation of the power consumed by the relay to its weight is not linear on a 
logarithmic scale, but follows a miore complicated law. The power consumed by the relay also depends on the 
ratio of the core length to its diameter. 


At small values of nominal work (up to 10 G- cm) for Jd = 7.8 and a relay weight between 100 and 3,000, 
the relation of the power consumed by the relay to its weight, is almost a straight line on the logarithmic scale. 
Hence in the limits mentioned above, this relation ¢an be approximated by the formula similar to (17): 


P = = (17b) 
From the curves of Fig. 5 we find that with nominal work of 2 G- cm, and a ratio of 1/d = 7.8, the relation 


between the power consumed and the weight of the relay in the limits of 100 to 3,000 g can be expressed by the 
formula; 
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P = 0.199080 1018 aon 
Vou “VQ 


where Q is the weight of the relay. 


Hence at small load (up to 5 G- cm) the power consumed during the operation of the relay is roughly in- 
versely proportional to the cube root of the relay weight (within the limits of 100 to 3,000 g). With a falling relay 
weight and rising nominal work the index of Q increases. Only when the weight of the relay is smaller than 50 g, 
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Fig. 5 


the ratio [/d = 7.8 and the nominal work is greater than 50 G- cm, does the power consumed by the relay become 
proportional to the cube root of the square of the weight. With 1/d= 7.8 and nominal work between 0.5 and 

5 G-cm, the relation of the power consumed to the weight of the relay in the limits of 100 to 3,000 g can be ex- 
pressed by the approximate formula: 


0.08A 


P ~ 0.095 AQ-°38 = — Ws 8b 


It is possible to obtain from the latter expression an approximate formula for the conditions of relay operation 
in the limits of 0.5 to 5 G.cm at I/d = 7,8 and the weight of the relay between 100 and 3,000 g: 


A=125P VQ. (19) 


Thus, the nominal work of a relay at small loads (up to 5 G-cm), a ratio ZAl = 7.8 and its weight between 


100 and 3,000 g is roughly proportional to the cube root of the relay weight and proportional to the power consumed 
by its winding. 
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BRIDGE RECTIFIERS WITH MAGNETICALLY COUPLED LOADS 


O. G. Malkina 


(Moscow) 


It is shown that bridge rectifiers with a negative inductive coupling between their 
loads, which consist of large resistances, can be used for measuring the R and X compo- 
nents of en impedance Z. With suitable selection of parameters it is possible to obtain 
an independent measurement of R and X . 


As an example, current changes are examined in identical circuits and in circuits 
whose bridge ac diagonals are connected to impedance, which differ from each other by 
a small quantity jAx. 


Bridge rectifiers with magnetically coupled loads can serve as control circuits for magnetic amplifiers in 
various automatic tracking and control devices, One of the windings of such an amplifier receives rectified signals 
of the controlled variable,and the other, of the command variable. In order to obtain at the output of the mag- 
netic amplifier a current proportional to the difference of the signals, the control windings must produce reversed 
magnetization. In a simplified manner the control circuits are shown in Fig. 1, without the load and other windings, 

on the same core and fed from the same supply where: R is 


d resistance in the ac diagonal of the bridge rectifier, ly is 
4 R OP xy b the resistance in the dc diagonal, Rg is the resistance of the 
a, uur———8 P supply circuit, AZ is the measured impedance, L is the in- 
ath 4 Le ductance of the windings of the main and auxiliary circuits, 
| . ®, is the flux produced by the ampere-turns of the main 
c l circuit winding, and 4, is the flux produced by the ampere- 


a % turns of the auxiliary circuit winding. 
Pp, 
{*Aux. circuit” 


It will be shown below that such circuits with varying 
"Supply circute parameter relations can have a varying relative sensitivity 
mite (~) to changes in the resistances and reactances across the ac 
diagonals of the bridge and produce but a small distortion 
of the sinusoidal shape of the currents. These properties of 
. the circuits can be utilized in industrial automatic equip- 
ment for controlling production by means of measuring the R and X components of an impedence Z, 


The calculation of currents in these circuits with nonlinear elements and periodic emfs can be reduced to 
the examination of the stable state condition, which can be obtained in the form of several consecutive transient 
states [1]. The processes in these circuits are at first examined with ideal rectifiers, i.e. with rectifiers which 
have a forward resistance Rg = 0 and a reverse resistance Ry = 0, 


Above measuring circuits (Fig. 1) in the absence of the measured impedance AZ are identical and their 
impedances are (even without rectifiers) resistive, since the inductance of the one winding is eliminated by the 
mutual inductance of the other. 
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Fig. 2. Configuration of magnetically coupled 
circuits at different time intervals. 

a) During the rectification interval of a posi- 
tive (negative) half-wave of the supply voltage; 
b, c, and d) during the transition intervals. 


It will be shown below that the inclusion of the 
impedance AZ under test will produce changes in the 
equivalent resistances of the ac and dc diagonals of the 
bridge rectifiers of either circuit, 


This change occurs because, with the impedance 
under test included in the circuit, the coupled windings 
(even if the coupling is tight) will no longer produce 
a complete compensation of induction. 


A difference between inductions in the dc diago- 
nal of the bridge rectifier will produce free state cur- 
rents, These currents will arise during the transition 
from rectifying the positive half-wave of the supply 
voltage to rectifying the negative half-wave and vice 
versa, 


Owing to the small difference between the self 
and mutual inductances and the large resistance across 
the ac and dc diagonals of the bridge, the duration of the 
transient state is considerably smaller than a half-period 
of the supply voltage. In order to determine the measur- 
ing circuit currents, it is sufficient, therefore, to ex- 
amine the operation of the circuit during one half- 


period of the supply voltage. 


The free current, which appears during the change 
of the resistance from its forward to its reverse value 
and vice versa, produces distortions in the sinusoidal 
shape of the measuring circuit currents, The size of 
the distortions will depend on the phase difference 
between the current and voltage Uap in each circuit 
shown in Fig. 2, a. The circuit in Fig. 2, a will cor- 
respond to the one shown in Fig. 1 for the time in- 
terval, when one pair of rectifiers ad and be (or ac and 
db) is passing current and their resistance is equal to. 
zero, and the resistance of the other pair of rectifiers 
ac and bd (or ad and bc) is equal to infinity. 


The applied sinusoidal voltage will produce in the circuits of Fig. 2, a, ac currents, which can be found by 


solving the following simultaneous equations: 


1,Zo + 1,Z, lum 


I, oZo + 1,2, (1) 
I, + i,, 


Zo = Ry, Z,=R+ ry + jzry + AZ, 
Z,=R+ry (2) 
= Jiu = 


The currents in the measuring circuits will be: 


B, (a+ 
+ Zola + + (3) 
(more) 


p 
| 


Lon E, E, + 


+ Lala + Lids + (%) 


The voltage between points a and b of the circuit in Fig. 1 is across the diagonal of the bridge rectifiers, 
These voltages are determined by the expressions: 


Oy =) yry + — = 1,Z,st eq, 


1eq 


It will be seen from formulas (3) and (5) that the phase differences between the currents and voltages de~- 
pend on the circuit parameters, i.e., on impedance AZ and on the relation between impedances Xy» R and Ty: 


1, Case of AZ =Q, For AZ =0 the currents will be equal to each other, since Zteq will be equal to Zgeq and | 
the inductive components in the dc diagonals of both bridges will be produced by stray fluxes only. q 


With small stray fluxes and large resitances ty, the phase difference between the currents and voltages will q 
be so small that the change in the value of all the rectifier resistances will occur practically simultaneously, q 
Since the change occurs when the current is equal to zero, no free state currents are produced. 


The currents in the circuits will remain sinusoidal, since the self-inductance of the one winding is elimi- ) 
nated by the mutual inductance of the other. 


Distortions of the sinusoidal shape of the current will be due in this case only to the deviations of the recti- 7 
fiers from the ideal operation. 


If the volt-ampere characteristics of the rectifiers are approximated to short stretches of a straight line, the 
instantaneous voltage across the rectifier can be found from the expression: 


Us =U +ieRa, 


where Uat is the threshold forward voltage and Rq is the dynamic forward resistance of the rectifier. 


In order to obtain a practically sinusoidal current, the threshold voltage of the rectifier should not exceed 
3% of the source emf, 


2. Case whenAZ# 0, The current in the circuits at AZ 4 0 can be found by examining the consecutive q 
changes in the operating conditions of circuits with a changing resistance of the rectifiers. ty 


As anexample,let us examine the process taking place in the circuits at AZ = +jAx, where Ax << R, 


The phase differences between the voltages and currents across the bridge diagonal in the circuit of Fig. 2a, . 
will vary in this instance, They will depend on the relation of the circuit parameters, but the voltage across q 
the auxiliary circuit diagonal will always lead the voltage across the main circuit diagonal and the currents in | 
both the circuits, The phase difference between these currents and voltages leads toa time difference in the 4 
changing of the resistance values in the rectifiers from the forward to the reverse and vice versa, The resistance 
of the auxfliary circuit rectifier will change from the reversed condition to the forward one first, This change will 
take place when the voltage across the diagonal of this circuit bridge passes through zero, 


a 
q 
the 
od 
| 
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The measuring circuits with the new conditions of the rectifiers resistance are shown in Fig. 2,b. Currents 
and voltages after this change in the operating conditions of the circuit can be found from the following simulta- 
neous equation: 


igRe +4 (R-+ ry) + (ay + — 


= ry + ty — 


E g max sin (wt + 6,). 


Graphs of the currents and voltages variations in a circuit with ry = R, Ry = 2R, xy = 0.95xy, Ax = 0.1R 
are shown in Fig. 3. 


Sinusoidal currents and voltages of the circuit in Fig. 2, a are shown in a dotted line on the graph. 


It will be seen from the drawing that the second change of the operating condition of the circuit will occur 
at the instant voltage Up; passes through zero, when the main circuit bridge diagonal will be short circuited. 


The circuit of the new operating condition is shown in Fig, 2, c. Currents in these circuits can be found from 
the simultaneous equations: 
dis di, 


di dis 


+ At = &, lo = ig, @g= Ey max sin (wt + 


For the example under consideration, these currents are shown in Fig, 3 to the right of pointwt, = 0, 


The next change in the working conditions is ttre restoration of the rectifying action of one of the bridges. 
It occurs at the instant when the currents in the diagonals equalize and the current in the two rectifiers (with 
reversed resistance values) is equal to zero [2], 


It follows from Fig. 3 that the rectifying action of the auxiliary bridge circuit is restored first when i, = ~ fy. 


In any further consideration of currents it is necessary to remember that the equality of the currents occurred 
with a negative value of i, and that the rectffiers change the negative direction of the alternating current into a 
positive one with respect to i, and {,. 


The condition of the measuring circuits after the restoration of the rectifying action is shown for the example 
under consideration in Fig. 2, d. 


The instantaneous values of the currents in the circuits under the new conditions can be found from the equa- 


tions: 


+ + = eg, 


+ in (R + By — = 


— = 0, ig = + inser = Ey may Si (Wt + 0). 


The curves of the current variations in the circuits for the example under consideration are shown in Fig. 4. 


igRy + ™ 6g, 
i di, di, 0 
Wy + Ly = (1) 
+ — 0, | 
+ = &, (8) 
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Fig. 3. Graphs of the currents and voltages for the first and second 
change in the operating conditions of the rectifier circuits for the 
case under consideration. 


For the parameter relations under consideration the last change in the operating conditions of the rectifier 
circuit in the given cycle (half=period) will occur when current i, =i,. 


With large values of Ax the restoration of the rectifying condition of the bridges can occur at negative stable 
state values of voltages across the diagonals and the given change in the operating conditions of the circuit may 
not be the last. In the given half-period of the supply voltage there can appear repeated shorting of the bridge 
diagonals owing to additional points at which the diagonal voltages pass through zero. 


The condition after restoration of the rectifying action of both bridges will be similar to that shown in 
Fig. 2, a. 


The instantaneous values of circuit currents can be found from the following simultaneous equations; 


+ is (R + ry) + (ty + Az) — ay meg 
igho + ig (R — ty meg (10) 


ig = i, + ig, 
eg= Ey max sin (wt — 6,). 


Thus, with AZ = Othe transition from the rectifying condition with a positive half-wave of the supply 
voltage to the rectifying condition with a negative half-wave will be accompanied by consecutive short-circuit 


ing of the rectifying bridge diagonals. 


The free state currents and the time during which the diagonals are short-circuited, will be a function of 
AZ. But with small values of AZ and small stray fields, when the ratio of the resistance in the ac diagonal 
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Fig. 4. Graphs of currents for the third and fourth change in the 
operating conditions of rectifier circuits for the case under con- 
sideration, 


(including the resistance of the supply source) to that in dc diagonal does not exceed three, the time for which the 
bridges are short circuited will be limited to a small portion of a period, and the disruption of the rectifying ac- 
tion of the bridges will lead to an insignificant distortion of the sinusoidal shape of the currents in the circuit. 


Magnetic amplifiers respond to the direct component of the difference of the rectified currents instantaneous 
values, a component which can be found from the equation: 


(ig — ig) dot. (11) 


In the circuits under consideration, the rectifying action of the bridges during the transition from the positive 
to the negative half-period and vice versa is disrupted. for a small portion of the period only, and the direct com- 
ponent of the rectified currents difference is determined in the main by the real component of the equivalent 
stable state current difference Ai in the normal stable state rectifying operation of the bridges. This current dif- 
ference in its complex form can be expressed by the formula: 
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Fig. 5. Difference in the instantaneous values of rectified currents for the case 
under consideration 


For AZ = jAx the quantity Z,, in the circuits under consideration can be found from the expression; 


Za = Ro + Ra +2 + Roll + + — (13) 


The phase difference between the current Ai and supply voltage is found from the expression 


2RoRa + zy (1 —k) 
Az 
Ro + Ra + ry [Re + Ro (i + k,)] 


For AZ = AR, Z, is found in the circuits under consideration from the expression: 


2z 


The phase difference between current Ai and the supply voltage is equal to; 


+ [Ra + Re (1 + 


Ro + Ra + 


The differences in the instantaneous values of the rectified currents and the angle gA in the example under 
consideration are shown in Fig. 5, where the following notations have been adopted: ¢g is the half-wave of the 
supply voltage, AS» is the area formed by the curve of the stable state currents difference and the X~axis in the 
interval of the normal rectifying condition, AS, is the area formed by the curve of the current difference in the 
transient intervals and the X-axis, AS is the area formed by the curve of the difference of the stable state cur- 
rent and the curve of the actual current difference during the restoration time of the rectifying conditions in the 
bridges of both measuring circuits, ~; is the angle of advance of the auxiliary bridge diagonal yoltage, which 
determines the beginning of the transition from the positive to the negative half-wave rectification, $, + 
are angles of the total transition interval, gy is the angle of the phase difference between the supply voltage and 
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Fig. 6. Relation between the phase angle 
¢Y, and AZ = jAx for various inductances 
of the control winding. 
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Fig. 7. Relation between the phase angle 
%, and AZ = AR for various inductances 
of the control winding. 


the equivalent difference current Ai = ig— i, in the 
normal rectifying condition of the bridge circuit. The smal- 
ler the value of angle y,, the greater will the relative sensi- 
tivity of the magnetic amplifier be with respect to the vari- 
ations of resistance Ax (or AR). The value of angle gy de- 
pends on the parameters of the circuits. 


Figures 6 and 7 show the relation between angle gy, 
and the measured impedances AZ = jAx and AZ = AR for 
various inductance values of the magnetically coupled wind- 
ing. It will be seen from the diagram that the circuits with 
small values of inductance, have small angles g, with 
AZ = AR and large angles gy, with AZ = jAx, t.e., theyhave 
a larger relative sensitivity for variation of AR. 


Increased inductance leads to a decreased relative 
sensitivity to variations of AR and increased relative sensi- 
tivity to variations of Ax. 


Thus, in independent tracking of the quantity AR of 
impedance AZ = AR + jAx can be accomplished in circuits 
with a small induction of the control windings, i.e., with 
a small induction unbalance, The small value of inductive 
unbalance is easily attained by using additional short-circuited 


windings. 
Independent (uncoupled control of the quantity Ax 
should be carried out with large values of inductance in the 


control windings or with a negative feedback in the AR con- 
trol circuits, 
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POWDER CLUTCH ELECTRIC DRIVE AND ITS APPLICATION 
IN AUTOMATIC SYSTEMS 


G. F. Konovalov and Ya, I. Flid 


(Moscow) 


Investigation results of powder clutches with a solid filling are given. Expres- 
sions for static characteristics and transfer functions of clutches and clutch units are 
derived. A block schematic and a characteristic of a servo-system with powder clutch 
drives is cited. 


The automatic control system's speed of operation and accuracy standards, which depend to a certain ex- 
tent on the parameters of the drives, are continuously rising and require the development of new types of drives. 
One of these new drives is the powder clutch drive, which possesses many good qualities, The drive has a high 
ratio of torque to the moment of inertia of the driven part and a short duration of transients in the clutch. The 
transmitted torque in the powder clutch drive is a linear function of the control current, it has a high control ratio, 
small control power and a small weight per unit of the transmitted torque. Above qualities of the powder clutch 
drives give an idea of the possibilities of their application in automatic systems. 


The Torque of a Powder Clutch 


Investigations of powder clutches have shown that the relation of the torque transmitted by the clutch to 
the control current shown in Fig. I consists of three portions, In the first portion at currents i, < {,, the driving 
moment is a function of the control current and of the difference in the speeds of rotation of the driving and 
driven parts of the clutch, i.e., 


M = (nm — ng) + Ciiy, (1) 


akg -m where f is the coefficient of viscous friction, which depends on the geom- 
| etrical dimensions of the clutch, the composition of the powder and the 

material of the core, C; is the coefficient depending on the geometrical 
dimensions of the clutch and material of the core, n, and ng are the 
speeds of roation of the driving and driven parts of the clutch, 


With control currents thy > ly > fp, the driving moment of the 
ata © clutch is completely determined by the quantity iy. For this portion of 
6 540 5 O15 H4,ma the characteristic (Fig. 1) the driving moment transmitted by the clutch 
is determined by the expression: 


M =(C;,(iy — iy). (2) 


Neglecting stray fluxes, whose value,owing to the screening of the control winding, is small, it is possible 
to write; 


M => Ce (D, ®,), (3) 
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where C,. is a constant coefficient, , and , are the fluxes in the gap of the clutch corresponding to currents 
ty ond by. 

If the control of the clutch is exercised around some mean value of current {', expressions (2) and (3) can 
be written in terms of increments; 


AM =C,Aiy, AM =C,A9,, (4) 


where AM is the increment of the driving moment corresponding to the increment of the control current Aly, and 
Ay is the increment of flux corresponding to the current increment Aly. 
In the third portion of the characteristic M = f (ty) is limited by the saturation of the clutch core. 
Henceforth we shall examine the control of the clutch in the second portion of the characteristic M = f (ty) 
only, and will omit the notation A in the formulas (4). 


The Transfer Function of Powder Clutch 


The lag of the magnetic flux across the clutch gap depends on the induction of the clutch winding and on 
the eddy currents in the clutch core 


If the stray fluxes in the clutch are neglected, the effect of eddy currents and the winding inductance on the 
duration of the transients in the clutch can be represented in the electrical balance equation by an equivalent 
time constant T. In this case the equation representing the variations of flux © in the gap of the clutch will 
assume the form: 


d® 
®=0,+7%. 


Bu multiplying equation (5) by Cc we obtain: 
M=M,+T%%, 


Let us write down the equation of the driven part of the clutch: 


dn, 
Jz=M—M., (7) 


where J is the moment of inertia referred to the driven part of the clutch, and Mc is the moment of resistance, 


In servo and automatic control systems the drive is usually chosen in such a manner that the driving moment 
exceeds considerably all the braking moments. In such a case it can be considered that M, = 0. Having deter- 
mined from (7) the moment of the clutch and substituted it in (6) we obtain: 


where k = Cy/J is the clutch transmission coefficient, v = dn,/ dt is the acceleration of the driven partof clutch. 
For a stable state equation (8) becomes: 


v= (9) 


It follows from (9) that the acceleration of the output clutch axle varies proportionately to the control 
voltage, With Uy = const,the speed of rotation increases with time, Its maximum speed is limited by the speed 
of rotation of the driving part. Thus, characteristic ng = f (Uy) has the form of a relay characteristic. Applying 
the Laplace transformation to equation (8) and considering that v = d’p /dt*, where ¢ is the angle of rotation of the 
driven part of the clutch, we can obtain an expression for the transfer function of the clutch with respect to the 
angle of rotation: 


k 
= 


Double Connection of Clutches, Rotation Speed Control 


In order to obtain reversible movement,double connection of clutches is used, whose kinematic diagram is 
shown in Fig, 2, Clutches connected in this manner shall be called henceforth clutch units, Figure 3 shows the 
driving moment characteristic curve for a double connection of clutches, The clutch unit characteristic isbased 

on the initial coil currents igy = igg whose value is selected so as to 
. ve obtain a linear characteristic with small control currents, 


With My, > Mg, clutch M, Operates with a slip S< 1 and clutch 
Mg with a slip S > 1, Moreover the flux in the gap of clutch M, is 
larger than that in clutch Mg, Under these condition clutch Mg can 
produce a certain amount of damping and the transfer function of 
the clutch unit with respect to the rotation angle will be: 


ks 
pad qt) 


W (p) = 


where kg = C, / fs Tm=J/f, and fg isthe damping coefficient of the 
Fig. 2. Schematic diagram of a clutch unit, 

double connection of clutches, 
My, Mg are the clutches, M is the 
motor and n is the transmission 
axle. 


According to calculations and tests of a number of servo-sys~ 
tems, however, the damping coefficient of a clutch unit is very small 
and it can be assumed that in equation (11), for the frequency range 
of interest to us, T}4 >> 1. Then the transfer function for a powder 
clutch unit can be expressed with an accuracy sufficient for practical 


purposes as: 


ks, 
Ws(P) = (12) 


sion of g/sec**v. 


Since the powder clutch unit is a link with astatism of the 
second order, it is advisable for purposes of speed control to intro- 
duce a fixed feedback across the clutch unit, then the transfer func- 
tion with respect to the angle will be: 


W (p) 


Fig. 3. Relation between the driv- 
ing moment transmitted by the where kg is the transfer constant of the feedback circuit. 

clutch unit and the control current. By substituting equation (12) in (13) and performing simple 
1) For the first clutch; 2) for the operation we obtain: 

second clutch; 3) for the clutch unit. 


k 

where ker = 1/kg> transfer const, T,, = VT7kg Ky is the time constant and § ="/5(T/ Ty) is the index of oscillation, 


From expression (14) follows that tn the state (p © 
linear relation to the input voltage, i.e., 


n= Uin- (15) 
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The linearity of this characteristic is limited by the speed of rotation of the driving portion of the clutch, 


Experimental Investigation 


The object of the experimental investigation was to check the above expressions, to determine the time 


constant T of the clutch, to obtain the frequency and static characteristics of the clutch units with a fixed velocity 
feedback. 


Fig. 4. Clutch frequency characteristic 
with respect to the driving moment; 

2) the same for a clutch unit with res- 
pect to the speed of the output axle. 


Fig. 5. Block schematic of a servo- 
system. Ages is the error signal ampli- 
fier, F is the filter, CL is the correct- 
ing link, A is the amplifier, CU is 

the clutch unit, R is the reduction gear, 
TG is the tachometer generator, P is 
potentiometer and SF is the selective 
feedback. 


0 4 & 2 UE 


‘Fig. 6. Frequency characteristic of the 
servo-system. 1) The amplitude~fre- 
quency characteristic; 2) the error 
characteristic, 
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Investigations were conducted with cylindrical clutches 
filled with a ferromagnetic powder, consisting of a mixture of 
40partsin weight of ferrocarbonyl and one part of zinc oxide. 
The clutch parameters were determined at a constant speed of 
1,500 rpm. 


The relationship M =F (ly) for the clutch under investi- 
gation is given in Fig. 1. The characteristic determines the 
maximum driving moment M,, transmitted by the clutch and 
the control properties of the clutch with respect to its moment, 
The initial moment Mg is provided by the friction of the powder 
against the effective surfaces of the clutch and the friction in 
the bearings, its value does not exceed 1-4% of Mjyy. The 
hysteresis shaped characteristic {s not due to the magnetic 
properties of the core but to mechanical causes [2]. 


The clutch time constant was determined from the 
frequency characteristic of the clutch (Fig. 4) plotted against 
the driving moment over a linear portion of the characteristic 
M =f (i,) with a sinusoidal driving signal on the control wind- 
ing of the clutch, The instantaneous value of the moment was 
obtained by measuring the strain in the beam which held the 
clutch axle in a fixed position. The strain was measured by a 
bridge transducer, which had a strain gage connected in one 
of its arms. 


The frequency characteristic of the clutch unit was de- 
termined with a velocity feedback connected across the clutch 
unit, The amplitude-frequency characteristic of the clutch 
unit (Fig. 4) was determined as the ratio of the tachometer 
generator voltage amplitude to the amplitude of the fault 
signal and the phase~frequency characteristic as the phase dif- 
ference between these two voltages. 


A clutch unit with a velocity feedback ensures a large 
range for a smooth speed regulation. Tests have shown that 
the speed control ratio amounts to 1,000. Moreover the band- 
width of the drive is equal to 35-40 cps, which shows its good 


dynamic properties, 


A Servo-System with Powder Clutch Drives 


Figure 5 shows a block schematic of a servo-system 
with a powder clutch drive. The seryo-system has two paral= 
lel feedbacks, a fixed velocity feedback and a selective feed- 
back with a transfer function: 


pTi:pT2 
W 


& U 
| 


Moreover, in order to expand the bandwidth of the drive, a correcting filter with a transfer function given 
below is fitted; 


1+ 
W(p) = key 


Figure 6 gives the experimental amplitude-frequency characteristic and the frequency error characteristic 
with asystem quality factor of 250 sec! (full lines), When the quality factor was raised to 400 sec the band- 
width of the system at halfpower level increased to 25-30 cps, Dotted lines in Fig. 6 give the calculated values 


of the frequency characteristics. 
SUMMARY 


1, The paper provides expressions for a single clutch and clutch-unit transfer functions,which can be used 
for designing and investigating automatic control systems with power clutches, 


2. It is shown that the acceleration of a clutch or a clutch unit is proportional to the control signal. When 
a velocity feedback is connected to the clutch unit, a reversible speed control system is obtained which has good 
dynamic properties and ensures a large control ratio, 


3. Tests have confirmed the accuracy of the derived expressions, which represent the static and dynamic 
characteristics of the powder clutches. 


The application of powder clutches to automatic control systems provides these systems with good dynamic 
properties. 
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CURRENT STABILIZING DEVICES 


OE. Kubyshin and A. N. Milyakh 


(Kiev) 


LC devices are described which produce stabilized current in the load inde- 
pendently of variations in the load resistance, 


With the realization of various types of automatic circuits in electronic devices, in metering technology, 
and in other cases, the need arises to maintain automatically a constant value of load current independently of 
variations in the load resistance, To do this, ballast resistors and ballast lampsare used. These elements possess 
short lives and can not stabilize large currents. The efficiency coefficient of these devices, as well as devices 
which are based on the use of resistances and reactances in nonlinear regions is low. In addition, they are charac- 
terized by large time constants (with exception of the ballast lamps). Frequently, to get practically a nonfluc- 
tuating current, ballast resistors are used in series with a variable load. The efficiency coefficient in this case also 
is very low. The defects just enumerated limit the possibilities of applying the current stabil ization principle 


Fig. 2 


Fig. 1. T-circuit stabilizing device. U, = network voltage, Ug = voltage across 
the primary windings of the matching transformer, I}; = stable current with load of 
varying magnitude Zj; L, M, R = inductance, mutual inductance, and resistance, 
respectively, of the coil; C = capacitance of the condenser. 


Fig. 2, I-circuit stabilizing device. Designation same as in Fig. 1. 


At the Institute of Electrical Engineering of the Academy of Sciences of the USSR, LC devices with reso~ 
nance circuits were developed, which possess many useful properties, including that of stabilizing current in a 
load with widely varying load resistance. These devices do not have the drawbacks pointed out above, To re- 
alize these devices,the current stabilization principle was used as a basis on which to construct the idealized 
circuits of Boucherot, The devices built according to these circuits did not find wide application since they pos- 
sess low efficiency, bad power coefficient at the input, and excessively large dimensions, On the basis of de- 
velopment of Boucherot circuits and the introduction into the device of complementing elements and intercoupl ing 
between inductances, LC devices with very good characteristics may be worked out, 


| 
to automatic devices [1, 2]. 
! 
Matching 
transformers 
Fig. 1 
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In Figs. 1, 2, and 3 circuits of LC devices are shown, which were 
worked out for various powers and verified experimentally under laboratory 
conditions, The effect of stabilizing current with power is obtained by choos- 
ing the inductance of the coil and the capacitancesuch that at the working 
frequency of the alternating current, the capacitative and inductive resist- 
ances are equal, With this condition satisfied, the current in the load at 

the output of the device does not change if the voltage of the network re~ 
mains constant, Current at the input is obtained in proportion to the voltage 
across the load, thus preserving the required power output to the device, An 
LC device with a resonance circuit, which is connected between an alter- 
nating current network and a load maintains a stable current to the same ex~ 


tent as it does the voltage which drives the device, with accuracy within 
Fig. 3. Bridge circuit sta- ¥%. 


bilizing device. 


The introduction of inductive coupling between the coil windings 
makes it possible to reduce the capacitance, the dimensions of the coil, and in addition, to increase the effective 
Qof the coil and consequently, the accuracy in stabilizing currents by means of the device described. 


Proper choice of the resistances in the legs of the inductances makes it possible to lower the Q of the in- 
ductance to the value which is necessary to preserve the required frequency characteristics of the LC device. 


In the prototypes made for these LC devices, a change of load current between 0,2 to 0.3% correspond to a net~ 
work frequency variation of 1%. 


$ Network 
rac= 
e also 
r 
152 
EE 
435 
4 
_ Fig. 4. Three~phase bridge Fig. 5. Three-phase T~cir- 
. circuit stabilizing device. cuit stabilizing device. 
“ Proper choice of the resistance value makes it possible to construct devices to correspond to a given ac~ 
pos curacy for maintaining a constant current, 
; Placing into the circuit of the device a special transformer which matches the resistance at the output of 
upl ing the device with the load resistance which itself varies within wide limits preserves high power and efficiency 


coefficients, 
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At the Institute, theoretical and experimental investigations were carried out on devices for stabilizing 
alternating current, consisting of T- and II-circuits (Figs. 1 and 2), a bridge circuit (Figs. 3 and 8), and three- 
phase circuits (Figs.4 and 5), Devices for stabilizing constant current were also investigated, which were made 
from the above circuits with the addition of rectifiers (Figs. 6 and 7), as was a device for stabilizing load cur~ 
rent with varying input voltages (Fig. 9) for a power range of 51 w to 1.5 kw [3]. 


The performance of these LC devices was verified under working conditions, These devices were designed 
for supplying stabilized alternating current to the input terminals of track circuits, which must have small vari+ 
ation corresponding to variation in the leak resistance of the track circuits, The situation is that the leak resist- 
ance in these circuits is not very large, rising in good weather, falling many times in bad weather. In order to 
preserve the necessary value of the current which feeds into a track circuit in any weather and to protect the 
windings of the step-down transformer from overloading in the event of derailing, ballast resistors are used, which 
are periodically adjusted by service personnel. The application of ballast resistors to the supply and receiving 
terminals of track circuits leads to considerable power consumption for supplying 1 km of track circuits (up to 
150-250w/km). The LC devices, which are used for driving one-kilometer sections of track circuits (Fig. 8), 
worked well during the course of a year, maintained reliable signalization, and lowered the consumption of elec- 
trical energy per kilometer track circuit by a factor of 10 to 15. In this way, the application of LC devices to 
railroad eliminates the need for adjusting manually the currents at the supply end of the track circuits and signifi- 
cantly reduces nonproductive power consumption. 


LC current stabilizing devices may be used for driving trans- 
former-and autotransformer-coupled magnetic amplifiers and relays, 
| Matching magnetic triggers in electronic computing devices, various circuits 
See transformer in automatic devices for the normal operation of which ballast re- 
sistors are commonly used, It must be pointed out that because of 
the large power consumption in ballast resistors, noncontact trans- 
former-coupled relays, which are good in their own right, have not 
received wide application. In a unit with voltage regulators worked 
out either according to the ferroresonance principle or the principle 


of magnetically excited devices, LC devices may be used in the 
instrumentation technology and for various automatic devices in 
which a high degree of current stability is required. One of the pos~- 
sible means of making a load-current stabilizing device with varying 
input voltage is a magnetically excited matching transformer (Fig. 9). 


Fig. 6. Circuit for connecting a 
rectifier bridge to a single-phase The arrangement at the output of the semiconductor recti- 
current stabilizing device (Figs. 1, fier, which is connected with a bridge circuit (Fig. 6) or a three- 
2, 3). phase device connected to a Larionoy circuit (Fig. 7) makes it pos- 
sible to obtain a stable, rectified output current, In this way, it 
is possible to replace ballasts, which are used in constant current 
circuits ,to obtain a stable, constant current of great strength and 
power. LC devices with rectifiers may be used to obtain a stable 
current for galvanic electrolytic shops. 


In alternating current circuits, LC devices may be used 
with rectifiers in battery chargers to regulate current in charging 
storage batteries. According to weight, size, and cost, LC devices 
will not be widely uséd as battery chargers with rotary converters, 
but, due to reliability of performance and the possibility of regulat- 
ing within a wide range the value of the charging voltage and out~ 
put current, they will significantly surpass them (the battery chargers). 
In addition, the voltage which is needed for charging will be set 
automatically, depending upon demand,without a system of special 


regulating devices, 
Fig. 7. Circuit for connecting a rec- The dimensions of LC devices and their weight-to-output- 


tifier bridge to a three-phase current power ratio are, at the present, small. They are characterized by 
stabilizing device (Figs. 4 and 5). the following examples: the device which is supplied with three- 


Larionor 
circuit 
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Matching 


track circuit 


Fig. 8 Fig. 9 
Fig. 8. Bridge circuit stabilizing device for driving the track circuits in railway 
signaling. 


Fig. 9. Principle circuit of a magnetically excited matching transformer for cur- 
rent stabilization with fluctuating network voltage (U, = var). The stabilizing 
part of the circuit is indicated by the ballast resistor Rj, 3 Ro, Re are the resistances 
for a set of values of the stabilized current ly W.., Woc, ald Woy are the control 
windings of the magnetic transformer, Ry, R}, and Rf are the rectifier bridges, 


phase current with frequency 50 cps and network voltage 220 v, gives the possibility of obtaining 1.5 kw power with 
a stabilized current of 4amp. The weight of the device is approximately 40 kg/kw. 


The devices described are extremely simple in construction, They consist of a small quantity of reliable 
components; coils with air gaps, resistors, capacitors, and matching transformers, The transformers and coils are 
calculated and manufactured according to the nature of the device, but capacitors and resistors may be chosen 
for available lists. 


The analysis of LC devices by a method which has been worked out is not involved and giveshigh accuracy, 
All the properties pointed out above and other properties, to be brought out below, are a consequence of 


minimizing the self-impedance at the input and output. For all practical purposes, there remains only the mutual 
resitance zm, and, as a result, the behavior of the device may be described approximately by the equations 


4 
U, = 2mls, I, = (1) 


where U, and Ug are the voltages, I, and I, are the currents at the input and output, respectively. 
Equations (1) make it possible to determine the property of LC devices of transforming impedances in the 
following form: 
gt 
a= (2) 


In this way, the resistance which is connected at the secondary terminal of the LC device is transformed 
to the primary terminal as admittance, the value of which is proportional to the load, A break in the secondary 
circuit is equivalent to a short made at the secondary side equivalent to the open state of the primary circuit. 
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The relations introduced show also that an inductive load is transformed by the device into a capacitive 
load, while a capacitive load is transformed into an inductive load. A linear or nonlinear inductance is trans- 
formed proportionally to the linear or nonlinear capacitance, respectively, and conversely, An LC device may 
be used to convert a continuously regulated inductance into a continuously regulated capacitance, and vice versa, 
This property may be used to obtain noncontact and regulated inductances with the aid of variable capacitors 
(variometers) and to obtain by means of variometers or magnetically excited coils continuously regulated capaci- 
tances (to get regulated capacitors with large capacitance). 


For the proper choice of components in an LC circuit for current stabilization with an active load at the 
output of the device (Z,,; = Ry) the input power fact « 1, The inductive resistance at the output of the device 
(Z.44 = WL4y) decreases the input power factor and makes it negative, converting the device itself into a capacitive 
load for the network, 


The investigations of LC devices with resonant circuits have demonstrated the expediency of applying 
them in various automatic devices, most of all due to their property of stabilizing load current for variations 
within wide bounds of the load resistance, The methods which have been used up to this time for stabilizing cur- 
rent, are extremely uneconomical; they do not allow stabilization of very large and very small currents for driv~ 
ing large and small loads and are not sufficiently reliable. 


It is natural that the defects enumerated above have severely restricted and still restrict application of the 
current stabilization principle to automatic devices. Electrical engineers have become accustomed to dealing 
with a constant voltage (more or less stable), dismissing as not worthwhile those possible solutions of technical 
problems which might have advantages for use as energy sources, but which work with constant current (current 
sources), 


In this fashion, LC current stabilizing devices fill a considerable gap in the means and possibilites of 
automation. They may be considered as converters of a voltage source into a current source which do not con- 
tain moving parts or contacts and which, therefore, work for an indefinite period of time. It is possible to con- 
struct these devices practically for any power and for any current strength. 
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A STABILITY CRITERION FOR NONLINEAR CONTROL SYSTEMS 


Chang Ssii-ying 


(Moscow) 


1, Posing the Problem 


Let us examine a system of differential equations which represent a disturbed movement of an automatic 
control system, 


m 
Me = + (k = 4, 2,..., m), 


Vea + Wit Sp = f*(e), 


m 
» Pate — 


In order to find a stability criterion let us reduce system (1) to A, 1. Lur‘e’s [1] canonical form: 


— + (0) (k= 1, 2,..., 844), 
n+1 


om tate 
kel 


. 
a= >) r’f(o), 


Rel 
where p,, is the root of the determinant: 


Next with Rep, > 0, A. I. Lur’e obtained Lyapunov'’s function [1] 


n+in+1 
= 


and found the conditions of stability 
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1 
B, +2Vr'a, + 2a, 


(k= 1,..., 4). (5) 


The criterion of stability is represented by the inequalities 
Fy Ba gy) (k= 4, 2,...), (8) 


whose implementary solves equation (5). But with k > 3, the development of the inequalities (6) becomes very 
difficult. 


For the case when P,, are real quantities, A. M, Letoy [2] obtained another Lyapunov's function 


n+1 
i 
k=1 
and found the criterion of stabflity 
n+1 
(1+ RB,) 
(8) 


where R is a positive constant. 


2. Additional Investigation of Criterion (8) 


Let us examine the inequality (8) where R is a positive arbitary quantity. Each value of R which satisfies 
(8) corresponds to a criterion of stability. If in (8) we change to sign of inequality by a sign of equality, we ob- 


tain the expression: 
2 n 1(4 + RB,)* 


which represents a boundary of the stability region, 
Let us introduce new variables: 
(B,, Per’) y= y (By, 


4 so that (9) assumes the form: 
F (z, R) = 0. 


' This form represents a family of curves of plane x, y for which R is a parameter. For such value of R which 
satisfies (8), there corresponds a definite curve of the family of curves, which in turn determines a region of stability 


in the plane x, y. 


If there is an envelope [3] which embraces all the regions of stability, there will be one value of R which 
will correspond to that envelope. By substituting this value of R in (8) it is possible to find one criterion of sta~ 
bility which will include all the criterions determined by the inequality (8). 


Let us prove the existence of the envelope. 
Developing (8) we obtain: 


Sant Pe (9) 

k=l Pk (10) 
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and rewrite (10) in the form 


Let us introduce new variables x, y so that this inequality has the form: 


(y+ R)>M, 
Where M {is a constant and x and y are determined by the equations 


In the case of equality in (11) we have: 


R)(y+ R)=M. (12) 


This is an equation of a hyperbola. | 
Let us derive an equation for a tangent through an arbitrary point (a,b) to this hyperbola, This equation 
has the form: 
y—b=— (ea). (13) 


If a point on the tangent very close to point (a, b), does not lie in the region of stability, then the region 
on the concave side of the hyperbola is a region of stability. In this case the required envelope exists, In the 
contrary case it does not exist. 


Let (x, yy) be points on the tangent where 


Mh 


and h is a very small quantity, 
After simple calculations we obtain: 


(21 — + R)= M 


where h?/(a — R)* is always positive, 
Hence 


By comparing the last expression with (11) we find that point (x4, y;) lies outside the region of stability and 
hence the envelope exists. 


Now let us calculate R, Let: 


Let us denote: 
o+ 198 n 4 
— RU— Rm—s>0. 
| 
aF (2. y, R)_ 
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n+i1 


k=l 
R= n+1 


D> (82 /ex) 


Substituting R in (8) we obtain a single criterion of stability: 


3. Examples of Applying Criterion (14). 
A. Bulgakov's First Problem [2]. Let us reduce the system of equations: 


= f* (0), 


= f(s), 
— pets + f (0), 
= + Bars — f (0). 


In this case p; = 0, r*° = 1 and k = 1,2. 
According to (14) we obtain the criterion of stability: 


Ba 
Let us express it by means of the initial parameters 


>e. 


This result coincides with that of Lur'e, Bulgakov, Andronov and Bautin. 
B, Bulgakov's Second Problem [2]. Let us reduce the system of equations 


= f* (6), 


= — pits + f (0), 
— pats + f (0), 
= Biz + Bare — f (6). 


From (9) we find: 
2 
ay 
\Pk 
to a canonical form: 
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4) 


According to (14) we obtain: 
Since pypy > 0 and 8 )* > 0, we obtain: 


1—(E +P) 50, 


By substituting Y* in (15) we have: 
(Bs — Ba) V + 2 >0. 
Inequalities (6) in this problem have the form [2} 


Gab) 


It will be seen that the first inequality coincides with (16) and the second has certain differences ‘rom (17). 
The application of inequality (14) shows that the calculation of the stability criterion by this method ts 
much simpler than by A, I. Lur’e's method, 
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TRANSISTOR COMPARISON DEVICES 


Vv. M. Polyakov 


(Moscow) 


Comparison devices of the null~element type are considered which are con- 
structed using transistors and which have high reliability, small dimensions, and low 
power requirements, 


Electronic comparison devices find wide application in many branches of technology. The purpose of these 
devices is to obtain an output signal at the moment at which a variable input voltage reaches a definite given 
level. This voltage level may be zero or close to zero in magnitude, In the latter case, the comparison devices 
are called null elements, 


_ Null elements find use in automation and remote control and, in particular in telemetering technology, for 
example in multicomponent devices, which operate by means of remote pickups, etc. 
In apparatus of this type the use of a large number of null elements is sometimes required. In these cases 
the null elements must have the smallest possible dimensions, high reliability, and low power requirements. 


_. Vacuum~tube comparison devices cannot completely satisfy these requirements for many reasons — in particu- 
lar, because of inadequate reliability of vacuum tubes themselves. 


Input 1 


Fig. 1 


The use of transistors for these purposes makes it possible to obtain null elements which satisfy the require- 
ments listed above. It is possible to obtain such comparison devices either by partially replacing vacuum tubes 
with transistors in known circuits or by working out new circuits, using only transistors, 


As example of a comparison device in which vacuum tubes have been partially replacee by transistors, Fig. 
1 shows a circuit of a device which is a variation of the tube circuit of Fig. 2. 


As is evident from Figs, 1 and 2, input to the circuits is to two stages of comparison devices, which work in 
parallel on one load, The principle of operation of the circuits is illustrated by the oscillograms shown in Fig. 3. 
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Fig. 2 


The oscillograms describe voltages at separate points of the circuits for a sinusoidal input voltage of amplitude 
5 v, frequency 50 cps. 


The principle of operation of these circuits is as follows, During the negative half-cycle of the grid voltages 
of tubes T,, at which time these tubes are cut off, in the circuit of Fig. 1, the charge of condensers C, and Cy 
(Cy and C, in the circuit of Fig. 2) rises), The condensers build up charge to a voltage close to the plate power 
supply of tubes T, (6 v). By varying the voltage on the grids of tubes T,, at the moment of their operation, condensers 
Cy and Cg (Cy and C,) discharge ~ a process which is accelerated due to the action of a strong positive feedback 
coupling, which is obtained with the aid of the amplifier T, in the circuit of Fig, 2 and transistors 1 and 2 in that 


of Fig. 1. The input impedances of transistors 1 and 2 in the circuit of Fig. 1 are used as resistive differentiating 
networks with condensers C, and Cy. 


Transistor 3, which is connected at the output with transistors 
Va 1 and 2, produces summation of pulses from the separate stages of 


beh bf te the comparison devices, thus providing full independence in their 
operation, 


The large amplification in transistor amplifiers, with the 
same power supply, makes it possible to obtain better results with 


We ia Ww api Ow a, the circuit of Fig, 1 than those obtained with the circuit of Fig. 2. 


7 For example, a voltage pulse taken at the output of the 
, ae circuit of Fig. 1 has period 2~3 psec with amplitude 6-7 v; in the 
! circuit of Fig, 2, the period is 10-15 ysec with amplitude approxi- 
Fig. 3 mately 0.05 v. 


The circuit of Fig. 1 operates stably with input voltages of 
1 v and higher, frequency 10 to 700 cps, and filament voltage for 
f, : =20v tubes T, of 2.5 v. 


Input 1 . A comparison device assembled according to Fig, 1 has 
smaller dimensions than one using the circuit of Fig, 2; it is also 
Output more economical in power requirements. 


The common defect of the circuits of Figs, 1 and 2 is the 


ey presence of a large number of components, This , as well as the 


Yb 


presence of vacuum tubes, leads to lessening of reliability and 
1 increases in size of the devices, 


Fig, 4 These defects are absent in comparison devices assembled 
according to the circuit of Fig. 4, 


The operation of the circuit of Fig. 4 is based on the fact that transistors 1 and 2 start conducting for a defi- 
nite value of the input voltage, This leads to a rapid voltage drop across the inductive loads Ry and Ry. This sud- 
den change of voltage leads to forced oscillations in the circuit consisting of the primary windings of the pulse 


Input 2 
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Fig. 5 


transformer Tr, condenser C, (or C3), and the collector resistance — the transistor base, The frequency and ampli- 
tude of the forced damped voltage oscillations may be determined by the parameters of this circuit. The period 
and amplitude of single pulses separated from these oscillations may also be determined by the circuit parameters. 


With the circuit parameters indicated in Fig. 4, and with the sinusoidal input signal of amplitude 5 v, fre- 
quency 50 eps, output pulses were received of period 3-5 psec and amplitude 0.5 v. 


The circuit normally works in a temperature range of —50 to + 45°C. 


An example of a multistage circuit based on the circuit shown in Fig. 4 is given in Fig. 5, which shows a 
ten-stage device, 
As is evident in this circuit, the transformers are not connected to more than 3 to 4 transistors. 


Such a restriction is explained by the shunted opeiation of the transistors, which decreases the pulse ampli- 
tudes. 


The advantages of the device of Fig. 5 are its small dimensions (50 x 70x 70 mm) and low power require- 
ment, The current required does not exceed 1 ma per 20 v collector voltage, 
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INDUCTION TACHOMETER AS AN ANGULAR ACCELERATION 
TRANSUCER 


S. T. Kazaryan 


(Moscow) 


An angular acceleration transducer in the form of an induction tachometer 
with de excitation is examined, 


Induction tachometers are widely used in industry and a large literature has been devoted to their study [1-7]. 
In automatic control circuits, however, the induction tachometer can find a new application as a transducer of 
angular acceleration generating a stabilizing signal, 


In fact, if the excitation winding of an induction tachometer is fed with a direct current, and if the speed 
of the rotor remains constant, the flux will remain constant and there will be no emf induced in the transverse 
winding. If, however, the speed of the rotor varies, the flux will change and its transverse component will induce 
an emf in the transverse winding. The value of this emf will bedirectly proportional to the angular acceleration 
of the rotor. The sign of the emf will depend on the direction of rotation. 


The magnetic field of the induction tachometer 
arises owing to the magnetizing force of the stationary 
windings of the stator and the magnetizing force of the 

* rotor currents, This field can be divided into two com- 
ponents — the longitudinal and transverse, whose axis 
coincide with the electrical axis of the tachometer wind- 


ing. 


In order to be able to conduct investigations inde- 
pendently of the number of poles in the machine, the 
angles, axis, angular velocities and speeds of rotation are 
everywhere taken as electrical, which differ from the 
mechanical ones by a factor of p, the number of poles in 
the machine, 


Only for the two-pole machine will these quantities 
coincide with the actual values. 


If the saturation of iron and hysteresis be neglected, 
there will be a linear relation between the magnetizing forces and the corresponding fluxes, 


It is possible to deduce balance equations for the magnetizing forces both of the longitudinal and transverse 
fluxes and balance equations for voltage of the winding of the machine, 


For the excitation winding placed along the longitudinal axis Od (Fig. 1), we have: 


Fig. 1 


di, 
Lea ar t+ * wae a 
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For the transverse (right angle) winding placed along the transverse axis Od we have: 
di 
a 
(rq + ig + (Lyq + Ly) =— ky = eq: 


Here rd, tq and r,, are the resistances of the stator windings and the load, Igq and Ls, are the stray inductances 
of the stator windings, Ly is the inductance of the load, wq and Wq are the number of turns in the stator windings, 
kwq and kwq are the winding coefficients, @q and @q are the longitudinal and transverse fluxes, ¥q is the full inter- 


linkage of the excitation winding, y and ig are the winding currents, ug and e, are the supply voltage and the emf 
induced in the transverse winding. 


For the longitudinal magnetizing forces we have: 


where Rq is the reluctance to the longitudinal flux @q and Fgp is the longitudinal magnetizing force of the rotor 
which is equal to the sum of two magnetizing forces, the transformer one and the rotational: 


Fap =F ap 
For the tranverse magnetizing forces: 


i 
where Fpgq is the transverse magnetizing force of the rotor, which is also equal to the sum of two forces: 


=F ap p. 


In order to determine the emf e, induced in the transverse winding it is necessary to establish the law of 
variation of the transverse flux: @q = q(t). 


The above equations determine the flux $g in the most general case, An exact solution of these equations 
presents considerable difficulties. Hence it is advisable first to simplify the problem, considering the actual oper- 
ating conditions of an induction tachometer. 


In the first place under normal operating conditions the induction tachometer has such a load that the second- 
ary winding currents can be neglected. Its operating condition approaches that of an open circuit: 


di 
i, =O and = 0. 


In the second place,when a direct current source is connected to the excitation winding and its voltage 
ud = Ug does not vary, the current flowing in the excitation circuit remains practically constant with the actual 
values of acceleration of the rotor: 


Ue di, 


in the third place, with varying fluxes both in the rotor and stator windings transformer emfs will arise. 
However, with the acceleration encountered in practice the transformer emfs arising in the rotor will be small 
compared with the rotational emfs,and the former can be neglected, The same applies to the transformer emf 
induced in the excitation winding, since it will be very much smaller than the supply voltage. This condition 
has already been expressed by ig » Ig, but now we can also add the expressions: 


_, 


ual 


ll 
nf 


and for the magnetizing rotor currents: 
The emf in the transverse winding by the variations in the rotor speed can now be easily determined if all 


the above conditions are taken into consideration, According to the above equations and conditions we can now 
write: . 


ralg= Uy 
(1) 
d®, 
(2) 
U4 


(4) 
Fau= Fey = Fapp 


The transverse flux @q produces in the gap between the stators a sinusoidal distributed flux density equal 
to the expression below if the harmonic components of the magnetic field are neglected: 
Bae = Bay Sin E, 


where Bdm is the amplitude of the flux density along the longitudinal axis Od (Fig. 1). 


Let us assume that the hollow rotor cylinder in the limits of its effective length 1 consists of an infinite 


number of adjoining elements, rods, placed along the generating line of the cylinder, The longitudinal field will 
induce in these elements a rotational emf equal to: 


where q = IDBgm is the longitudinal flux. 


If the rotor stray flux is neglected, which is permissible for an actual induction tachometer with a thin-walled 
elongated hollow rotor, the resistance of a rotor element will be: 


kyl kyl 
as ™ D 
a 
Here dS =A(D/ 2) d§ is the cross-sectional area of the rotor element, shown shaded in Fig. 1,4 is the thick- 


ness of the rotor, D/ 2 is the mean rotor radius, ky = 2a/I> 1 is a coefficient accounting for the length?, of the 
actual rotor currents path being longer than the effective length of the rotor, 


The conductivity of a rotor element will then be equal to 


AD 
= 


where gy = yAD/2kjz! is the conductivity of a portion of the rotor, bounded by a double-sided axial angle of 1 
radian, 


The rotational emf edp,R will produce elementary currents in the rotor elements, Assuming that these 
currents are directed in the effective part of the rotor along the generatrix of the rotor cylinder, we shall obtain 
the following expression for the currents: 


| 

= 
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| 
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di, PR = godee, PR? = F sowD, sin 


These currents produce a magnetizing force standing wave: 


The transverse flux ©, produces in the gap of the induction tachometer a sinusoidal distributed flux density, 
which establishes a standing wave in the rotor magnetizing force: 


2 
FapR = \ dign R = — 


Substituting in (3) and (4) the magnetizing forces Fqp p and Fay pby the expression obtained for them and 
eliminating we obtain: 


Uy (5) 


Here kwq/Rq and k,q/R, are the permeances along the axes Od and Oq, whichtretaken to be equal to G 
for a symmetrical machine. 


The emf induced by this flux is: 


— 4 


d® 


After substituting w* and w' = dw/ dt by n®and = dn/dt, and assuming = 0,79, we obtain the 
design formula for the emf in the transverse winding of an induction tachometer with dc excitation at any rotor 
speed: 


A An? —i 
where 
w? 
A = = 11-10-%g3G? (1) 
is a constant, 


We shall obtain expressions for the fluxes in the form: 


wyGVAn (8) 
7,’ 


i U 
(9) 


= \ dignn = 
Go kw kw &o 
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These fluxes for the given machine only depend on n, 
With an increasing n the longitudinal flux 4 decreases con- 
Region n< 0 tinuously, The transverse flux %q increases at first and at a 
speed determined by the equation 04,/0 n= 0, ie., atn= 
= +Vi7A = £9,5/ g9G it reaches a maximum, which is 
0 equal to half the longitudinal flux $q at n = 0, and then 
Yamal One decreases continuously. 


If gq and , are represented conventionally in the form 
Region n > 0 of vectors,the resultant flux: 


q 
On 


Fig. 2 will equal precisely the geometrical mean of fluxes g and 
%q,.max, Hence the locus of the end of the resultant flux @ 
vector is a circle (Fig. 2). The upper portion of the circle in Fig, 2 correspqnds to fluxes with a changed direction 
of rotation of the rotor, i.e., a changed sign of n. 


At high rotor speeds when n -> 00 or at finite values of n, but a superconducting rotor (gy = «0, A = oo) the 
rotor field completely compensates that of excitation and the resultant field tends to zero, i.e., the rotor screens 
the excitation winding field, 


6H 


Fig. 
Fig. 3 gives the curves of fluxes and #q,the emfs and + €q, plotted for one machine, according 
to calculations by means of above formulas, The curve for emf eq is based on formula (9) for flux 4, The same 


figure gives an oscillogram for emf eq for comparison with the calculated curve,” Curves and oscillograms have 
been obtained for the same value of; 


4000/2 — = const. 


* In Fig. 3 for convenience of plotting and comparison, a curve and an oscillogram of —e, ,1.¢.,¢q with a reversed sign 
— , are given instead of = 
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Theoscillations of the emf €, in the oscillogram are caused by the asymmetrical position of the rotoreach 
oscillation corresponds to one revolution of the rotor, By adjusting the position of the bearings these oscillations 
can be reduced to a minimum. 


The increase in the oscillation amplitudes eg at a certain speed n 1s due to the resonance, By changing 
the resonance frequency of the mechanical oscillations of the rotor, for instance by increasing the diameter of its 
axle, it is possible to shift the resonance outside the operating frequency range of the rotor, 


The coincidence of the theoretical curve of the emf eg with its mean values on the oscillogram is easy to 
see, Moreover the coincidence is not only qualitative but also quantitative. 


When an induction tachometer is used for acceleration measurements in the ideal case the emf eg should not 
depend on the speed n but only on the acceleration n', moreover, the relation should be linear, The relation between 
and n° is linear, but a complicated relationship between e&q and n places certain conditions on the use of an 
induction tachometer as an indicator of acceleration. These aio are determined by equation (6) or in a 
more graphic manner by the curves in Fig, 3, According to these curves, an induction tachometer can serve as 
an acceleration indicator in the following cases, 


a) If the values of conductances gp and G are large, the critical speedn at which eg = 0 is small, then 
permissible speeds exceed the critical, The inevitable errors caused by the relation between eg and n depend on 


the shape of the curve eq = eq(n) in the operating frequency range of n. 


b) If the values of gy and G are small, the critical speed n at which eg = 0 exceeds the values of the oper- 
ating range frequencies of the machine, Here any frequencies are permissible from n = 0. In this case the 
values of e, are sinall and the errors large. 


A successful solution of the problem may be obtained by taking into consideration the actual operating 
conditions of the tachometer in the circuit, 


In conclusion, it should be noted that when the tachometer is dc excited it is not necessary to have a trans- 
verse secondary winding, i.e., winding at 90°,out of phase with respect to the excitation winding. The wind- 
ings can be placed at any angle and even with their axes coinciding and then the emf will vary according to the 
curve eq in Fig. 3. 


This circumstance provides wide possibilities for constructional variation and simplifications and for design- 
ing new types of machines with the de excitation, 
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SIXTH YOUNG SCIENTISTS* SCIENTIFIC AND TECHNICAL CONFERENCE 
‘OF THE USSR ACADEMY OF SCIENCES AUTOMATION AND 
REMOTE CONTROL INSTITUTE 


The sixth scientific and technical conference of young scientists of the Institute of Automation and Remote 
Control of the Academy of Sciences of the USSR on problems of automatic control, dedicated to the opening of the 


21st Congress of the CPSU, was held on January 19-21, 1959. At this conference, organized by YCL Committee 
of the Institute, work completed in 1958 was discussed. 


The ‘conference was attended by more than 300 people, including representatives of 37 organizations of 


Moscow and the Moscow region, members of higher educational organizations, scientific research institutes, develop- 
ment and design organizations and factory personnel, 


At the conference 65 papers and 11 reports were read which dealt with a wide circle of scientific and tech- 


nical subjects in the sphere of automatic control. The papers aroused great interest at the conference, Over 50 
people took part in the discussion, 


In his opening remarks,the director of the Institute of Automation and Remote Control of the Academy of 


Sciences, USSR, V. A. Trapeznikov stressed the importance of solving problems of automatic control for the develop- 


ment of our national economy and called upon the young scientists to link up closely their work with production 
problems and to introduce intrepidly new technical ideas, 


At the plenary session a paper was read by Prof, V. S.Pugachev, D.Sc, on the “Basic problems of the informa- 


tion theory and the theory of statistical solutions and their significance for the development of the theory of auto- 
matic control”, 


The paper dealt with three types of automatic control problems which require statistical methods for their 
solution: the problem of filtering out signals from a noise background, problems of transmitting signals along 
channels with interference, and the problem of designing self-adjusting systems. It was shown that these three 
types of problems, differing in their posing, complexity and technical applications have a common basis, The 
author noted that all these problems represent a general type of problem of the theory of statistical solutions de- 


veloped hy Vold. The solution of the above problem opens up prospects of deriving a statistical theory of the pro- 
cesses of “self-tuition® in machines, 


Further work of the conference was conducted in sections; I— automatic control, I] ~ automatic checking, 
III — computers, IV — components and devices in automation and remote control, V ~ the theory of transmitting 
messages, VI— the theory of relay operated devices, and VII — the theory of automatic drives, 


l, Automatic Control Section 


The papers of L, I, Rozonoér and A, G, Butkovskii dealt with the general theory of optimum systems, In his 
paper, L, 1, Rozonoér showed that B. V. Bulgakov's problem of the accumulation of deviations in a dynamic system 
can be reduced to a variational problem of the type considered in the theory of optimum systems, L, S, Pontryagin's 
maximum principle provides a solution of the problem on the accumulation of deviations for a wide range of linear 
and nonlinear systems, This principle is applicable in all problems of theory of dynamic accuracy where it is as- 
sumed that the control system is subjected to external interference, information about which is limited, 


Having applied the classical method of infinitely small variations, A. G, Butkovskii proved L. S, Pontryagin’s 
maximum principle for the case when the controlled system is represented by a differential equation of the type: 
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f(x, 2’,..., + u (2). 


A large proportion of the papers dealt with self-adjusting and extremal systems. 


V. A. Ivanov's paper dealt with self-adjusting systems of automatic control of electric pipe-welding machines, 
On the basis of an analysis of the work of such machine the author made a choice of the best system for control- 
ling automatically the operating conditions of welding and found a design and a circuit for the compensating and 
correcting computers, which maintain the welding temperature at a given value, The results of the study of 
self-adjusting processes in the presence of sudden and harmonic disturbances were cited and test results appended. 


In S, A, Doganovskii's paper the results of analysing a self-adjusting automatic control system of the thick- 
ness of rolled metal on an electronic analog were given. Compensation for the rolled metal thickness variations 
was made by a special computer, The correction of the parameters of the computer was made by means of a 
special optimizing device. 


N, V. Grishko's paper dealt with determining the optimum (in the sense of the maximum efficiency of 
an engine) law of supplying fuel to internal combustion engine according to its load. Optimum operation was 
attained on the basis of the mean low pressure value in the suction pipe, since it was shown that minimum pres- 
sure corresponds to maximum efficiency. 


A. G. Shubin's paper dealt with control unit circuits of multi-channel automatic optimizers and provided the 
specific requirements for the control units, The author gave a number of concrete recommendations and a circuit 
for the control units of automatic scanning systems, 


A. V. Kalinina examined in her report electronic circuits of a 12-channel automatic optimizer operational 
unit and electronic relay circuits of an operational unit of a 10-channel automatic optimizer. 


E, M, Braverman provides in his paper formulas for calculating small variations in adjustment, which im- 
prove transient processes for industrial control systems, Consecutive application of these variations in adjustment 
with the determination each time of the actual position on a complex plane of the three roots nearest to the 
imaginary axis leads to an optimum adjustment of regulators, Several papers dealt with the theory of pulse control, 


In his paper, entitled "Automatic control systems which attain a steady state of operation in a finite time" 
Wang Sing-ming shows the advisibility of using filters in time delay units and provides a method of calculating a 
link with a transfer function inverse to one of a noncompensated closed system, 


V. I, Gukov's paper deals with a method of correcting linear systems by means of discrete filters. A correc- 
tion of transient processes of any system can be easily obtained by means of discrete filters which produce com- 
binations of transient process. derivatives. This method of correction and differentiating could be applied for cor- 
recting linear open line systems, moreover, the transient process of the corrected system stabilizes itself precisely 
in the specified time. 


I, V. Pyshkin's paper entitled "Time delay compensation in pulse-width automatic control systems" dealt 
with a control system consisting of a series-connected width-modulated pulse element, a regulator, a delay link 
and the object of control, The idea of compensation consists in forming by means of an extrapolator and discrete 
filter a signal at the input of the pulse element which would have been there in a system without time delay, The 
paper included formulas for calculating parameters of such devices. 


In his paper, M. M. Simkin dealt with periodic operating conditions in systems using digital computers with 
level quantization, gave the conditions for the existence of periodic states and provided an approximate solution 
to the problem of finding the periodic states, 


L, A. Gusev examined a method of determining periodic modes of operation in automatic control systems, 
which contained elements with random and piecewise-linear characteristics, The periodic states were determined 
by the author accurately in the form of complete Fourier series without neglecting the harmonics. 


In his paper, £, M. Solnechnyi attempted to develop the general synthesis of a multilink control system, The 
author examined the problem of designing a control system which would provide the required transfer function from 
each of the considered disturbances and each controlled coordinate. The paper provides the required transfer func- 
tion from each of the considered disturbances and each controlled coordinate. The paper provides the required 
conditions for making this problem soluble. 
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V. S. Kislyakov's paper dealt with certain problems of the theory of delayed automatic control represented © 
by linear differential equations with a delayed argument, 


In N. P. Kozlov's paper the dynamics of a hydraulic oscillating circuit were examined by means of a phase 


plane. The hydraulic vibrator examined by the author served for booster stabilization and vibrational linearization, 


1, A, Mikhlin dealt in his paper, entitled "Automation of the evaporation equipment of a sugar plant*, pro- 
vided the results of investigating typical evaporation installations with the object of automation, The author 


examined the problem of finding the optimum solution from the point of view of using the minimum amount of 
steam with maximum productivity. 


Il. Automatic Checking Section, 


The first session of the section examined papers dealing with automatic checking of composition and 
quality. 


In a paper entitled “Checking of the resistivity and geometrical dimensions of products made of nonmag- 
netic metals by means of eddy currents", V. P. Grabovetskii examined the basic relations between resistance intro- 
duced by the coil and the parameters of samples made of nonmagnetic. metals of cylindrical and prismatic shape 


placed in the field of this coil. On the basis of this analysis, the author recommended circuits for devices using 
this method of checking. 


In a paper entitled "Contactless indicator of the concentration of two-compoment mixtures", A, A, Desova 
dealt with a contactless instrument for determining concentrations of solutions, by measuring their electrical con- 


ductivity by means of eddy currents, The paper proposed a system of temperature compensation of the instrument, 


In a paper entitled “Checking of the composition of multi-component systems", V. N, Skugoroy examines 
methods of determining and controlling concentrations of multi-component mixtures on the basis of measuring 
their physical and chemical properties. 


In a report entitled “Certain methods of increasing the speed of operation of mixers for physical constant 
analyzers of multi-component mixtures", B. A. Pereverzey dealt with the development of a mixer for measuring 
and controlling parameters of solutions. The new mixer has a number of advantages as compared with the exist- 
ing types. 


The two following papers dealt with the analysis of composition by radioactive methods, 


A. A. Naumov dealt in his paper with the theoretical basis of using re-radiation of beta-radiations for the 
automatic control of complex processes, The theoretical conclusions were checked by experiments, The author 
recommended the design of automatic control instruments based on this method, 


In his paper, N. A. Pivovarov examined the possibilities of using the absorption of gamma-radiations for 
automatic checking of the composition of multi-component media, The paper included the theoretical basis 
of the method, the results of preliminary tests and practical recommendation for the implementation of the method, 


The second session of the section examined papers dealing with the development and investigation of auto- 
matic control systems, 


The paper by I, N, Panasenko, M, B, Rybashov and I, A, Tsaturova on “An automatic potentiometer with a 
dynamic correction of primary converters" dealt with the very important problem of automatic checking and cor- 
rection of dynamic transducer errors, The paper dealt with the technique of correcting dynamic errors of primary 
converters, approximated by a linear link of the first order and described the automatic potentiometer circuit in 
whose feedback was connected an operational amplifier, This potentiometer lowers considerably the dynamic 
errors of temperature measurements made by means of themocouples, In the discussion of this item, suggestions 
were made for speedy production of this potentiometer, 


In her report entitled “Decreasing the time constant of an output device of a mass-spectrometer™, I, A, 
Tsaturova recommended the selection of parameters for the correcting device of an electrometric amplifier of 
a mass-spectrometer and provided a method for precise calculations of these parameters. 


In the paper entitled “Some dynamic problems of automatic ac bridges and compensators" V, Yu, Kneller 
dealt mainly with problems of the effect of coupling between the balancing circuits on the dynamic properties of 
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the system, The results obtained hold for a wide range of bridges and compensators and provide a rational choice 
of parameters and adjustments. 


In his paper entitled, “Measurement of continuous quantities with digital registration", S. P, Adamovich 
examines the problem of selecting the measuring frequency for a given accuracy and certain properties of the 
measuring instrument, The result obtained may be used for designing discrete measuring systems. 


In her paper entitled "Hall's effect and the galvanomagnetic effect and their technical application", L. V. 
Sentyurina reviewed the experimental data on the electrical properties of semiconductor materials and evaluated 
the accuracy of transducers based on Hall's effect. The author outlines the possibilities of using Hall's effect for 
measuring magnetic fields, current and power, for detecting and modulating ac signals, etc, 


E, P, Svitina's report entitled "A multiplier based on Hall's effect", dealt with preliminary data of a voltage 
multiplier proposed for calculating devices. 


III, Computer Section 


In his report entitled "A functional transducer with several outputs", M. V. Rybashov provides a criterion of 
the possibility of producing functions of two or more variables provided analytically by means of diode functional 
converters based on the principle of approximating given curves by piecewise linear ones. 


The paper by A, A. Maslov and Yu. G. Purlov on "A functional converter based on the principle of square- 
law approximation” dealt with establishing an analytical and graphical method of approximating functions by 
stretches of parabolas and with a concrete circuit of a converter of that type. 


Zh, A. Kimmel* dealt in his paper with the possibility of using a controlled delay line as a universal analog 
link and with the evaluation of the errors of such a device. 


In his report entitled "Memory units for mathematical analog computers", K. B. Norkin suggested a circuit 
for a memory unit based on capacitors with a long memory time. 


__ In his paper entitled "Devices for recognizing the movement of a quantized scale", A. D, Talantsev examined 
a new circuit for determining the direction of scale movement obtained on the basis of applying the theory of 
logical circuits, 


In his paper entitled *A device for plotting automatically root hodographs and determining the boundaries 
of automatic control system.s* regions of stability", N. N, Mikhailov proposed a new design for an instrument which 
would automatically plot conformal images by means of polynomials. Such an instrument provides root hodo- 
graphs for plotting stability region boundaries for automatic control systems with respect to one and two parameters. 


In his paper entitled “Determination of the accuracy in evaluating variance at the output of automatic moni- 
toring systems when the variance is calculated by means o1 a canonical representation of random functions", L. A. 
Tel'ksnis derived a general formula for the evaluation of the error in the value of the variance, when the random 
function is represented by a finite series, 


In his paper entitled “Digital analogs", A. S. Shileiko examined the possibilities and principles of designing 
digital differential analyzers. 


In his paper entitled “Solution of algebraic and transcendental equations by means of analogs", M. V.Rybashov 
dealt with the development of the algorithm for transforming the initial equations to a form ensuring structural 
stability, 


In his paper entitled * Certain methods of designing elements of digital computers with transistors” B. N. 
Strelkoy examined a number of transistor circuits for elements of digital computers. 


IV. Components and Devices in Automation and Remote control System 


V. M. Dvoretskii dealt in his paper with devices for obtaining small stable flows of liquids in hydro-auto- 
matic systems, This device provices a flow of mineral oil one-tenth for the existing types of throttling devices at 
the same supply pressures, An arrangement was suggested, which made the use of this device possible as a throttling 
element for an auxiliary equpment of a hydraulic governor. 


i- 


T. K. Berends described components for a discrete action computer which were developed by 2 group of 
scientists of the Pneumatic and Hydraulic Automation Laboratory of the Institute of Automation and Remote Control 
of the Academy of Sciences, USSR, and which performed elementary logical operations and time delays, Control 
devices incorporating these elements and suitable for use in industrial automation were demonstrated, 


Chou Ching-leng analyzed in his paper the effect of the variations of the effective area of a membrane and 
their possible skewing on the operation of pneumatic compensating devices, 


V. 1. Chernyshev described pneumatic nonmembrane devices which square and obtain the square root of 


quantities in the form of air pressures, The operating range of the deviceslies between 0-lat. These devices 
can work in industrial automatic control systems, 


A. I, Semikova and V.I. Chernyshev described in their paper systems of new jet pneumatic automatic ele- 
ments (comparison elements, pressure and power amplifiers), On the basis of these elements a theoretically novel 
solution of many problems in pneumatic automation becomes possible. 


I, A, Lyubinskii, V, A. Milyutina and N, V, Pozin in their paper, entitled "A transistor pulse-frequency high 
sensitivity telemetering device", described an instrument designed to measure dc voltages from 0 to 100 my. 


Contrary to existing devices of this type the instrument is very speedy in operation and reliable, 


V. N. Silaev described in his paper various types of pulse generators for remote control devices with non- 
synchronous supplies. The author suggested a circuit of an LC oscillator which provides a highly stable frequency 
by means of high Q tuned circuit and a voltage limiting amplifier in the feedback circuit, 


In his paper entitled "Multi-channel tape recording of programs for continous automatic control systems* 
V. I. Shadrin described the basic principles underlying the various methods of separating channels in multi-channel 
recording and compared them critically, The author showed that the phase-pulse method has all the advantages of 


a space separation of channels, possesses high resistance to interference and can be recommended for industrial ap- 
plication, 


In his paper entitled "Certain problems of designing command devices in digital systems of program control", 


S. P. Khlebnikov described the results of the work in designing a universal command device for controlling milling 
machines, 


In his paper entitled "The problem of calculating dynamic conditions of electromagnetic mechanisms", 
V. B. Gogolevskii described an analytical method of calculation suitable for dc and ac electromagnetic mechanisms, 
The author produced expressions for calculating displacements, speeds, accelerations and traction efforts. 


B. P. Petrkhina reported on the design of an attachment to a dc oscillograph type ENO-1 , which provides 
the measurements of static characteristics of transistors and diodes in common base and common emitter con- 


nections, The basic advantage of the instrument is its simplicity of construction, achieved without lowering the 
quality of reproduction, 


I, V. Prangishvili described a contactless relay consisting of a transistor amplifier in conjunction with a 
magnetic element, The relay has a high operating speed and sensitivity, a large output power, and good stability, 


In his paper entitled "An oscillator with a small second harmonic content", S, F, Stepanov dealt with the 
supply source for a magnetic modulator. 


V. Transmission of Messages Section 


In his paper, V. M, Baikovskii described a discrete difference method of transmitting information based on 
prediction, The author demonstrated the possibility of discrete-difference modulation in telemetering devices and 
described a model of such a device, The results of resistance to interference tests are included and its connection 
with the statistics of the transmitted messages established, 


In his paper, Yu. I. Chugin examined,on the basis of correlation methods, the resistance to interference of 
a single channel, frequency modulated telemetering system in the presence of strong fluctuating noises. The author 
determined mean-square and systematic measurement errors at the output of the frequency selector and showed 
that with a smaller signal to noise ratio the latter error predominates, The author proposed a method of compen- 
sating for the systematic error by means of a minor elaboration of the receiving device circuit, 
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S. N. Diligenskii's paper dealt with differentiation of slow changing signals on the background of noise. For 
a better mean evaluation of the derivative of a signal, which is a smooth function and mixed with noise, the author 
uses a system with a finite memory. 


On the basis of a general criterion of the Bayes type, M. Yu, Gadzhiev provided in his paper a solution of the 
problem of discovering a batch of radar signals in the presence of white noise. Various types of optimum reception 
circuits were examined both for independently fluctuating signals and for the case of a statistical interconnection 
between them, 


In his paper, L. B. Venchkovskii described a photographic method of investigating one and two-dimensional 
probability distribution laws of random process amplitudes, The method is based on photographing the image produced 
on the screen of a cathode-ray oscillograph when its input (or both inputs) is fed with the signal under test and sub- 
sequent examining of the negative by the means of a microphotometer, This method is a development of the 
photometric method and possesses several advantages: simplicity, high-speed analysis, and its documentary nature, 


The report of V. D, Zenkin, S. 1. Ogandzhanyants and N. V. Pozin dealt with the description of a model of 
a pulse~code modulated system, which illustrated the possibility of an automatic exchange of accuracy for a higher 
speed of transmission in the course of operation of a telemetering system, The authors based their development 
on the consideration that in many cases as, for instance, in power networks, it is possible to reduce accuracy when 
rapid changes of one of the parameters are transmitted, which improves the transmitting efficiency of the channel. 
Development experience has shown that the additional complication due to pulse-code modulation is not great. 


VI. Theory of Relay Operated Devices Section 
The section discussed problems of the analysis and synthesis of relay systems. 


In his paper V. P. Didenko dealt with the synthesis of noniterated circuits by a method which is simpler than 
the existing Trakhtenbrot'’s method,and which,as it would appear, can be extended to quasi-noniterated circuits. 
The suggested method is convenient for programming by means of computers, 


O, P, Kuznetsov proposed synthesizing multi-contact circuits by means of a new method in which the circuit 
Operating conditions are registered by a system of equations excluding predicates instead of being registered similarly 
to existing methods on a table of connections, For circuits which require a large number of connection tables, this 
method in many instances considerably simplifies the synthetic process, 


A. D, Talantsev examined in his paper the analysis and synthesis of pulse potential circuits and suggested the 
use of transition operators with Boolean functions which thus represent more adequately the behavior of these cir- 
cuits, The calculus developed by the author produces interesting analogies with a mathematical analysis. 


P, P, Parkhomenko described the basic principles of a machine for automatic checking of the quality of 
production in the electrical industry. The principles suggested by him are suitable for machines checking com- 
ponents (resistances, diodes, etc.) and for checking complicated units (for instance relay units, radio receivers), 


V. M. Kharlamov proposed in his paper a structural classification of telemetering signal types. The author 
developed two multi-pulse systems with combined characteristics for remote control of concentrated objects carried 
out with protective non-operation from any single static error and with a large traffic capacity and he investigated 
the application of the discrete time method in remote control, This method improves the utilization of the time 
characteristic by a factor of 1.5-2. 


B. 1, Timofeev's paper dealt with the mechanization of the finding of minimum-number relay contact cir- 
cuits in the class of disjunctive normal alternatives. The author reviewed the existing method of minimization, 
evaluated them from the point of view of their mechanical applicability, and reported on the schematic circuit 
of a machine he developed for 6 variables according to M.A. Gavrilov's method. 


In his paper, V. D. Kazakoy described the technique of finding the aspect of a function of N variables with 
a maximum number of implicants and gives an approximate evaluation of that maximum for an arbitrary n, This 
result is among other applications useful for calculating the size of the minimization machine memory, 
VIl, Automatic Drives Section 


O, A. Kossov examined the existing circuits of crystal-megnetic amplifiers and suggested a new circuit, He 


examined reversible and non-reversible crystal-magnetic amplifiers with dc and ac outputs and the operating con- 
ditions of transistors in such amplifiers, 


A. G. Cherkashina dealt in her paper with the design of a transistor four stage 50 cps amplifier connected to 
an induction motor type DA-1, 


A. S, Volodin made two reports, In the first one he described a two stage push-pull amplifier with transistor 
elements and a power up to 1.1 kw. In the second one he dealt with a device for an experimental regulator, The 
author showed the possibility of using a magnetic modulator in conjunction with a magnetic amplifier supplying 
an integrating two-phase induction motor for the determination of extremums of functions. 


M. A, Boyarchenkov showed in iis paper the possibility of making a high-speed, cheap and economic system 
for controlling the speed of a dc motor with a magnetic amplifier. The mechanical characteristics of such a drive 
were analyzed and methods for increasing their rigidity indicated, 


Chao Chou-lun's paper dealt with the determination of the technical conditions formaking a controlled 
high-speed flying cutter working in high-speed rolling. The author noted the peculiarities of controlling such an 
object and outlined a method of practical implementation of his suggestion by means of an electromagnetic clutch, 


V. G, Volik's paper dealt with test results of an induction motor with Schenfer's rotor and choke-coil control, 
which allowed one to choose an appropriate power for the motor according to the control range, 


A. A, Yanshin's paper dealt with a system of starting and controlling the speed of an induction motor with a 
short circuited rotor by means of a single-armature frequency converter, analyzed statistical problems of control 
and quoted experimental results obtained. 


S. P, Adamovich, V. M, Baikovskii, V. Yu, Kneller, 
O, P. Kuznetsov, A. A, Maslov, D, P, Petlin, B, P, Pet~ 
rukhin and A, B, Shubin 
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